Abstract
In this paper, the authors study the bifurcation problems of the composite laminated piezoelectric rectangular plate structure with three bifurcation parameters by singularity theory in the case of 1:2 internal resonance. The sign function is employed to the universal unfolding of bifurcation equations in this system. The proposed approach can ensure the nondegenerate conditions of the universal unfolding of bifurcation equations in this system to be satisfied. The study presents that the proposed system with three bifurcation parameters is a high codimensional bifurcation problem with codimension 4, and 6 forms of universal unfolding are given. Numerical results show that the whole parametric plane can be divided into several persistent regions by the transition set, and the bifurcation diagrams in different persistent regions are obtained.
1. Introduction
The bifurcation equations with multiple variables can be reduced to the bifurcation equation with single state variable by elimination method, but the work [1] showed that some bifurcation characteristics of system would be vanished. In addition, the bifurcation property of system with multiple bifurcation parameters is more complex than single parametric bifurcation feature [1]. Therefore, considering the system with all variables and parameters is close to the real physics problem. When an original steady-state system is subjected to small perturbations, all possible bifurcation behaviors can be revealed from the universal unfolding of bifurcation equations.
The singularity theory unifies the treatment of many diverse problems in steady-state bifurcation; such unification has the obvious advantage of elegance. The singularity theory was introduced to study the bifurcation problems by Golubitsky and Schaeffer [2] in 1985; the singularity theory has been greatly developed. Golubitsky and Keyfitz [3] studied the steady-state solutions for a continuous flow stirred tank chemical reactor by the singularity theory and obtained the qualitative bifurcation diagrams. Schaeffer and Golubitsky [4] applied the singularity theory to investigate the bifurcation phenomena of a model chemical reaction. Futer [5] and Sitta [6] considered symmetry of bifurcation parameters and studied the classification of bifurcation problems of codimension not greater than 1. Lavassani et al. [7] studied equivariant multiple parameter problems by singularity theorem and gave finite definite theorem and normal form of more parameters problem. Chen and Colleagues [8, 9] firstly proposed the C–L method to study the bifurcation of the periodic solution of nonlinear dynamical equation undergoing parameter excitation by combining the L–S method and the singularity theory. Seyranian and Mailybaev [10] gave multiple parameter stability theory and discussed the influence of system parameters on stability when parameters change. Jin and Zou [11] applied singularity theory to a restrained pipe conveying fluid and obtained the dynamical behavior in different persistent regions. Chen et al. [12] gave bifurcation analysis of an arch structure of parametric and forced excitation with codimension 5. Li et al. [13] extended Lyapunov-Schmidt reduction to fractional ordinary differential systems (FODSs) with Caputo derivatives. Li et al. [14, 15] gave bifurcation analysis of the Kuramoto-Sivashinsky equation in one-spatial dimension with three kinds of boundary value conditions and studied symmetry breaking bifurcation in -symmetric nonlinear large problems. In the applications of the two state variable singularity theory, the galloping of transmission conductor was studied. Nigol [16, 17] gave torsional galloping mechanism under galloping excitation. Yu [18] gave torsional feedback mechanism galloping excitation.
Piezoelectric materials are new type functional materials in engineering applications and can be utilized as the actuators and sensors in engineering structures; for example, it was used to construct morphing structures or morphing wings which can change the aircraft shape [19–21]. Many researchers have given considerable attention to the nonlinear dynamics of composite laminated plates. Pai and Nayfeh [22] presented a general nonlinear theory for research on the nonlinear dynamics of elastic composite plates undergoing moderate-rotation oscillations by considering the geometric nonlinearities. Shukla et al. [23] gave an analytical approach to analyze the nonlinear dynamic responses of a laminated composite plate with spatially oriented short fibers in each layer of the composite. Ye et al. [24, 25] dealt with the nonlinear dynamic characteristics of a parametrically excited, simply supported, symmetric cross-ply composite laminated rectangular thin plate, and a simply supported antisymmetric cross-ply composite laminated rectangular thin plate under parametric excitation. Lee and Reddy [26] investigated the nonlinear dynamic responses of laminated composite plates under thermostatically loading. Zhang et al. [27] investigated the nonlinear oscillations and chaotic dynamics of a parametrically excited simply supported symmetric cross-ply laminated composite rectangular thin plate with the geometric nonlinearity and nonlinear damping. Hao et al. [28] analyzed the nonlinear oscillations, bifurcations, and chaos of a functionally graded materials plate. Zhang et al. [29] established the governing equations of motion for the nonlinear oscillations of a simply supported symmetric cross-ply composite laminated piezoelectric plate subjected to the transverse, in-plane excitations, and the excitation loaded by piezoelectric layers and studied the periodic and chaotic dynamics of the composite laminated piezoelectric plate. Zhang et al. [30] studied the motion equations of simply supported flexible beam with five nonlinear terms under parametric excitation and analyzed the different forms of bifurcated response curve when the parameters are located in different regions. Zhang et al. [31, 32] studied the global nonlinear dynamics of thin plates under parametric excitation. The research on mechanic-electric couple responses of piezoelectric composite laminates plates has important theoretical significance and engineering application value. In the study of one-way mechanic-electric couple responses, Crawley et al. [33] presented the development and experimental verification of the induced strain actuation of plate components of an intelligent structure. Equations related to the actuation strains, created by induced strain actuators, to the strains induced in the actuator/substrate system were derived for isotropic and anisotropic plates. Plate strain energy relations are also developed. Wang et al. [34] analyzed laminated plates with distributed piezoelectric actuators by classical laminated plate theory. In studying of two-way mechanic-electric couple responses, Allik et al. [35] derived the finite element equations of piezoelectric body vibration. Lim et al. [36] studied the dynamic characteristics of isotropic plates with piezoelectric sensors and actuators by a 20-node block element, a 13-node transition element, and a 9-node plate element. Wang et al. [37] studied the optimization of thickness and embedding position of a piezoelectric actuator in composite laminates and analyzed an isotropic plate with isotropic actuator and a composite plate with anisotropic actuator.
The singularity theory of asymmetry dynamical systems with multiple state variables and single bifurcation parameter had been well developed by Golubitsky and Schaeffer. The universal unfolding of bifurcation equations was obtained [2], and it showed that the number of codimensions was equal to the number of auxiliary parameters. However, in asymmetry dynamical systems with multiple bifurcation parameters, the number of auxiliary parameters is less than codimension, which is contradiction to the previous conclusion. In this paper, the sign function is introduced to solve this problem, which leads to the fact that nondegenerate conditions of universal unfolding of the bifurcation equations are satisfied. It may determine the summation of the number of sign functions and the number of auxiliary parameters is equal to codimension. For the dynamical equations of the composite laminated piezoelectric rectangular plate structure in the case of principal parametric resonance and 1/2 subharmonic resonance for the first-order mode and primary resonance for the second-order mode [29], the four-dimensional nonlinear averaged equation is obtained by the Galerkin approach. By employed sign function, the universal unfoldings of the 1:2 internal resonance bifurcation equation are derived. The transition sets in the parameters plane are calculated and the bifurcation diagrams are depicted.
2. The Motion Equations and Perturbation Analysis
In this section, we summarize on modeling of composite laminated piezoelectric rectangular plate simply supported at four-edge in the literature [29]. The mechanical modal of the composite laminated piezoelectric rectangular plate is shown in Figure 1, where the edge lengths are and , respectively, and thickness is . The composite laminated piezoelectric rectangular plate is considered as regular symmetric cross-ply laminates with layers. Some of these layers are made of the PVDF piezoelectric materials as actuators, while others are made of fiber-reinforced composite materials. It is assumed that different layers of the symmetric cross-ply composite laminated piezoelectric rectangular plate are perfectly bonded with each other and with piezoelectric actuator layers embedded in the plate. A Cartesian coordinate system is located in the middle surface of the composite laminated piezoelectric rectangular plate. In this study, that and are assumed to represent the displacements of an arbitrary point and a point in the middle surface of the composite laminated piezoelectric rectangular plate in the , , and directions, respectively. It is also assumed that the in-plane excitations of the composite laminated piezoelectric rectangular plate are loaded along the direction at and the direction at with the form of and , respectively. The transverse excitation, which loads to the composite laminated piezoelectric rectangular plate, is represented by . The dynamic electrical loading is expressed as .

