Dynamic Analysis of Multilayered Piezoelectric Quasicrystal Three-Dimensional Sector Plates with Imperfect Interfaces

: Piezoelectric quasicrystals have a tt racted extensive a tt ention due to their unique physical and mechanical properties. This paper studies the dynamic response of multilayered two-dimensional decagonal piezoelectric quasicrystal sector plates with imperfect interfaces. Based on the quasicrystal linear elasticity, partial di ﬀ erential state equations along the thickness direction are derived by using the state-space method. Then, by virtue of the di ﬀ erential quadrature method and the Fourier series expansions, this boundary-value problem with mixed boundary conditions and imperfect interfaces is solved. In addition, via the joint coupling matrix, the ﬁ eld quantities in the interior of the structure are connected to those on the external surfaces with numerical instability. Finally, parameter studies on the e ﬀ ects of angular spans, imperfect interfaces, and mixed boundary conditions are numerically investigated where the dimensionless frequencies and modes are exhibited.


Introduction
Quasicrystal (QC) was discovered by Shechtman [1] in the electron diffraction pattern of rapidly frozen Al-Mn alloy in 1982.QCs have excellent physical and mechanical properties such as high hardness, high conductivity, low thermal conductivity, and photoconductivity due to their unique atomic configuration [2,3].These properties make QCs widely useful in engineering as composite strengthening phase, thermal barrier coating, and solar thin film [4][5][6].As such, QCs have attracted great attention of researchers in many fields, such as mathematics, material science, crystallography, and physics.QCs can be divided into one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) QCs according to the numbers of the quasiperiodic directions [7].In the 2D QCs studied in this paper, the arrangement of atoms is quasiperiodic in two directions and periodic in the other direction.To understand the deformation behavior of QCs and describe their deformation process, QC elasticity is needed.The QC elasticity theory originated from condensed matter physics rather than traditional mechanics.According to Landau's density wave theory, there are two displacement fields in QCs: the phonon field and the phason field.Based on Landau's elementary excitation theory and the symmetry-breaking principle of condensed matter, Bak [8] developed the QC elasticity theory, whilst Lubensky [9] developed the elasticity of QCs based on the elementary excitation theory.
As for the dynamics of QCs, different with the certainty of statics, due to its unclear physical phenomenon of the phason field, there are some arguments on the dynamic behavior of QCs.Lubensky et al. [10] claimed that the phonon field and the phason field played different roles in the dynamic behavior of QCs because the phason field was insensitive to spatial translations, and the phasons were diffusive with very large diffusion times rather than propagating.Based on this, Lubensky et al. [9,10] proposed the hydrodynamics model of QCs.On the other hand, Bak [8,11] pointed out that the phonon field and the phason field behaved similarly in dynamics, and both followed the acoustic modes.Based on this point, Ding et al. [12] first discussed the motion equation of QCs where the law of momentum conservation held for both phonons and phasons, i.e., both the phonons and phasons are represented by wave propagations, named the elastodynamic model of wave type [13], and the equation of motion for both phonons and phasons followed the wave type.However, there was another voice that argued that the phonons followed the wave type while the phasons followed the diffusion type.Based on this point of view, Fan et al. [14] and Rochal and Lorman [15,16] compromised, respectively, the hydrodynamic model and the elastodynamic model to generate the elasto-hydrodynamics model of QCs.In this model, the equation of motion for phonons followed the wave type, while those for the phasons followed the diffusion type by replacing the mass density with the kinematic coefficient of the phason field.After that, Agiasofitou and Lazar [17] unified the elastodynamic model of wave type and the elasto-hydrodynamic model to establish the elastodynamic model of the wave-telegraph type.In this model, they claimed that the phonons were represented by undamped waves, while the phasons were represented by waves damped in time and propagating with finite velocity.Thus, the equation of motion for phonons followed the wave type, while those for the phasons followed the telegraph type.To sum up, these models based on different hypotheses were used by many researchers to collect the dynamic behavior of QCs.However, due to the lack of proper experimental data, all of these existing models possessed their reasons.In this paper, we used the elastodynamic model with the wave type to study the dynamic behavior of QCs.
Layered composites can be used to design the strength and stiffness of the involved structures according to the load requirements.As such, layered composites have become the main structural form of composites in industrial applications, while QCs, due to their material properties, usually serve as the thin film, coating layers, or reinforcement phase on advanced structures.Their unique properties can be used to improve the hardness and stiffness of the composites [18,19].As such, layered QC composites have become an attractive topic of research.By extending the pseudo-Stroh formalism, Pan and Heyliger [20] and Yang et al. [21] found exact solutions of multilayered multiferroics plates and 2D decagonal QC plates with simply supported boundary conditions.Guo et al. [22][23][24] obtained the analytical solutions for free and forced vibrations and static bending deformation of 1D QC nanoplates.Huang et al. [25] used the state-space method in their dynamic analysis of multilayered 2D piezoelectric QC (PQC) cylindrical shells filled with compressible fluid.Li et al. [26] deduced QC elasticity theory and elastodynamics for wave propagation.However, the pseudo-Stroh formalism [27,28] and the state-space method [29,30] are limited to the QC plates with simply supported boundary conditions.Plates with clamped or mixed boundary conditions cannot be solved by these two approaches.
Typically, the differential quadrature method (DQM) is employed to address static response and free vibration problems in plates subjected to arbitrary boundary conditions.In past decades, DQM was proven to be very efficient in solving differential equations with controlled boundaries or initial conditions [31].Chen et al. [32] applied the DQM to the structural free vibration analysis of composite beams under different boundary conditions.Also via the DQM method, Lu et al. [33] carried out the static analysis of composite beams with complex boundaries and introduced a new propagator matrix method to ensure numerical stability.Zhou et al. [34] used DQM along with the state-space method to derive the natural frequencies of multilayered plates under multi-field couplings with complex boundary conditions.The dynamic and static problems of various crystal plates and crystal beams under different boundary conditions were also solved [35][36][37][38].
Recently, Feng et al. [39] analyzed cubic QCs with imperfect interfaces under mixed boundary conditions.While the traditional state-space method can only solve problems with simply supported boundaries, the DQM along with the state-space method used in this paper can solve the problems with arbitrary boundary conditions.
Piezoelectric material is capable of converting electrical energy into mechanical energy and vice versa.As such, they are widely used in the design of piezoelectric sensors in biomedicine and other fields [40].Piezoelectric quasicrystal (PQC) materials combine the advantages of piezoelectric and QC materials and therefore offer superior performance compared to conventional piezoelectric materials.Huang et al. [41] analyzed the interlaminar stress of free-edge piezoelectric composite laminates.Khalid et al. [42] introduced a new loading condition when considering thermal-electric coupling and analyzed the free-edge interlaminar stress of smart composite laminates.Li et al. [43] constructed the exact solutions of 1D hexagonal QC with a piezoelectric effect.
A common assumption in the aforementioned works is that the extended displacement and traction vectors are continuous on the interface of two adjacent layers.However, most interfaces are imperfect as various damages could occur there (e.g., debonding, sliding, or cracking across the interface).So far, various imperfect interface models were proposed on heat conduction [44], piezoelectricity [45,46], and elasticity [47][48][49][50][51].In this paper, the imperfect interface with mechanical flexibility and weak (or high) conductivity proposed by Pan et al. [52] is used and verified in the 3D analysis of multilayered 2D decagonal PQC sector plates under mixed boundary conditions.Because the interface is not perfect in elasticity and dielectric properties, it is not easy to handle.Chen et al. [53] deduced the exact solution of a simply supported multilayered orthotropic piezoelectric rectangular plate with mechanical compliance and weak dielectric conduction interface by using the state-space method under isothermal conditions.The vibration analysis of circular, annular, and fan-shaped layered structures in the generalized boundary value problem is one of the classic problems [54].However, the reported study mainly focuses on the static response of QC laminated circular plates or cylindrical shells, while the vibration problem of the open sector plates under multi-field coupling with different angular spans is still limited.
In this paper, the dynamic analysis of a layered 2D PQC sector plate with mixed boundary conditions and imperfect interfaces is carried out.In Section 2, the governing equations for 2D PQCs are established based on QC linear elasticity, and the state equations are formulated as partial differential equations via the state-space method.These state equations are modified by using the DQM and the Fourier series expansions to meet different boundary conditions.In Section 3, the joint coupling matrix is derived to calculate the field quantities with numerical stability.In Section 4, the imperfect connection forms of the elastic field and the dielectrically weakly (highly) conducting interface are derived.In Section 5, the effects of different angle spans, imperfect conductive interfaces, and boundary conditions on 2D PQC laminated sector plates are studied.Finally, Section 6 draws the conclusions.

