Magneto-Electric Effect on Guided Waves in Functionally Graded Piezoelectric–Piezomagnetic Fan-Shaped Cylindrical Structures

Functionally graded piezoelectric–piezomagnetic (FGPP) material simultaneously consists of piezomagnetic and piezoelectric phases, which are able to convert energy among mechanical, electric, and magnetic fields. The magneto-electric effect on waves in FGPP fan-shaped cylindrical structures is studied by exploiting the double Legendre orthogonal polynomial method. By means of the Heaviside function, the initial conditions are brought into wave motion equations. Dispersion properties, electric and magnetic potential, and the Poynting vector are calculated. Subsequently, the effect of the graded variation and geometric size on wave characteristics is analyzed. The FGPP fan-shaped cylindrical structures are of complex geometrical shape and material inhomogeneity, so their influences on the magneto-electric effect are the focus of discussion. Results reveal that the cut-off frequencies have a negative relationship with the cross-section area of the structure. The magneto-electric effect could be adjusted via altering the geometric size of the cross-section. These results can be utilized to design and optimize piezoelectric–piezomagnetic fan-shaped transducers.


Introduction
It is well known that the piezoelectric material has the piezoelectric effect and converse piezoelectric effect, i.e., energy transformation between mechanical and electric fields. Therefore, it has been widely used in various fields, such as semiconductors [1], transducers [2], actuators [3], micromotors [4], and piezoelectric charge-coupled devices [5]. These devices' performance could be optimized by changing the composition of piezoelectric material [6]. With the development of material science, a class of composite material, the functionally graded piezoelectric-piezomagnetic (FGPP) material, simultaneously consisting of piezomagnetic and piezoelectric phases, was developed, which exhibited a coupled mechanical, magnetic, and electric field. It is able to convert energy among mechanical, electric, and magnetic fields, which is known as the magneto-electric effect [7]. Accordingly, FGPP has promising applications in different areas, such as transducers, sensors, and storage devices [8,9].
To design and optimize the FGPP transducers, wave characteristics in FGPP structures were investigated by utilizing different kinds of models and methods. The propagative behaviors of guided waves in the FGPP plates were investigated by Cao et al. [10], utilizing the power series technique. SH waves, propagating in multi-layered FGPP structures, were studied by Singh and Rokne [11]. For the FGPP cylindrical structures, the gradient direction is along the r direction. Accordingly, material parameters are compactly expressed as the following form, via the least square method.
Here, the function I(r, θ) [32] is introduced into Equation (1), to automatically satisfy the tractionfree and open-circuit boundary conditions, which can be denoted as follows:  For the FGPP cylindrical structures, the gradient direction is along the r direction. Accordingly, material parameters are compactly expressed as the following form, via the least square method.
where p i m,j (i = 1, 2, 3, 4, 5) denote expansion coefficients and where P m and P j represent the Legendre polynomials with the order m and j, respectively. m and j are evaluated from 0 to ∞ in theory. However, as a matter of fact, the Equation (5) are convergent at some finite values m = M and j = J.
Multiply Equations (4a)-(4e) by Q n (r) × Q p (θ), with n running from 0 to M and p running from 0 to J, respectively. Subsequently, integrating over r from a to b, and over θ from 0 to β. Then, we can obtain the following system by utilizing the orthonormality: where l M n,p,m,j and l A n,p,m,j αγ (α, γ = 1, 2, 3, 4, 5), which can be calculated based on Equation (6), are the elements of the non-symmetric matrices. Equation (7) can be transformed as the following matrix system, and the detailed transformation process is given in the Appendix A.
where the dimensions of matrices l A n,p,m,j αγ (α, γ = 1, 2, 3) and l M n,p,m,j are all (M + 1) × (J + 1). Consequently, Equation (8) is an eigenvalue problem about angular frequency ω. The profiles of displacement components are calculated according to the eigenvectors. At last, the phase velocity is obtained via the equation V ph = ω/k.

Numerical Results
In this section, the equivalent parameters of the FGPP cylindrical structures are calculated utilizing the Voigt-type model [29].
According to the abovementioned equations, the computer programs in the light of the double orthogonal polynomial method are written using the software "Mathematica" to calculate the dispersion curves and the Poynting vectors.
In the present paper, such a kind of FGPP fan-shaped cylindrical structure is considered: the inner layer material is CoFe 2 O 4 , and the outer layer material is BaTiO 3 . Their material constants are listed in the Table 1 in the reference [34].

