Abstract

Piezoelectric ceramics are used more in the active control of noise and vibration with sensors and actuators. As a result, a one-dimensional homogeneous thermoelastic fixed piezoelectric rod subjected to a moving heat flow has been considered. The heating rod is affected by three fields, namely, thermal, mechanical, and electric potential. In order to design a reliable smart structure, this study is necessary and the thermopiezoelectric behavior of piezoelectric ceramics must be understood clearly. The Laplace transform method is applied to acquire the distributions of stress, deformation, and temperature fields. The effect of heat source velocity and phase lag parameters on the considered physical fields is displayed graphically, and the results are compared with other models of thermoelasticity. The present findings could find possible applications for designing receiving portions of transformers of Rosen type for voltage magnification.

1. Introduction

Normally, in recent decades, more consideration has been given to the thermoelasticity field due to its applied aspects in various fields, especially geology, geophysics, structures, biology, physics, acoustics, plasma, and so forth. The generalized thermoelasticity theory is one of the improved forms of the classical thermoelasticity and has been modified to remove the defects of the classical coupled thermoelasticity. The modified thermoelasticity models were introduced by Lord and Shulman [1], Green and Lindsay [2], and Green and Naghdi [3] as well as by Tzou [4], Tzou [5], and Tzou [6]. On the other hand, many works, in literature, have applied generalized thermoelasticity theories to investigate the thermoelastic response of solid materials [712].

Piezoelectric materials are an important kind of smart material, as they are able to cause strain or generate electrical energy. Piezoelectric materials are essential materials that have the ability to convert electrical energy into mechanical energy to stimulate strain and act as an engine when these materials are powered by an electric field. Mostly, piezoelectric devices operate in a thermoelastic coupling medium. Often, piezoelectric devices and piezoelectric materials operate in a medium with pyroelectric and thermoelastic coupling.

To illustrate the thermopiezoelectric interaction, various works have been carried out [1325]. Abo-Dahab et al. [26] discussed the electromagnetothermoelastic medium with diffusion and voids in the context of Lord-Shulman or dual-phase-lag models under influence of rotation and gravity.

A reverse problem of including interfaces in a piezoelectric structure was proposed by Nanthakumar et al. [27] in order to overcome this algorithm. The material interfaces are implied by the level sets defined by regularization with total variation penalty words. The reverse issue is iteratively solved, and, for the study of each iteration, the extended finite element approach is used. For three-dimensional structures, the wording is presented, and inclusions of different materials are recognized by multiple level sets.

In order to mathematically model various phenomena in science and engineering, partial differential equations (PDEs) are important. To solve these problems is a key step in achieving accurate understanding of the conduct of natural and engineered structures. In general, analytical methods are typically not sufficient to solve PDEs representing real systems to an appropriate degree. One must use methods of choice. The most popular alternative is likely the finite element approach for engineering problems (FEM). But strong alternatives are also available, such as mesh-free methods and the Isogeometric Analysis (IGA). The basic idea is to approximate the PDE solution with specially constructed functions with desired properties. Samaniego et al. [28] have explored an alternative for approximation Deep Neural Networks (DNNs). In that work, the emphasis is on mechanical problems, and the PDE energy format is analyzed. For a computer learning to solve a mechanical problem, the energy of a mechanical system seems to be a normal loss function.

Hamdia et al. [29] have used sensitivity analysis to define crucial parameters affecting flexoelectric materialsʼ energy conversion factor (ECF). The governing flexoelectricity equations are modeled on a NURBS-based IGA formula that uses their higher-order continuity and thus avoids a complex combined formulation. The input parameters examined include model and material characteristics, and sampling has taken place in the likelihood region using the Latin hypercube sampling (LHS) method. Khalil et al. [30] investigated electromagnetic field and initial stress on a photothermal semiconducting voids medium under thermoelasticity theories. Abo-Dahab et al. [3337] studied the effect of relaxation times in different mediums with several variables. Alotaibi et al. [38] used fractional calculus of thermoelastic P-Waves reflection under influence of gravity and electromagnetic fields. Abd-Alla et al. [39] illustrated the effect of several fields on a generalized thermoelastic medium with voids in the context of Lord-Shulman or dual-phase-lag models.

The Laplace transform method is simple to describe. Given an initial value problem (IVP), apply the Laplace transform operator to both sides of the differential equation. This will transform the differential equation into an algebraic equation whose unknown, F (s), is the Laplace transform of the desired solution. Once you solve this algebraic equation for F (s), take the inverse Laplace transform of both sides; the result is the solution to the original IVP.

In addition, the method used in this research is conceptually simple and easy to program. This method is effective in meeting the requirements of automatic digital computation in an efficient manner for two reasons. First, it is fast (economical) on digital computers now available in contrast to other more complicated methods of numerical quadrature. Second, the method is conceptually simple and requires a minimum of programming effort because the resultant inverse function is given as a Fourier cosine series with the coefficients being appropriate values of the forward transform. Therefore, one avoids the use of involved algorithms which are encountered in methods using orthogonal functions.