Problem Description
Shown in Figure 1 is the structure of a layered QC sector plate with imperfect interfaces in terms of the cylindrical coordinate system (r, θ, z).The atomic arrangement of 2D PQC is quasiperiodic along the r-θ plane and periodic along the z-direction.The total thickness is h, the inner diameter is ra, the outer diameter is rb, and the angle span is α.The thickness of any layer p (p = 1, 2, 3, …, M) is hp = zp−zp − 1 with the generalized spring model being used to simulate the imperfect interface between the layers.In this imperfect interface model, the tractions are continuous across the interface, whilst the displacements are not.According to the "spring factor" interface parameters, it is further assumed that the jump of the displacement component is proportional to the corresponding interface traction.For the dielectrically weakly conducting interface, the electric displacement is continuous, whilst the electric potential is not, with the jump in the electric potential proportional to the electric displacement.For the dielectrically highly conducting interface, the electric potential is continuous, whilst the electric displacement is not, with its discontinuity being a function of the electric potential.We further assume that the boundary conditions of the two opposite sides (θ = 0 and θ = α) of the sector plate are simply supported, while r = ra and r = rb can be simply supported or clamped.Additionally, for the case of free vibration, it is assumed that both elastic traction and vertical electric displacement are zero on the top and bottom surfaces of the plate.