Comparison with Rectangular Bar
As far as we know, rare references for the two-dimensional FGPP fan-shaped cylindrical structures are available. Consequently, we simplify the model as a purely elastic structure, and assume its geometric size is a = (10 −4 − 1) mm, b = 10 4 mm, and β = 10 −4 rad, respectively, to compare with a square steel bar using the two dimensional Rayleigh-Ritz method [35]. Its cross-section area is 1 mm 2 . Their phase velocity curves are shown in the Figure 2, where the lines are the authors' results, and the dotted lines are results from reference [35]. c p represents the phase velocity, and c s represents the shear velocity. These results of two methods are overlapped completely. Accordingly, the correctness of the present method is confirmed. Equation (7) can be transformed as the following matrix system, and the detailed transformation process is given in the Appendix A.

Numerical Results
In this section, the equivalent parameters of the FGPP cylindrical structures are calculated utilizing the Voigt-type model [29].
According to the abovementioned equations, the computer programs in the light of the double orthogonal polynomial method are written using the software "Mathematica" to calculate the dispersion curves and the Poynting vectors.
In the present paper, such a kind of FGPP fan-shaped cylindrical structure is considered: the inner layer material is CoFe2O4, and the outer layer material is BaTiO3. Their material constants are listed in the Table 1 in the reference [34].

Comparison with Rectangular Bar
As far as we know, rare references for the two-dimensional FGPP fan-shaped cylindrical structures are available. Consequently, we simplify the model as a purely elastic structure, and assume its geometric size is a = (10 −4 − 1) mm, b = 10 4 mm, and β = 10 −4 rad, respectively, to compare with a square steel bar using the two dimensional Rayleigh-Ritz method [35]. Its cross-section area is 1 mm 2 . Their phase velocity curves are shown in the Figure 2, where the lines are the authors' results, and the dotted lines are results from reference [35]. cp represents the phase velocity, and cs represents the shear velocity. These results of two methods are overlapped completely. Accordingly, the correctness of the present method is confirmed. Subsequently, a comparison with an FGPP square bar from reference [24] is made, and the influence of the radius-thickness ratio is analyzed. To ensure the area of cross-section has very little difference, their geometric parameters are η = 10 and β = 1/9.5 rad, η = 100 and β = 1/99.5 rad, and η = 1000 and β = 1/999.5 rad, respectively. Their thicknesses are both 1 mm. Their phase velocity curves Subsequently, a comparison with an FGPP square bar from reference [24] is made, and the influence of the radius-thickness ratio is analyzed. To ensure the area of cross-section has very little difference, their geometric parameters are η = 10 and β = 1/9.5 rad, η = 100 and β = 1/99.5 rad, and η = 1000 and β = 1/999.5 rad, respectively. Their thicknesses are both 1 mm. Their phase velocity curves for the fourth mode are illustrated in the Figure 3. The cases of the other modes are similar. Vp is a dimensionless parameter (Vp = (c p *d)/c s1 , where c p is the phase velocity value, and c s1 is the shear velocity value of CoFe 2 O 4 ). Here, Figure 3b is the partial enlarged drawing of Figure 3a. We can note from these figures that the phase velocity for the fan-shaped cylindrical structures is getting closer to that of rectangular bar as the radius-thickness ratio increases. Moreover, the dotted line and the line are overlapped completely as the radius-thickness ratio η = 1000. for the fourth mode are illustrated in the Figure 3. The cases of the other modes are similar. Vp is a dimensionless parameter (Vp = (cp*d)/cs1, where cp is the phase velocity value, and cs1 is the shear velocity value of CoFe2O4). Here, Figure 3b is the partial enlarged drawing of Figure 3a. We can note from these figures that the phase velocity for the fan-shaped cylindrical structures is getting closer to that of rectangular bar as the radius-thickness ratio increases. Moreover, the dotted line and the line are overlapped completely as the radius-thickness ratio η = 1000.