According to literature [29], applying the Reddy’s third-order shear deformation theory (TSDT), the Hamilton’s principle, the Galerkin method, and the boundary condition, a two degree-of-freedom nonlinear ordinary differential equation of the composite laminated piezoelectric rectangular plate is obtained as follows:where all coefficients in (1a) and (1b) are presented in Appendix.
To obtain a system suitable for the application of the method of multiple scales, the scales transformations may be introduced aswhere is a small parameter.
Substituting (2) into (1a) and (1b), we obtain the following dimensionless two degree-of-freedom nonlinear vibratory systems:
The method of multiple scales is used to find the uniform solutions of (3a) and (3b) in the following form:where and
The time derivative used in the method of multiple scales is given aswhere and .
We consider the case of principal parametric resonance and 1:2 internal resonance. In this resonant case, there are the following resonant relations:where and are two different linear natural frequencies, and are two detuning parameters.
Substituting (4a) and (4b)–(6) into (3a) and (3b) and balancing the coefficients of corresponding powers of on the left-hand and right-hand sides of equations, the differential equations are obtained as follows:
Order :
Order :
The solution of (7a) and (7b) in the complex form can be written aswhere and are the complex conjugates of and , respectively.
Substituting (9a) and (9b) into (8a) and (8b) yieldswhere represents the parts of the complex conjugates of the function on the right-hand side of (10a) and (10b) and represents the terms that do not produce secular terms.
Eliminating the secular terms from (10a) and (10b) yields
The functions and may be expressed in the polar formwhere and are the phase angles. Substituting (12) into (11a) and (11b) and separating the real and imaginary part, the four-dimensional averaged equation in polar form is obtained as
Making the left-hand of (13a), (13b), (13c), and (13d) equal to zero, eliminating in (13a) and (13b), and eliminating in (13c) and (13d) by using the relations between trigonometric functions, we obtained the frequency-response functions of the system as follows:
Expanding (14a) and (14b) leads to the following bifurcation equations:where
3. Singularity Analysis
We know that the universal unfolding of the bifurcation equations can reveal all possible bifurcation behaviors when the original system is subject to small perturbations. In this section, we will use the introduced sign function in Proposition 4 to establish the universal unfolding of the bifurcation equations (15a) and (15b).
Let where
In the following discussion, is adopted in .
3.1. The Restricted Tangent Space
Proposition 1. The restricted tangent space of the germ is given
Proof. According to Proposition 1.4 (see pages 169 in Golubitsky and Schaeffer II (1985)), is generated (as module over ) by the fourteen mappings: is generated by the twenty mappingsBy the following proof that there exists a reversible matrix between (20) and (21) when , we observe thatwherewhere In (15a) and (15b), if , bifurcation analysis is meaningless, so will be studied in this paper. When , substituting into (24a), (24b), (24c), (24d), (24e), and (24f), it can be obtained thatThis matrix equalswhereIn the matrix , it is easy to observe that the following rows and columns are linearly dependent: the 1th and the 5th rows; the 2th and 5th columns; the 10th, 11th, 12th, 15th, 16th, and 17th columns are identically zero. Therefore, the 5th row and the 5th, 10th, 11th, 12th, 15th, 16th, and 17th columns are eliminated. The original matrix has rank 13. The proof is completed.
3.2. The Recognition Problem
Proposition 2. Let where
If is a higher order term, then is strongly equivalent to .
Proof. By Proposition 1, a higher order term of is For polynomial , it is observed that Thus By applying Theorem 1.3 (see pages 168 in Golubitsky and Schaeffer II (1985)), we can obtain that is strongly equivalent to . The proof is completed.
According to Proposition 2, polynomial will be replaced by polynomial in the following discussion.
Proposition 3. The nondegenerate conditions of satisfy is equivalent to
Proof. According to , , and (see pages 402 in Golubitsky and Schaeffer I (1985)), we haveandSince , , we obtainWe rescale and ; the normal form of the germ can be expressed asThe proof is completed. It can be observed from (38) that its universal unfolding needs to complement the linear term about the state variables and constant term.
We improve Theorem 3.1 (see pages 409 in Golubitsky and Schaeffer I (1985)) by introducing a sign function . Proposition 4 is given as follows.
Proposition 4. Letbe a four-parameter unfolding of bifurcation problem in two state variables which satisfies the nondegeneracy conditions Proposition 3. Then is a universal unfolding of if and only ifandwhere , , and are bifurcation parameters and is auxiliary parameter when evaluated at .
Proof. Matrix can be expressed by six vectorsin whichSubstituting (34) into (43a), (43b), (43c), (43d), (43e), and (43f),In (44), only two vectors and are linearly independent. Since is a four-parameter universal unfolding of , there exist four complementary linearly independent vectors for . It is easy to find a complementary subspace to , spanned by where the values of , , and are one of , , and , which are the values different from each other, where are sign function, .
The introduced can satisfy (40) and (41). It is obtained from (45) Obviously, when , we haveSubstituting (46) into (41), when and , we can getConversely, only two vectors and in (44) are linearly independent. Because , the linearly independent vectors that complement need four. It can be adopted thatThus will be obtained. The proof is completed.
3.3. The Polynomial Space
Proposition 5. The polynomial space is given by
Proof. To derive the above equality, it should be noted thatwherewhereIt is found that can be expressed by three vectors , , and . Hence it can be obtained thatThe proof is completed.
3.4. The Universal Unfolding
Proposition 6. The codimension of in is 4, and the universal unfolding of the normal form (34) is given aswhere is auxiliary parameter and are sign functions, ; , .
Proof. According to Theorem 2.1 and equation 2.7 (see pages 211 in Golubitsky and Schaeffer II (1985)), it can be obtained thatIt is not difficult to see that the largest intrinsic ideal contained in is justandwhose dimension is .
Clearly, the dimension of is ; hence there exists a basis for a subspace of which is complementary to .
Consequently, it is easy to find a complementary subspace to , spanned by According to Proposition 4, simplified (59),It can be observed that is generated by the twelve mappingsEquation (60) makes two auxiliary parameters and for and , then and will occur. According to Proposition 4 and (61), and one of three vectors , , and are linearly independent. and one of three vectors , , and are linearly independent. Therefore, and will be eliminated.
Three sign functions , , and and an auxiliary parameter in universal unfolding of in (60), such as , , , and , will occur in universal unfolding of . Since , , , and are all parameters, according to Proposition 4, substituting four vectors , , , and into (41), only two vectors are linearly independent. By (41), it is easy to find will necessarily occur in matrix .
Thus, or or will be chosen as one of complementary bases of in .
In order to simplify the notation, is denoted as and is denoted as .
According to Proposition 4 and (61), there exist two linearly independent vectors as complementary basis of in . They are or when and are chosen as one of complementary basis of in . They are or when and are chosen as one of complementary bases of in . They are or when and are chosen as one of complementary basis of in .
There exist two sign functions and , such that or will occur when and are chosen as one of complementary basis of in . or will occur when and are chosen as one of complementary basis of in . or will occur when and are chosen as one of complementary basis of in .
According to the analysis above, the universal unfolding of is given byWhen sign functions and auxiliary parameter equal zero, namely, (; ) and , (62a), (62b), (62c), (62d), (62e), and (62f) can be obtained such thatUsing Proposition 4, when , it is obtained from (62a), (62b), (62c), (62d), (62e), and (62f) thatThe proof is completed. It can be found from (55a), (55b), (55c), (55d), (55e), and (55f) that codimension of in is equal to the summation of the number sign functions and the number of auxiliary parameter.
3.5. The Transition Set
The bifurcation behaviors of (55a), (55b), (55c), (55d), (55e), and (55f) are discussed in the following analysis. The transition set of the universal unfolding is formulated as follows (see pages 409 in Golubitsky and Schaeffer I (1985)).
The bifurcation variety set is
The hysteresis variety set is
The double limit point variety set is
Then the transition set is .
The equations of bifurcation points for (55a) are obtained as
By the third equation of (66c), it can be obtained thatand
Substituting (68) and (69) into (67), bifurcation occurs when
By the first and second equation of (66c), bifurcation occurs when
The equations of hysteresis points for (55a) are obtained as follows:where .
Now we assume and . Then, it can be taken thatand it is noted that if and only if
Equation (74) implies that
Multiplying (75) by and substituting it into (72c) we obtain
Using (72c) and (76), we obtain
Substituting (76) and (77) into (67), hysteresis occurs whenand
The equations of the double limit points for (55a) are obtained as
Using (70), double limit points occur when
, , and have equations as follows:
In the same way, transition set of the universal unfolding (55b), (55c), (55d), (55e), and (55f) may be obtained.
For (55c), a short calculation shows that bifurcation occurs when
By the first equation of (83c) and (84), a short calculation shows that bifurcation occurs when
By the second and third equation of (83c), bifurcation occurs when
The equations of hysteresis points for (55c) arewhere .
Now assume and , then it may be taken thatand note that if and only if
Equation (89) implies that
Next multiply (90) by and substitute (87c) to obtain
Finally, substituting (91) and (92) into (84), hysteresis occurs whenand
The equations of double limit points for (55c) are
By (85), a more involved calculation implies that double limit points occur for
We define , , and as follows:
Thus are a singular modification of cylindrical coordinates in space.
In these coordinates, substituting (97) into (85), (93), and (94), , , and have equations as follows:
From the above analysis, we deduce that (82a), (82b), (82c) and (98a), (98b), (98c) represent the relationship between two detuning parameters and a transverse excitation when the composite laminated piezoelectric rectangular plate structure subjected to small perturbations shows the bifurcation, hysteresis, and double limit points. If the disturbed system is equivalent to the original system, the space made up by these perturbations is the restricted tangent space. We know that the items and of complemented to in (55a) belong to the restricted tangent space and the variation of parameters and make no difference in the bifurcation property of the original system, while the variation of parameter of complemented to in (55a) has a qualitative influence on the bifurcation property of the original system. The similar analysis can be made in (55b), (55c), (55d), (55e), and (55f).
4. Numerical Results
The transition sets of the composite laminated piezoelectric rectangular plate structure with respect to small perturbations of the bifurcation parameter are depicted by (82a)–(82c) and (98a)–(98c) in Figures 2, 3, and 6, where implies bifurcation set, implies hysteresis set, and implies double limit points. It is observed from Figures 2, 3, and 6 that the amplitudes has a solution in region (1), which means that frequency-response equation of the original system has only one solution and the composite laminated piezoelectric rectangular plate structure subjected to small perturbations is in stable state. In region (2), there are multiple solutions, corresponding to the multiple solutions band of frequency-response equations of the original system. The composite laminated piezoelectric rectangular plate structure is in an unstable state under small disturbance, and there is a jumping phenomenon. There is no solution to region (3). It means that there is no solution to the frequency-response equations of the original system. The composite laminated piezoelectric rectangular plate structure with small disturbance is in a stable state.