Basic Equations
Based on the QC linear elasticity, the strain-displacement relationship of 2D PQCs is , , , where εij (i, j = r, θ, z) and wkl (k, l = r, θ) are the strains in the phonon and phason fields, respectively; ui and wk are the displacements in the phonon and phason fields, respectively; Ei stands for the electric field; and φ represents the electric potential. , where σij and Hkl are the stresses in the phonon and phason fields, respectively; Di stands for the electric displacements; Cij, C44, and C66 = (C11-C12)/2 are the elastic constants in the phonon field; K1, K2, and K4 represent the phason elastic constants; R1 represents the phonon-phason coupling elastic constant; ξii represent the dielectric coefficients; e15, e31, and e33 represent the piezoelectric coefficients in the phonon field; and d112 represents the piezoelectric coefficient in the phason field.
Based on Newton's second law and Bak's theory [8,11], the equations of motion of 2D PQCs in the absence of body force and electric charge densities following the elastodynamics model [12,13] can be expressed as follows: where ρ is the mass density and t is the time.

State Equations of a Homogeneous QC Sector Plate Layer
Substituting Equation (1) into Equation ( 2) and combining it with Equation (3), the partial differential state equations are where state variable ξ1 = [ur, uθ, uz, wr, wθ, φ, σzz, σθz, σrz, Hθz, Hrz, Dz] T , in which the superscript "T" denotes the transpose of a vector or a matrix; and D stands for the state matrix.
For the electro-elastic problem of a 2D PQC sector plate, the simply supported and electrically grounded boundary conditions of the two opposite sides θ = 0 and θ = φ are 0, 0, 0.
Then, the general solutions of displacements, stresses, electric potential, and electric displacements satisfying Equation ( 5) are (for the time-harmonic case) as follows: , sin where ω is the angular frequency, and imaginary i = √-1; p = nπ/α, in which n is the halfwave number.
In the radial r-direction, we assume more general boundary conditions, including the clamped or mixed boundary conditions.Therefore, in the r-direction, the DQM is used to solve the partial differential equations.Assuming that the function f(r) is smooth and continuous in the interested r-interval, in terms of the DQM method, the approximate value of the m-order derivative of the function f(r) at ri is ( ) where X ij (m) are the differential quadrature weight coefficients and N is the number of discrete points in the r-direction.
The form of Chebyshev-Gauss-Lobatto discrete points in a cylindrical coordinate system is given as Then, the stresses, displacements, electric potential, and electric displacements in Equation ( 4) are expanded by the Fourier series and DQM, resulting in the following state equations: (1) 13 1 , (1) 10 12 13 14 14 2 (1) 14 14 14 14 2 2 1 1 (1) 14 1 1 , ( ) 1  (2) (1) where i = 1, 2, …, N, and coefficients am (m = 1, 2, …., 18) are listed in Equation (A1) in Appendix A. Equations ( 9) and (10) represent the state equations at the endpoints when i = 1 and N. Therefore, the simply supported and clamped boundary conditions can be expressed directly by displacements, stresses, and electric potential: where d = 1 and N. To simplify presentations, the following abbreviations are introduced to express the boundary conditions of the QC sector plate: SSSS, SCSS, SSSC, and SCSC, where for example, 'SCSS' indicates that the boundary condition is simply supported (S) at r = rb, clamped (C) at r = ra, and simply supported (SS) along the edges θ = 0 and θ = α.Therefore, in general, for the given four sets of boundary conditions: (SSSS, SCSS, SSSC, and SCSC), we can arrive at the state equation for any given layer.Taking a 2D PQC sector plate with SCSS boundary conditions as an example, the state equations are shown in Equations (A2) and (A3) in Appendix A.
As such, Equations ( 9) and (10) of Layer p can be simplified to matrix form: where ( ) , , , , , , , , , , , u u u w w φ σ σ σ H H D is the overall state vector and each sub-vector is composed of state variables at all discrete points (u r T =  ri ); and P (p) is the coefficient matrix of Layer p.