Convergence Confirmation
To confirm the convergence of the present method, the phase velocities of a linear FGPP fanshaped cylindrical structure (β = π/6 and η = 2) with various M and J are calculated. Their phase velocity values are listed in Table 1 at kd = 1.01.
Here, as M or J varies, the relative error is It is convergent in the present paper. c represents the phase velocity. We can note that the first two modes are convergent as M = J = 5, and the first three modes are convergent as M = 6 and J = 5.

Convergence Confirmation
To confirm the convergence of the present method, the phase velocities of a linear FGPP fan-shaped cylindrical structure (β = π/6 and η = 2) with various M and J are calculated. Their phase velocity values are listed in Table 1 at kd = 1.01.
Here, as M or J varies, the relative error is It is convergent in the present paper. c represents the phase velocity. We can note that the first two modes are convergent as M = J = 5, and the first three modes are convergent as M = 6 and J = 5.

The Magneto-Electric Effect
The phase velocity curves of FGPP fan-shaped cylindrical structure, corresponding graded piezoelectric structure, and graded piezomagnetic structure, are shown in the Figure 4. Figure 4b is the partially enlarged figure of Figure 4a. They have the same geometry: η = 2 and β = π/6. The phase velocity values of the FGPP structure and piezoelectric structure are almost the same. We can only note that the phase velocity for FGPP structure is little higher than that of the piezoelectric structure in . Moreover, the phase velocity for FGPP structures is higher than that of the piezomagnetic structure. Therefore, the influence of piezoelectric effect is much more significant than that of the piezomagnetic effect, which is similar to other FGPP structures [33]. Figure 5 shows the phase velocity curves of the graded elastic and piezomagnetic structure. We can note that their differences are also very little. The piezomagnetic effect makes the phase velocity decrease only a little.
The phase velocity curves of FGPP fan-shaped cylindrical structure, corresponding graded piezoelectric structure, and graded piezomagnetic structure, are shown in the Figure 4. Figure 4b is the partially enlarged figure of Figure 4a. They have the same geometry: η = 2 and β = π/6. The phase velocity values of the FGPP structure and piezoelectric structure are almost the same. We can only note that the phase velocity for FGPP structure is little higher than that of the piezoelectric structure in Figure 4b. Moreover, the phase velocity for FGPP structures is higher than that of the piezomagnetic structure. Therefore, the influence of piezoelectric effect is much more significant than that of the piezomagnetic effect, which is similar to other FGPP structures [33]. Figure 5 shows the phase velocity curves of the graded elastic and piezomagnetic structure. We can note that their differences are also very little. The piezomagnetic effect makes the phase velocity decrease only a little. Subsequently, the electric and magnetic potential distributions are calculated and illustrated in Figures 6 and 7, at kd = 2.01 and kd = 120.01, respectively. Here, we just calculated the absolute values of the electric and magnetic potential. Amplitudes of electric potential are much bigger (about 10 3 times) than that of the magnetic potential, which has a relationship with the electric and magnetic parameters. The phenomenon that the piezoelectric effect is stronger than the piezomagnetic effect is confirmed again. Furthermore, regardless of the amplitude, there are obvious differences between the distribution shapes at small wavenumbers, but they are similar at bigger wavenumbers. The phase velocity curves of FGPP fan-shaped cylindrical structure, corresponding graded piezoelectric structure, and graded piezomagnetic structure, are shown in the Figure 4. Figure 4b is the partially enlarged figure of Figure 4a. They have the same geometry: η = 2 and β = π/6. The phase velocity values of the FGPP structure and piezoelectric structure are almost the same. We can only note that the phase velocity for FGPP structure is little higher than that of the piezoelectric structure in Figure 4b. Moreover, the phase velocity for FGPP structures is higher than that of the piezomagnetic structure. Therefore, the influence of piezoelectric effect is much more significant than that of the piezomagnetic effect, which is similar to other FGPP structures [33]. Figure 5 shows the phase velocity curves of the graded elastic and piezomagnetic structure. We can note that their differences are also very little. The piezomagnetic effect makes the phase velocity decrease only a little. Subsequently, the electric and magnetic potential distributions are calculated and illustrated in Figures 6 and 7, at kd = 2.01 and kd = 120.01, respectively. Here, we just calculated the absolute values of the electric and magnetic potential. Amplitudes of electric potential are much bigger (about 10 3 times) than that of the magnetic potential, which has a relationship with the electric and magnetic parameters. The phenomenon that the piezoelectric effect is stronger than the piezomagnetic effect is confirmed again. Furthermore, regardless of the amplitude, there are obvious differences between the distribution shapes at small wavenumbers, but they are similar at bigger wavenumbers. Subsequently, the electric and magnetic potential distributions are calculated and illustrated in Figures 6 and 7, at kd = 2.01 and kd = 120.01, respectively. Here, we just calculated the absolute values of the electric and magnetic potential. Amplitudes of electric potential are much bigger (about 10 3 times) than that of the magnetic potential, which has a relationship with the electric and magnetic parameters. The phenomenon that the piezoelectric effect is stronger than the piezomagnetic effect is confirmed again. Furthermore, regardless of the amplitude, there are obvious differences between the distribution shapes at small wavenumbers, but they are similar at bigger wavenumbers.  The abovementioned phenomenon can be detailed according to Equation (4) and Table 1 in the reference [34]. We can note from Equation (4) that the piezoelectric parameters and piezomagnetic parameters have the same influence. The piezoelectric and piezomagnetic parameters have a positive relationship with the magneto-electric effect, and the dielectric and magnetic permeability coefficients are negatively related to the magneto-electric effect. Besides, the average absolute values of piezoelectric parameters are about 50 times the piezomagnetic parameters, and the average values of magnetic permeability coefficients are about 5 × 10 6 times the dielectric coefficients. Hence, the influence of piezomagnetic is very weak, and the amplitude of electric potential is about 10 3 times of magnetic potential, due to the much higher magnetic permeability coefficients.