By the changing of three bifurcation parameters, to understand bifurcation properties of the composite laminated piezoelectric rectangular plate structure subjected to small perturbations, Figure 4 is the bifurcation diagrams of taken the plus sign in (55a) when . In Figure 4(a), the parameters are chosen as , , , in which there are three limit points. In Figure 4(b), the parameters are chosen as , , , in which there is one limit point. In Figure 4(c), the parameters are chosen as , , , in which there is one limit point. Figure 5 is the bifurcation diagrams of taken the minus sign in (55a) when . In Figure 5(a), the parameters are chosen as , , , in which there are two limit points. In Figure 5(b), the parameters are chosen as , , , in which there is one limit point. In Figure 5(c), the parameters are chosen as , , , in which there is one limit point. By comparing the numerical results of the disturbed system (55a) above, it is easy to find that the bifurcation properties of the composite laminated piezoelectric rectangular plate structure under small disturbances in region (2) is significantly affected by the positive and negative signs in (55a).

(a)

(b)

(c)

(a)

(b)

(c)

Figure 7 is the bifurcation diagrams of taken the plus sign in (55c) when . In Figure 7(a), the parameters are chosen as , , , in which there is one limit point. In Figure 7(b), the parameters are chosen as , , , in which there is one limit point. In Figure 7(c), the parameters are chosen as , , , in which there is no limit point.