General Solutions for a 2D QC Sector Plate
Since the state equations represent a system of homogeneous linear ordinary differential equations with constant coefficients, their general solutions can be obtained directly.From the general solutions, the propagation relationship of the state vectors on the top and bottom surfaces of sector plates can be obtained by combining the interface conditions between the layers.
The general solutions of the state equations Equation ( 13) can be expressed as where zp-1 and zp, respectively, represent the coordinates of the lower and upper interfaces of Layer p.If the state vector of the lower surface of Layer p is known, the state vector at any thickness position in this layer can be obtained by Equation (14).In particular, for z = zp, we have where , subscripts 0 and 1 represent the lower and upper interfaces of Layer p.Similarly, for Layer p + 1, we have If we assume that these two layers are perfectly connected, then all variables are continuous at the interface Substituting Equation (17) into Equations ( 15) and ( 16) to eliminate the state variables on their common interface zp we obtain We first assume that all the interfaces are perfect.Then, the state vectors at all interfaces can be eliminated in the same way to obtain the relationship of the state vectors on both the top and bottom surfaces of the layered plate ( ) ( ) where and Equation ( 20) is called the global propagator matrix.For convenience, Equation ( 19) can be written in a block matrix form as follows: ( ) ( ) where Then, from Equation (21), we have ( ) ( ) ( ) Since the tractions and electric displacements on the top and bottom surfaces of the layered plate are zero for the free vibration problem, Equation ( 23) is simplified as ( ) The above equations are a system of homogeneous linear equations, and the coefficient matrix is a function of the dimensionless frequency Ω.A nontrivial solution requires that the determinant of the coefficient matrix in Equation ( 24) be zero, i.e., For the solved dimensionless frequency Ω from Equation ( 25), the corresponding eigenvectors (namely, the eigenmode of each state vector) can be obtained from Equation (19).
For the forced vibration problem, the top surface of the QC sector plate is subjected to a uniform normal mechanical excitation, namely where  represents the uniformly distributed excitation on the top surface.Substituting Equation (26) into Equations ( 22) and ( 23), we obtain ( ) ( ) The displacements and electric potential on the bottom surface under forced vibration can be obtained from Equation ( 28) With these solutions, the state vectors at any thickness position in the laminated sector plate can be obtained using Equation (19).