In the current investigation, the response of a thermoelastic piezoelectric fixed rod is investigated. Also, the rod exposed to a heat source moves over its axis. The problem is expressed in the context of Tzou [5] model. By using the Laplace transform and the Laplace numerical inversion, we solved the problem. The influence of speed of the moving heat source and the phase lags on displacement, stress, and temperature is explained.

2. Basic Equations

The governing equations of thermopiezoelectric materials without body forces and free charge are as follows [16, 21]:Strain-displacement relations:Stress-strain-temperature and electric field relations:The motion equation:The heat equation proposed by Tzou [6]:Electric and Gauss equations:where represents the stress tensor components, represents the strain tensor components, represents the electric field vector components, represents the thermal conductivity tensor, represents the elastic constants, represents the thermal moduli, represents the function of electric potential, represents the piezoelectric moduli, represents the temperature increment, represents the initial temperature, represents the electric displacement components, represents the displacement vector components, represents the mass density, and represents the electric permittivity. Also, refers to the specific heat, refers to the applied heat source, refers to the pyroelectric constants, refers to the phase lag of heat flux, and refers to the dual-phase-lag parameter. Equations (1)–(5) are the field equations of the generalized linear thermoelastic-piezoelectric applicable for the following:

(i)Dynamical coupled thermoelastic-piezoelectric (CD theory), when (ii)Generalized piezothermoelasticity (LS), when , where is the relaxation time(iii)Generalized linear piezothermoelasticity without energy dissipation (GN model), when and (iv)Generalized piezothermoelasticity with two different phase lags (DPL model), when

3. Problem Formulation

We consider the direction of polarization of a piezoelectric rod finite length that is parallel to the axial direction of the rod. Also, the two ends of rod are assumed to be fixed and thermally insulated. Initially, the rod is not deformed and at rest. The -axis coincides with the axial of the rod. The temperature and the displacement component , for a one-dimensional problem, can be written as

The strain component is given by

We consider that there are no free charges within the piezoelectric rod. Then, Gauss’s law giveswhich can be rewritten as

This leads to . For simplification, we will consider D as constant along the piezoelectric rod. Equations (2), (3), and (5) may be reduced towhere is the thermal expansion coefficient, is the thermal conductivity, and is the piezoelectric constant. To obtain the dimensionless form of the basic equations, we introduce the following variables:

Introducing equation (14) into (11)–(13), we get (dropping the asterisks)where

The nondimensional form of the moving heat source can be taken aswhere is the delta function, is the heat source strength, and is the velocity of along the positive direction of the -axis.

4. Initial and Boundary Conditions

We consider the initial conditions to be

As the rod is fixed and thermally insulated, the boundary conditions are as follows:

5. Laplace Transform

The Laplace transform of any function is given by

After using the Laplace transform, we obtainwhere

The boundary conditions (19) and (20) are also transformed to

6. Solution in the Transformed Field

Eliminating between (23) and (24), we have

The characteristic equation of (28) takes the form

The roots and of equation (29) satisfy the relationswhere and are given by

Solving equation (28), we getwhere all are parameters and

Likewise, the function satisfies the following equation:

The solution to is given from (34) aswhere are some parameters that can be determined from the boundary conditions.

Substituting from equation (32) and from equation (35) into (23), we can find the following relationship:

In order to determine the parameters , we need to apply boundary conditions (26) and (27).

Solving equations (40)–(43), can be expressed as

Substituting and into equation (22), we obtain

7. Laplace Transform Inversion

We will apply a numerical technique to calculate the inversion of Laplace transform to get the physical fields of the problem in the actual field. To reverse the Laplace transform of various fields, we espouse a method based on the expansion of the Fourier series [31]. In this technique, any field in the inverse field of the Laplace transform is reversed in the physical domain using the relationwhere is a sufficiently large integer which was chosen such thatwhere is a small positive and the parameter is positive-free as stated by Honig [31].

8. Numerical Results and Discussion

To illustrate the above results, the material constants of the piezoelectric rod were taken as follows given in Table 1 [8].

The numerical calculation and the discussion of the problem are carried out in two different groups for different positions of and time . The first group is studying how the variations of displacement, temperature, and stress differ with various velocities of the heat source when the other parameters are constants. The second case shows the study of variations of different fields in models of thermoelasticity and the numerical results were also compared in different models. The last case is the study of the dependence of different physical distributions on the time when the parameter remains constant. The figures display that there are two separate phases in the history of electrical potential: before and after heat source leaves the rod. The effect of a heat source prevails in the first stage, while the relation between electrical power and displacement plays an important role in the second stage and contributes to the vibration of the electrical potential as shown in Figure 1.

Figures 24 display the variants of temperature, displacement, and pressure with various values of the parameter in which we detect the following.

Figure 2 displays that when the parameter increases, there is an increase in the profile of temperature.

The distribution of temperature begins at its greatest value at and gradually decreases until the zero wave reaches the direction of propagation of the heat wave.