(a)

(b)

(c)
It is observed that jumping phenomenon and hilltop bifurcation occur at the critical point from region (1) to region (2), and only hilltop bifurcation occurs at the critical point from region (3) to region (2).
Meanwhile, by choosing different bifurcation parameters, the changing of the width of the multiple solutions band is studied. For taken the plus sign in (55a) when , detuning parameter is chosen as the controlling parameter to study the changing of the width of the multiple solutions band. In Figure 8(a), the parameters are chosen as , , . In Figure 8(b), the parameters are chosen as , , . In Figure 8(c), the parameters are chosen as , , . In Figure 8(d), the parameters are chosen as , , . It is observed from Figure 8 that parameter has a great influence on the width of the multiple solutions band. The parameters and are almost negligible for the width of the multiple solutions band when is fixed. In Figure 9(a), the parameters are chosen as , , . In Figure 9(b), the parameters are chosen as , , . In Figure 9(c), the parameters are chosen as , , . In Figure 9(d), the parameters are chosen as , , . It is observed from Figures 9(a) and 9(b) that parameter has a great influence on the number of solutions when and are fixed. The parameter has a great influence on the width of the multiple solutions band when and are fixed from Figures 9(c) and 9(d).

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)
Taking the plus sign in (55c) when , the detuning parameter as the controlling parameter to study the changing of the width of the multiple solutions band. In Figure 10(a), the parameters are chosen as , , . In Figure 10(b), the parameters are chosen as , , . In Figure 10(c), the parameters are chosen as , , . In Figure 10(d), the parameters are chosen as , , . It is observed from Figure 10 that parameters and have a great influence on the number of solutions when is fixed.