Analysis of Imperfect Conductive Interfaces
An imperfect interface with mechanical flexibility is used in QC laminated sector plates.It is assumed that the traction at the interface is continuous, while the displacements are discontinuous.The imperfect condition on the interface for the elastic field is  The imperfect condition of the electric field can be divided into two cases [43]: dielectrically weakly conducting and dielectrically highly conducting.In the case of weakly conducting, it is assumed that the electric displacements on the interface are continuous, whilst the electric potential is not.In the case of highly conducting, the electric potential is continuous, whilst the electric displacements are not.In terms of mathematical forms, they are as follows: Dielectrically weakly conducting: .
where I is the identity matrix.
In the case of dielectrically weakly conducting In the case of dielectrically highly conducting If the plate is perfectly connected, G 12 (p) = G 21 (p) = 0.
Now as an illustration, we further assume that where K is defined as the interface compliance, and it is used to characterize the degree of imperfection of the interface.With the formulations above for the imperfect interface case, Equation ( 19) for the perfect interface condition can be concisely written in the following simple form for the general interface cases as follows: To summarize, the state-space method is employed to obtain the partial differential state equations from the governing equations of 2D PQCs.These state equations are transformed to an ordinary differential form by using the DQM and the Fourier series expansions.Then, the joint coupling matrix and the imperfect interface connections are considered and combined with the propagator matrix to obtain the state vectors at any thickness position in the laminated sector plate.

Numerical Examples
The tractions and electric displacements on the top and bottom surfaces of the 2D PQC sector plate are assumed to be zero.The effects of different angular spans and interface compliances on the dimensionless frequencies and modes of the sector plate in free vibration, as well as the effects of different boundary conditions on the displacement and stress under mechanical loads are studied.The dimensionless form of field variables and frequencies is defined as where Cmax, emax, and ρmax are the maximum values of the phonon elastic coefficient, phonon piezoelectric coefficient, and density of the 2D PQC sector plate, respectively.In this section, the aspect ratio of the three-layer QC sector plate is fixed at H/(rb-ra) = H/R = 0.03, with each layer having the same thickness.The sector plate is composed of two different 2D QC materials.The parameters of these two materials are given in Table 1, where R1 = R2 = R3 = R5 = R6.The discrete points in all examples in this section are taken as N = 13 to meet the needs of computational accuracy and convergence.The radial coordinates of data points in the resulting graph of this section are fixed at r * = 1.25.The ring coordinates of data points are all basic solutions, and the angle variation can be multiplied by the sinecosine factor from Equation ( 6).

Effect of Different Angle Spans
We first examine the influence of different angular spans (α = 60°, 90°, 120°, 150°, 180°) on the free vibration characteristics of QC1/QC1/QC1 sector plates with SSSS boundary conditions.The interface compliance of the QC sector plate is set as K* = 0 (for the perfect interface case).Table 2 shows the first nine dimensionless frequencies Ω of the QC sector plate with different angle spans.It can be seen from the table that with increasing α, the dimensionless frequencies Ω decrease, which indicates that the variation of α changes the stiffness of laminated sector plates and leads to a change in structural stability.Figure 2 shows the variation of the phonon and phason displacements and the electric potential of the QC sector plate along the thickness direction at the first dimensionless frequency for the five different angular spans.It is observed that u * (Figure 2a) is insensitive with varying α and it is antisymmetric with respect to the midplane.The closer to the top and bottom surfaces, the greater the absolute value of u * .u θ * , u z * , w * , w * , and φ * (Figure 2b-f) is.Also, u * is antisymmetric with respect to the midplane, with its absolute value at the top and bottom of the plate decreasing with increasing α. u * does not change with thickness, but changes with α, indicating again that a decrease in the angular span results in an increase in the stiffness of the QC plate.The absolute value of w * along the z-direction decreases with increasing α, and the absolute value of w * also increases with increasing α.The variation of φ * becomes less evident as α changes, indicating a decrease in electric potential energy due to free vibration.