As Figure 2 shows, although still the wavefront of the simplified solution travels to the right end of the pole, the wave front is on its way back to the left end of the rod. This means that the thermal wave speed is greater than the simplified solution in the current solution. More importantly, while in [32] the authors believed thermal wave fronts were present, they did not actually exist for a longer rod with L = 10, as shown in [32].

From Figure 3, we can observe that the distribution of displacement decreases with increasing the parameter . This results in a reduction in the density of thermal energy along the length of the rod.

It can also be found from Figure 3 that the displacement at is kept at a zero value, which corresponds to the boundary conditions of the problem where the piezoelectric rod is fixed.

The displacement variation increases monotonically and it reaches a peak value at another wave front.

From Figure 4, it is clear to us that the absolute stress decreases with increasing speed due to the decrease in the energy of the applied heat source along the axial of the rod.

The distributions of physical fields are graphically depicted in Figures 57 for different theories of thermal elasticity. This group illustrates the influence of two-phase-lag parameters on the numerical results, and we observed a slight difference in the value of the distributions in the three models.

Figures 810 are presented to illustrate the distributions of the physical fields at various values of time to show the influence of time in each field. We note the following:

The temperature profile increases with increasing value of time during certain periods.

The maximum values of the distributions of displacement and stress increase as time increases.

Finally, through Figures 210, it becomes clear that once the immediate time is given, the nonzero values of temperature, displacement, and thermal stress are only in a limited region and, beyond this area, gradually fade away. This explains why the heat dissipation speed is limited in the piezoelectric rod, as this is completely different as compared to the classical thermal models, which predict unlimited velocity. Due to the limited heat diffusion, the heat-disturbed region is limited when the appropriate time is given, resulting in the pressure and displacement caused by the heat also in a specific region.

From Figure 3, we can observe that the distribution of displacement decreases with increasing the parameter . This results in a reduction in the density of thermal energy along the length of the rod.

It can also be found from Figure 3 that the displacement at is kept at a zero value which corresponds to the boundary conditions of the problem where the piezoelectric rod is fixed.

The displacement variation increases monotonically and it reaches a peak value at another wave front.

From Figure 4, it is clear to us that the absolute stress decreases with increasing speed due to the decrease in the energy of the applied heat source along the axial of the rod.

The temperature field stabilizes as the heat source moves away from the rod, contrary to the acceleration, electrical potential, and stress fields. There can, however, be some variations before the temperature achieves a constant value, depending on the thermal rest time.

The distributions of physical fields are graphically depicted in Figures 57 for different theories of thermal elasticity. This group illustrates the influence of two-phase-lag parameters on the numerical results, and we observed a slight difference in the value of the distributions in the three models. The propagation of thermal waves is very sensitive to the speed of heat sources. Higher speed allows the thermal disturbance to spread more rapidly. The magnitudes of the nondimensional variables increase as the speed of the moving heat resource increases. As the largest values of curves indicate obviously, the effect of the thermal source transmission speed on all considered quantities is quite high. The influence of the heat source is an inherent consideration not to be overlooked in the evaluation of heat stress and fields at the heat source.

Figures 810 are presented to illustrate the distributions of the physical fields at various values of time to show the influence of time in each field. We note the following:

The temperature profile increases with increasing value of time during certain periods.

The maximum values of the distributions of displacement and stress increase as time increases.

Finally, through Figures 210, it becomes clear that once the immediate time is given, the nonzero values of temperature, displacement, and thermal stress are only in a limited region and, beyond this area, gradually fade away. This explains why the heat dissipation speed is limited in the piezoelectric rod, as this is completely different as compared to the classical thermal models, which predict unlimited velocity. Due to the limited heat diffusion, the heat-disturbed region is limited when the appropriate time is given, resulting in the pressure and displacement caused by the heat also in a specific region.

9. Conclusion

The response of a thermoelastic piezoelectric rod exposed to a heat source is studied using the model of thermoelasticity with phase lags. The coupled governing equations of piezoelectricity and thermoelasticity are expressed. Based on the Laplace transform and an approximation inversion technique, the system of equations is solved. The results of the physical variables are displayed graphically and discussed.

From the graphs and theoretical discussion, we can conclude the following: With a comparison between the different models of thermoelasticity, that is, coupled theory (CD), Lord and Shulman (LS), and Tzou (DPL) models, there are differences in the results between different models.

The phenomenon of limited propagation speeds is manifested in all these forms. The heat wave interface is moving forward at a limited speed over time. This indicates that the mechanism of heat wave propagation in the case of generalized theories of thermoelasticity differs from that in the classical theory.

The values of temperature, exclusion, and nondimensional stress increase with the increase of the moving heat source. The value of has a significant impact on all areas of distribution in the case of the three theories.

Data Availability

The data are supported from the authors’ published papers, and all the results are new.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. D-511-352-1441. The authors, therefore, gratefully acknowledge DSR for technical and financial support.