The Influence of the Graded Functions
For the FGPP material, the graded function, a significant index in the material design, has significant influence on the material performance. Consequently, it also has remarkable influence on the wave characteristics. To investigate the influence of the graded function, two FGPP cylindrical structures with power series gradient functions are considered, i.e., V1(r) = (r − a)/d and V1(r) = [(r − a)/d] 3 . They have the same geometry: η = 2, β = π/6. Figure 8 shows variation curves of material properties with different graded functions. Figure 9 illustrates their dispersion curves of the fourth to sixth modes. For a given mode, phase velocities with cubic function are higher than those with a linear function, since it is composed of more CoFe2O4, and the wave velocity of BaTiO3 is lower than that of CoFe2O4. The abovementioned phenomenon can be detailed according to Equation (4) and Table 1 in the reference [34]. We can note from Equation (4) that the piezoelectric parameters and piezomagnetic parameters have the same influence. The piezoelectric and piezomagnetic parameters have a positive relationship with the magneto-electric effect, and the dielectric and magnetic permeability coefficients are negatively related to the magneto-electric effect. Besides, the average absolute values of piezoelectric parameters are about 50 times the piezomagnetic parameters, and the average values of magnetic permeability coefficients are about 5 × 10 6 times the dielectric coefficients. Hence, the influence of piezomagnetic is very weak, and the amplitude of electric potential is about 10 3 times of magnetic potential, due to the much higher magnetic permeability coefficients.

The Influence of the Graded Functions
For the FGPP material, the graded function, a significant index in the material design, has significant influence on the material performance. Consequently, it also has remarkable influence on the wave characteristics. To investigate the influence of the graded function, two FGPP cylindrical structures with power series gradient functions are considered, i.e., V 1 (r) = (r − a)/d and V 1 (r) = [(r − a)/d] 3 . They have the same geometry: η = 2, β = π/6. Figure 8 shows variation curves of material properties with different graded functions. Figure 9 illustrates their dispersion curves of the fourth to sixth modes. For a given mode, phase velocities with cubic function are higher than those with a linear function, since it is composed of more CoFe 2 O 4 , and the wave velocity of BaTiO 3 is lower than that of CoFe 2 O 4 .   Subsequently, the influence of graded function on magneto-electric effect is discussed. Figure 10 shows the corresponding phase velocity curves. Comparing with the Figure 4, we can note that the piezoelectric effect obviously weakens. This is because the volume content for BaTiO 3 decreases as n increases, and the piezoelectric effect mainly results from BaTiO 3 . Subsequently, the influence of graded function on magneto-electric effect is discussed. Figure 10 shows the corresponding phase velocity curves. Comparing with the Figure 4, we can note that the piezoelectric effect obviously weakens. This is because the volume content for BaTiO3 decreases as n increases, and the piezoelectric effect mainly results from BaTiO3.