(a)

(b)

(c)

(d)
5. Conclusions
In this paper, the sign function was introduced, which leads to that the nondegenerate conditions of universal unfolding of the bifurcation equations of the composite laminated piezoelectric rectangular plate structure are satisfied. Then, we proved that bifurcation of the composite laminated piezoelectric rectangular plate structure in the case of 1:2 internal resonance is a high codimensional problem with codimension 4, and the universal unfolding of the bifurcation equations of the composite laminated piezoelectric rectangular plate structure was obtained.
The transition set divided the parameter spaces into three regions. The results indicated that the composite laminated piezoelectric rectangular plate structure with respect to small perturbations of the bifurcation parameter was in a stable state in regions (1) and (3) of Figures 2, 3, and 6. The composite laminated piezoelectric rectangular plate structure with respect to small perturbations of the bifurcation parameter was in an unstable state in regions (2) of Figures 2, 3, and 6, with jumping phenomenon or hilltop bifurcation. The results also showed that jumping phenomenon and hilltop bifurcation occurred at the critical point from region (1) to region (2), and only hilltop bifurcation occurred at the critical point from region (3) to region (2).
For this composite laminated piezoelectric rectangular plate structure, it was demonstrated that there exist abundant dynamic bifurcation patterns in a small perturbation sense. Some new dynamic bifurcation phenomena were also presented. Clearly, these results provided some inspiration and guidance for the analysis and dynamic designs of this structure.
Appendix
The coefficients of in literature [29] are presented as follows:in particular, for arbitrary constants , the expressionsare also solutions of in literature [29].
The coefficients of (1a) and (1b) are presented as follows:
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the support of the National Natural Science Foundation of China through Grant no. 11832002.