Effect of Different Interface Compliances
In this subsection, the influence of the interface compliance K* on the free vibration characteristics of QC sector plates with SCSC and zero electric potential boundary conditions is investigated.To simplify the model, we assume that (G1)33 = (G1)34 = (G1)43 = (G1)44 = 0, while other interface compliances are defined as (G1)11 = (G1)22 = (G2)55 = K* = 0, 0.3, 0.6, 0.9.In addition, the angular span of the QC sector plate is fixed at α = 30˚.Table 3 shows the first nine dimensionless frequencies of the free vibration of the QC2/QC1/QC2 sector plate with a dielectrically weakly conducting interface.It can be seen from Table 3 that K* has a significant effect on the eigenfrequency Ω.In the same mode, Ω decreases with increasing K*.Since increasing the interface compliance will reduce the stiffness and strength of the laminated QC sector plate, this could cause fatigue damage in the structure.Furthermore, a weakly conducting interface also influences the dimensionless frequency of the QC sector plate.
Figure 3 shows the variation of the phonon and phason displacements and the electric potential of the QC sector plate along the thickness direction at the first dimensionless frequency with different interface compliances when the interface is weakly conducting.It is noted from this figure that when K* = 0, which corresponds to the perfect interface case, the displacements and electric potential are continuous between layers.If K* ≠ 0, u r * , u * , w * , w * , and φ * become discontinuous between the layers.It is observed that the variation trends of u r * and u * (Figure 3a,b) for different K* are opposite to each other, but they are both antisymmetric with respect to z/R = 0. ) ) Table 4 shows the first nine dimensionless frequencies of the free vibration of QC2/QC1/QC2 sector plates when the interface is highly conducting.It can be seen from this table that interface compliance has a significant effect on each eigenfrequency Ω.For the frequency on the same mode, Ω decreases with increasing K*. Figure 4 shows the variation of the phonon and phason displacements and the electric potential of the QC sector plate along the thickness direction at the first dimensionless frequency for different interface compliances when the interface is highly conducting.The variation trends u r * and u * (Figure 4a,b) are opposite to each other for varying K*, and they are antisymmetric with respect to z/R = 0. Displacements and electric potential are continuous if K* = 0. u z * (Figure 4c) is nearly constant along the plate thickness for different K*, varying slightly around 0.27.It can be seen from Figure 4d,e that w r * and w * are antisymmetric with respect to the midplane, and the vibration mode shape changes with changing K*.The electric potential is continuous under the condition of the highly conducting interface, which is different from Figure 4f.