The Influence of the Geometric Size
Control wave characteristics via altering the geometric size is supposed as an efficient approach, such as the cut-off frequency. Accordingly, for the fan-shaped cross-section, the influence of variation in the angular measures and the radius-thickness ratios are investigated. Firstly, two linearly graded cylindrical structures with different angular measures (β = π/4, π/6) and η = 2 are considered. The dispersion curves are shown in the Figure 11. The first four modes start from 0 frequency. Moreover, for higher modes, the cut-off frequencies increase as angular measures decrease. For a given η (for example η = 2), as the angular measures decrease, the cross-section area decreases, and the cut-off frequencies increase. Therefore, we hold the view that cut-off frequencies have a negative relationship with the cross-sectional area of the structure.  Subsequently, the influence of graded function on magneto-electric effect is discussed. Figure 10 shows the corresponding phase velocity curves. Comparing with the Figure 4, we can note that the piezoelectric effect obviously weakens. This is because the volume content for BaTiO3 decreases as n increases, and the piezoelectric effect mainly results from BaTiO3.

The Influence of the Geometric Size
Control wave characteristics via altering the geometric size is supposed as an efficient approach, such as the cut-off frequency. Accordingly, for the fan-shaped cross-section, the influence of variation in the angular measures and the radius-thickness ratios are investigated. Firstly, two linearly graded cylindrical structures with different angular measures (β = π/4, π/6) and η = 2 are considered. The dispersion curves are shown in the Figure 11. The first four modes start from 0 frequency. Moreover, for higher modes, the cut-off frequencies increase as angular measures decrease. For a given η (for example η = 2), as the angular measures decrease, the cross-section area decreases, and the cut-off frequencies increase. Therefore, we hold the view that cut-off frequencies have a negative relationship with the cross-sectional area of the structure.

The Influence of the Geometric Size
Control wave characteristics via altering the geometric size is supposed as an efficient approach, such as the cut-off frequency. Accordingly, for the fan-shaped cross-section, the influence of variation in the angular measures and the radius-thickness ratios are investigated. Firstly, two linearly graded cylindrical structures with different angular measures (β = π/4, π/6) and η = 2 are considered. The dispersion curves are shown in the Figure 11. The first four modes start from 0 frequency. Moreover, for higher modes, the cut-off frequencies increase as angular measures decrease. For a given η (for example η = 2), as the angular measures decrease, the cross-section area decreases, and the cut-off frequencies increase. Therefore, we hold the view that cut-off frequencies have a negative relationship with the cross-sectional area of the structure. Subsequently, to confirm the above view, three linearly graded cylindrical structures, with the radius-thickness ratios (η = 2, 3, 4) and β = π/6, are also studied. Figure 12 shows their dispersion curves. We can note that the cut-off frequencies decrease as the radius-thickness ratios increase. For the given angular measure (for example β = π/6), as the radius-thickness ratios increase, the crosssection area increases, and the cut-off frequency decreases. 4.5 Figure 11. The dispersion curves with different angular measure: red lines-β = π/4, blue lines-β = π/6. Subsequently, to confirm the above view, three linearly graded cylindrical structures, with the radius-thickness ratios (η = 2, 3, 4) and β = π/6, are also studied. Figure 12 shows their dispersion curves. We can note that the cut-off frequencies decrease as the radius-thickness ratios increase. For the given angular measure (for example β = π/6), as the radius-thickness ratios increase, the cross-section area increases, and the cut-off frequency decreases. Figure 11. The dispersion curves with different angular measure: red lines-β = π/4, blue lines-β = π/6. Subsequently, to confirm the above view, three linearly graded cylindrical structures, with the radius-thickness ratios (η = 2, 3, 4) and β = π/6, are also studied. Figure 12 shows their dispersion curves. We can note that the cut-off frequencies decrease as the radius-thickness ratios increase. For the given angular measure (for example β = π/6), as the radius-thickness ratios increase, the crosssection area increases, and the cut-off frequency decreases. Then, comparing with 1-D cylindrical structures, there may be some different phenomena because of its complex geometrical shape. Therefore, the influence of the geometric size of the crosssection on the magneto-electric effect is also studied. Figure 13 illustrates phase velocity curves for fan-shaped cylindrical structures with η = 3 and β = π/6. Comparing with the Figure 4, a phenomenon is found that the piezoelectric effect becomes weak at the transition section (about 1.3-1.6 MHz-m), where the dispersion becomes strong from weakness. Figures 14 and 15 show the corresponding phase velocity curves for the fan-shaped cylindrical structures with β = π/4 and β = π/8, respectively. It can be seen from these two figures that, for the fan-shaped cylindrical structures with β = π/4, the influence of piezoelectric effect also becomes very weak at the transition section (about 1.6-2.1 MHzm), where the dispersion becomes strong from weakness. In summary, the magneto-electric effect could be adjusted via altering the geometric size.  Then, comparing with 1-D cylindrical structures, there may be some different phenomena because of its complex geometrical shape. Therefore, the influence of the geometric size of the cross-section on the magneto-electric effect is also studied. Figure 13 illustrates phase velocity curves for fan-shaped cylindrical structures with η = 3 and β = π/6. Comparing with the Figure 4, a phenomenon is found that the piezoelectric effect becomes weak at the transition section (about 1.3-1.6 MHz-m), where the dispersion becomes strong from weakness. Figures 14 and 15 show the corresponding phase velocity curves for the fan-shaped cylindrical structures with β = π/4 and β = π/8, respectively. It can be seen from these two figures that, for the fan-shaped cylindrical structures with β = π/4, the influence of piezoelectric effect also becomes very weak at the transition section (about 1.6-2.1 MHz-m), where the dispersion becomes strong from weakness. In summary, the magneto-electric effect could be adjusted via altering the geometric size. Figure 11. The dispersion curves with different angular measure: red lines-β = π/4, blue lines-β = π/6. Subsequently, to confirm the above view, three linearly graded cylindrical structures, with the radius-thickness ratios (η = 2, 3, 4) and β = π/6, are also studied. Figure 12 shows their dispersion curves. We can note that the cut-off frequencies decrease as the radius-thickness ratios increase. For the given angular measure (for example β = π/6), as the radius-thickness ratios increase, the crosssection area increases, and the cut-off frequency decreases. Then, comparing with 1-D cylindrical structures, there may be some different phenomena because of its complex geometrical shape. Therefore, the influence of the geometric size of the crosssection on the magneto-electric effect is also studied. Figure 13 illustrates phase velocity curves for fan-shaped cylindrical structures with η = 3 and β = π/6. Comparing with the Figure 4, a phenomenon is found that the piezoelectric effect becomes weak at the transition section (about 1.3-1.6 MHz-m), where the dispersion becomes strong from weakness. Figures 14 and 15 show the corresponding phase velocity curves for the fan-shaped cylindrical structures with β = π/4 and β = π/8, respectively. It can be seen from these two figures that, for the fan-shaped cylindrical structures with β = π/4, the influence of piezoelectric effect also becomes very weak at the transition section (about 1.6-2.1 MHzm), where the dispersion becomes strong from weakness. In summary, the magneto-electric effect could be adjusted via altering the geometric size.

Waves at High Frequencies
For the 1-D cylindrical structures [25], the phase velocity values of higher modes always increase as the frequency decreases, and are close to the shear velocity cs. However, the case for 2-D fan-shaped cylindrical structures is quite different. A homogeneous (V1(r) = 1) fan-shaped cylindrical structure with η = 2 and α = π/6 is taken into account. The phase velocity curves for the fifth and sixth mode are shown in Figure 16. The velocity values at high frequencies approximatively approach a value which is below cs. This is because wave motions at high frequencies mainly concentrate near corners of the fan-shaped cylindrical structures, while the other places remain almost motionless (see in the Section 3.8). The velocity for waves propagating near the edges at high frequencies is lower than that for waves propagating at a surface. Subsequently, two linear FGPP fan-shaped cylindrical structures with different angular measures are considered. Figure 17 shows the corresponding phase velocity curves at high frequencies. Phase velocity values of the fan-shaped cylindrical structures, with β = π/6, are higher