Effect of Different Boundary Conditions
In this subsection, we study the influence of four different boundary conditions on the modes of QC1/QC2/QC1 sector plates with perfect interfaces under the action of uniformly distributed mechanical excitation on the top surface.Other parameters of the QC sector plate are fixed as α = 60˚ and K* = 0.In Equation ( 26), the dimensionless external excitation frequency is fixed at σ 0 = 1, which is greater than its first dimensionless eigenfrequency as listed in Table 4.
Figure 5 shows the variation of the phonon and phason displacements and electric potential along the thickness direction under four different boundary conditions on the sides for θ = 0 and θ = α (with boundary conditions on the radial sides being fixed as SS): SS, CC, SC, and CS.According to Figure 5, the influence of simply supported and clamped boundary conditions on field variables is different.In Figure 5a, under the SC boundary, the closer to the top and bottom surfaces, the greater the amplitude of u * is, while there is nearly no change under the other three boundary conditions.In Figure 5b, the clamped boundary reduces the amplitude of u * , and the amplitude under the CC boundary is almost zero.In Figure 5c, u * is nearly constant along the thickness direction, and its amplitude under the CC boundary condition is nearly zero.In Figure 5d, e, the SC boundary condition has the greatest influence on the amplitude of w r * and w * , and the amplitude of phason displacements under other boundary conditions is almost zero.In Figure 5f, the electric potential is also mostly affected by the SC boundary condition with its amplitude being the largest in the midplane, while its amplitude under the other boundary conditions is very small (especially for the CC, it is zero).Figure 6 shows the variation of the phonon and phason stresses and electric displacements along the thickness direction under the four different lateral boundary conditions.The modes of stresses and electric displacements are zero on the top and bottom surfaces of the sector plate, which meets the assumption that the traction and electric displacements on the top and bottom surfaces are zero.It is observed from Figure 6 that among all boundary conditions, the influence of simply supported and clamped boundary conditions on field variables is completely different, with the SC boundary having the greatest influence on the change in field variables along the thickness direction.The amplitude of stresses and electric displacements under SS, CC, and CS boundaries change slightly (with the amplitude of the CC case being the closest to zero).The stresses and electric displacements caused by forced vibration under the SC boundary change most obviously and they further fluctuate as sinusoidal curves.In Figure 6a,d-f, the amplitude reaches the extreme value at the interface, and in Figure 6b,c, the amplitude reaches the maximum value in the midplane.

Summary
In this paper, the dynamic response of a layered 2D PQC sector plate with an arbitrary angular span is studied using a semi-analytical method.The state equations of a 2D PQC sector plate are obtained by using the state-space method based on the QC linear elasticity.The DQM and Fourier series expansions are used to satisfy the simply supported and clamped boundary conditions.The propagator matrix method is used to combine the ordinary differential form equations with the interface matrix to obtain the solutions of all field variables.The semi-analytical method proposed in this study is suitable for analyzing the vibration of 3D plates with arbitrary layers and sectorial shapes, which includes the special case of a layered rectangle.In the numerical examples, the effects of different angular spans, different interface compliances, and different lateral boundary conditions on the dimensionless frequencies and field variables distribution of QC sector plates are discussed.The main conclusions are listed as follows: 1.The variation of the sector angle α changes the stiffness of laminated QC sector plates, and consequently, it results in different modes for different α; 2.An increase in interface compliance will reduce the stiffness and strength of the laminated QC sector plate, and the displacements in the r-and θ-directions (in addition to u z * ) are discontinuous from one layer to the other due to the interface imperfection; 3. Out of all the boundary conditions, the impacts of simply supported and clamped boundary conditions on field variables are notably distinct.The mixed case where one side is simply supported and the other is clamped has the most significant effect on the amplitude variation of field variables along the thickness direction.(A1) The state equations under SCSS boundary conditions are                  , ,

Figure 1 .
Figure 1.The structure of a layered QC sector plate with imperfect interfaces.
" and "+" are the two sides of the common interface zp between Layer p + 1 and Layer p; and α a (p) and β b (p) (a = 1, 2, 3, b = 1, 2) are the interface-stiffness coefficients of the phonon and phason fields on the interface.It is noted that α a (p) , β b (p) ⟶ 0 indicates a completely separated interface, while α a (p) , β b (p) ⟶ ∞ denotes a perfectly connected interface.

Figure 3 .
Figure 3. Variation of the phonon and phason displacements and electric potential along the thickness direction for different interface compliances K* for the weakly conducting interface: (a) u * , (b) u * , (c) u z * , (d) w * , (e) w * , and (f) φ * .

Figure 4 .
Figure 4. Variation of the phonon and phason displacements and electric potential along the thickness direction for different interface compliances K* for the highly conducting interface: (a) u * , (b) u * , (c) u z * , (d) w * , (e) w * , and (f) φ * .

Table 2 .
First nine dimensionless frequencies Ω of the QC sector plate with different angle spans.

Table 3 .
First nine dimensionless frequencies Ω of the QC sector plate with a dielectrically weakly conducting interface.

Table 4 .
First nine dimensionless frequencies Ω of the QC sector plate with a dielectrically highly conducting interface.