Waves at High Frequencies
For the 1-D cylindrical structures [25], the phase velocity values of higher modes always increase as the frequency decreases, and are close to the shear velocity c s . However, the case for 2-D fan-shaped cylindrical structures is quite different. A homogeneous (V 1 (r) = 1) fan-shaped cylindrical structure with η = 2 and α = π/6 is taken into account. The phase velocity curves for the fifth and sixth mode are shown in Figure 16. The velocity values at high frequencies approximatively approach a value which is below c s . This is because wave motions at high frequencies mainly concentrate near corners of the fan-shaped cylindrical structures, while the other places remain almost motionless (see in the Section 3.8). The velocity for waves propagating near the edges at high frequencies is lower than that for waves propagating at a surface.

Waves at High Frequencies
For the 1-D cylindrical structures [25], the phase velocity values of higher modes always increase as the frequency decreases, and are close to the shear velocity cs. However, the case for 2-D fan-shaped cylindrical structures is quite different. A homogeneous (V1(r) = 1) fan-shaped cylindrical structure with η = 2 and α = π/6 is taken into account. The phase velocity curves for the fifth and sixth mode are shown in Figure 16. The velocity values at high frequencies approximatively approach a value which is below cs. This is because wave motions at high frequencies mainly concentrate near corners of the fan-shaped cylindrical structures, while the other places remain almost motionless (see in the Section 3.8). The velocity for waves propagating near the edges at high frequencies is lower than that for waves propagating at a surface. Subsequently, two linear FGPP fan-shaped cylindrical structures with different angular measures are considered. Figure 17 shows the corresponding phase velocity curves at high frequencies. Phase velocity values of the fan-shaped cylindrical structures, with β = π/6, are higher Subsequently, two linear FGPP fan-shaped cylindrical structures with different angular measures are considered. Figure 17 shows the corresponding phase velocity curves at high frequencies. Phase velocity values of the fan-shaped cylindrical structures, with β = π/6, are higher than that of structure with β = π/9 and β = π/12, i.e., phase velocity values at high frequencies decrease as with the decrease of β. This feature is similar to that of wedge waves.

The Stress, Electric, and Magnetic Displacement Distribution
The stress, electric displacement, and magnetic displacement of the first mode for a linear FGPP fan-shaped cylindrical structure (η = 2 and β = π/6) at kd = 2.01 are calculated and illustrated in Figure 18. Here, to save the space, only three representative figures are illustrated. It can be easily seen that the initial conditions are satisfied extremely well, i.e., T rr = T rθ = T rz = D r = B r = 0 at r = a and r = b; and T θθ = T rθ = T θz = D θ = B θ = 0 at θ = 0 and θ = β. than that of structure with β = π/9 and β = π/12, i.e., phase velocity values at high frequencies decrease as with the decrease of β. This feature is similar to that of wedge waves.

The Stress, Electric, and Magnetic Displacement Distribution
The stress, electric displacement, and magnetic displacement of the first mode for a linear FGPP fan-shaped cylindrical structure (η = 2 and β = π/6) at kd = 2.01 are calculated and illustrated in Figure  18. Here, to save the space, only three representative figures are illustrated. It can be easily seen that the initial conditions are satisfied extremely well, i.e.,  Figure 18. The stress, electric, and magnetic displacement of the first mode for a linearly FGPP cylindrical structure (η = 2 and β = π/6) at kd = 2.01.

The Poynting Vector
It is well known that the Poynting vector denotes the power flow density, which can be calculated by

The Stress, Electric, and Magnetic Displacement Distribution
The stress, electric displacement, and magnetic displacement of the first mode for a linear FGPP fan-shaped cylindrical structure (η = 2 and β = π/6) at kd = 2.01 are calculated and illustrated in Figure  18. Here, to save the space, only three representative figures are illustrated. It can be easily seen that the initial conditions are satisfied extremely well, i.e.,  Figure 18. The stress, electric, and magnetic displacement of the first mode for a linearly FGPP cylindrical structure (η = 2 and β = π/6) at kd = 2.01.
distributions are mainly concentrated near the outer side with more BaTiO 3 . This is because elasticity modulus of BaTiO 3 is smaller than that of CoFe 2 O 4 . Accordingly, the stiffness of BaTiO 3 is smaller, the displacement is bigger, and the energy mainly transmits in this region. Moreover, the Poynting vector distributions are concentrated near the boundaries, especially near the corner.
For the guided waves propagating in the axial direction, the Poynting vectors in the radial and circumference directions are 0. Figures 19-21 illustrate the Poynting vector distributions of the first two modes for the linear FGPP cylindrical structure (η = 2 and β = π/6), at kd = 2.01 and kd = 120.01, respectively. The detailed distributions at kd = 120.01 in the Figure 20 cannot be clearly seen. Hence, the partial enlarged drawings are made, as shown in the Figure 21. They seem to be discontinuous, owing to the existence of platforms. However, they are continuous, in fact, and truncation of the larger values leads to platforms. Furthermore, for the big wavenumber case, the Poynting vector distributions are mainly concentrated near the outer side with more BaTiO3. This is because elasticity modulus of BaTiO3 is smaller than that of CoFe2O4. Accordingly, the stiffness of BaTiO3 is smaller, the displacement is bigger, and the energy mainly transmits in this region. Moreover, the Poynting vector distributions are concentrated near the boundaries, especially near the corner.    For the guided waves propagating in the axial direction, the Poynting vectors in the radial and circumference directions are 0. Figures 19-21 illustrate the Poynting vector distributions of the first two modes for the linear FGPP cylindrical structure (η = 2 and β = π/6), at kd = 2.01 and kd = 120.01, respectively. The detailed distributions at kd = 120.01 in the Figure 20 cannot be clearly seen. Hence, the partial enlarged drawings are made, as shown in the Figure 21. They seem to be discontinuous, owing to the existence of platforms. However, they are continuous, in fact, and truncation of the larger values leads to platforms. Furthermore, for the big wavenumber case, the Poynting vector distributions are mainly concentrated near the outer side with more BaTiO3. This is because elasticity modulus of BaTiO3 is smaller than that of CoFe2O4. Accordingly, the stiffness of BaTiO3 is smaller, the displacement is bigger, and the energy mainly transmits in this region. Moreover, the Poynting vector distributions are concentrated near the boundaries, especially near the corner.    For the guided waves propagating in the axial direction, the Poynting vectors in the radial and circumference directions are 0. Figures 19-21 illustrate the Poynting vector distributions of the first two modes for the linear FGPP cylindrical structure (η = 2 and β = π/6), at kd = 2.01 and kd = 120.01, respectively. The detailed distributions at kd = 120.01 in the Figure 20 cannot be clearly seen. Hence, the partial enlarged drawings are made, as shown in the Figure 21. They seem to be discontinuous, owing to the existence of platforms. However, they are continuous, in fact, and truncation of the larger values leads to platforms. Furthermore, for the big wavenumber case, the Poynting vector distributions are mainly concentrated near the outer side with more BaTiO3. This is because elasticity modulus of BaTiO3 is smaller than that of CoFe2O4. Accordingly, the stiffness of BaTiO3 is smaller, the displacement is bigger, and the energy mainly transmits in this region. Moreover, the Poynting vector distributions are concentrated near the boundaries, especially near the corner.

Conclusions
According to the 3D linearly magneto-electric-elastic theory, waves in FGPP fan-shaped cylindrical structures are studied via exploiting the double Legendre orthogonal polynomial method. The magneto-electric effect is detailed. Based on the above numerical results, the following conclusions can be obtained: (1) If the radius-thickness ratio is bigger than 1000, fan-shaped cylindrical structures could be treated as a rectangular bar.
(2) The variation in geometric size of the cross-section has remarkable influence on wave characteristics. The cut-off frequencies have a negative relationship with the cross-section area of the fan-shaped cylindrical structures. (3) For the FGPP fan-shaped cylindrical structures, the magneto-electric effect could be adjusted via altering the geometric size. (4) The phase velocities for higher modes at high frequencies approximatively approach a value which is below the shear velocity, and they also increase with the increase of β. (5) For the big wavenumber case, the Poynting vector distributions are concentrated in the region with more material of smaller elasticity modulus, and the Poynting vector distributions are also concentrated near the boundaries, especially near the corner.