Abstract

Stability and bifurcation behaviors for a model of simply supported functionally graded materials rectangular plate subjected to the transversal and in-plane excitations are studied by means of combination of analytical and numerical methods. The resonant case considered here is 1 : 1 internal resonances and primary parametric resonance. Two types of degenerated equilibrium points are studied in detail, which are characterized by a double zero and two negative eigenvalues, and a double zero and a pair of pure imaginary eigenvalues. For each case, the stability regions of the initial equilibrium solution and the critical bifurcation curves are obtained in terms of the system parameters which may lead to Hopf bifurcation and 2D torus. With both analytical and numerical methods, bifurcation behaviors on damping parameters and detuning parameters are studied, respectively. A time integration scheme is used to find the numerical solutions for these bifurcation cases, and numerical results agree with the analytic predictions.

1. Introduction

Functionally graded materials (FGMs) are extremely excellent engineering composite materials. Because of their specially tailored thermomechanical properties, FGMs are widely applied to large space station, shuttle, aircraft, automotive, and many others [13]. With the increasing use of FGM plates in engineering fields, research on the nonlinear dynamics of the FGM plates plays a significant role in applications. Holmes [4] studied flow-induced vibrations, bifurcations, and the chaos of thin plates. Yang and Sethna [5] used the averaging method to study the local and global bifurcations for parametrically excited vibrations of nearly square plates. Yu et al. [6] investigated the nonlinear vibrations and bifurcations of a thin plate with the aid of computation of normal form. Ye et al. [7] investigated the nonlinear vibrations and chaotic dynamics of an ant symmetric cross-ply laminated composite rectangular thin plate under parametric excitation. Zhang et al. [8] studied the bifurcation and chaotic motions of a composite laminated piezoelecttric rectangular plate with 1 : 2 internal resonance. Chang et al. [9] investigated the bifurcations and chaos of a rectangular thin plate with 1 : 1 internal resonance. Zhang et al. [10] studied the nonlinear dynamics of a parametrically excited simply supported laminated composite plate and found that there exist the multiple steady bifurcation solutions under the certain conditions.

Using asymptotic perturbation method based on the Fourier expansion and the temporal rescaling [11, 12], Zhang et al. [13] investigated the nonlinear oscillations and chaotic dynamics of a simply supported FGM plate. Yu and Huseyin [14, 15] studied the stability and bifurcation behaviors of the systems which were characterized by double zero eigenvalues, a simple pair of pure imaginary eigenvalues, and two pairs of pure imaginary eigenvalues. Yu [16] used a perturbation technique and computer algebra to compute the normal form of the nonlinear dynamical systems with two pairs of pure imaginary eigenvalues and analyzed double Hopf bifurcation. Hao et al. [17] analyzed the nonlinear dynamics of a simply supported FGM plate subjected to the transverse and in-plane excitations in a thermal environment. Zhang et al. [18] studied the stability and bifurcation behaviors of FGM which were characterized by a simple zero and a pair of purely imaginary eigenvalues and two pairs of pure imaginary eigenvalues. By the analytical approach, the present work investigates theoretically the stability conditions of the equilibrium solutions and their stability boundaries.

The objective of this paper is to investigate the local dynamic behaviors of a simply supported rectangular plate subjected to the transversal and in-plane excitations in a thermal environment. The resonant case considered here is 1 : 1 internal resonance and primary parametric resonance. Both analytical and numerical approaches are employed to consider the bifurcation and stability of this system. Two types of degenerated equilibrium points are studied in detail, which are characterized by a double zero and two negative eigenvalues and a double zero and a pair of pure imaginary eigenvalues, respectively. The stability regions of the initial equilibrium solution and the critical bifurcation curves are obtained in terms of the system parameters. Bifurcations leading to Hopf bifurcation and 2D torus and their stability conditions are also investigated. All numerical results agree with the analytic predictions.

This paper is organized as follows. In Section 2, the averaged equations of transverse motion of the FGM plate are given and the stability conditions of initial equilibrium solution are obtained explicitly. Section 3 is devoted to the studies on the dynamical behaviors of the system in the vicinity of the critical points: a double zero and two negative eigenvalues and a double zero and a pair of pure eigenvalues. Finally some conclusions are drawn in Section 4.

2. Formulation of the Problem

The paper focuses on the stability and bifurcation behaviors of a functionally graded material rectangular plate subjected to the in-plane and transversal excitations. The model is shown in Figure 1.

The paper focuses on the stability and bifurcations behaviors of a simply supported FGM rectangular under combined parametric and external excitations. The model is shown in Figure 1. The plate is of length , width , and thickness . Let and be the displacements of an arbitrary point and a point in the middle plane of the plate in the , , and directions; and represent the mid-plane rotations of a transverse normal about the and axes. The plate is subjected to an in-plane excitation uniformly distributed along the plate edges ,   , and a transverse excitation . Here and are the excitation frequencies of the transverse and in-plane excitations, respectively.

Applying the Galerkin procedure, the nonlinear dimensionless governing differential equations of transverse motion of the FGM rectangular plate are obtained as follows [17]:

It is assumed that the width-to-length ratio of the FGM rectangular plate is . Therefore only the primary parametric resonance and 1 : 1 internal resonance are considered in this paper. In this resonant case, there are the following relations: where is a small perturbation parameter, and are the first-order and second-order linear frequencies, and and are two detuning parameters. Using the asymptotic perturbation method, the averaged equation in the Cartesian form is obtained as follows [17]: where , , , and , where all coefficients can be found in [17]. The nonlinear functions are presented in the appendix.

The Jacobian matrix of (3a)–(3d) evaluated at the initial equilibrium solution is as follows: From which one may obtain the characteristic polynomial where

By the Routh-Hurwitz criterion, the equilibrium solution is stable, if the following conditions are satisfied: Conditions in (7) imply that all the eigenvalues of the Jacobian matrix (4) have negative real parts. When conditions in (7) are not satisfied, the initial equilibrium solution is unstable, and bifurcation may occur. In the next section, the detailed analysis will be given when conditions in (7) are violated.

3. Stability and Bifurcation Behaviors

3.1. Stability and Bifurcation Behaviors on the Damping Parameters

In this Section, the stability and bifurcation analysis on the parameters and are investigated, which can be divided into two cases due to degenerated equilibrium points.

3.1.1. The Case of Double Zero and Two Negative Eigenvalues

In this part, the characteristic polynomial (5) is supposed to have a double zero and two negative eigenvalues, ,   , which implies that , . One choice of the parameters that satisfy these conditions is , , , , , , ,   , , . Choosing and   as perturbation parameters and using the parameter transformation , and the state variable transformation one may transform (3a)–(3d) into a new system as follows: where the nonlinear functions are exhibited in the appendix.

The Jacobian matrix of (9a)–(9d) evaluated at the initial equilibrium solution at critical point is the following canonical form: The local dynamic behaviors of system (9a)–(9d) are characterized by the critical variables and . Furthermore, the bifurcation solutions for the noncritical variables and may be determined from (9a)–(9d) up to leading orders terms as

In order to study the bifurcation and stable properties of system (9a)–(9d) in the vicinity of the critical point, one only needs to analyze the following two-dimensional system: To find the stability conditions of the initial equilibrium solution , one may evaluate the Jacobian matrix of (11a) and (11b) at and obtain where

The characteristic polynomial is

The stability conditions for the initial equilibrium solution are

It is easy to see that unless , so if , the initial equilibrium solution is stable. Then a critical bifurcation curve is obtained:

When , Hopf bifurcation may occur. In fact, the eigenvalues of the Jacobian matrix at the initial equilibrium solution are

Let

It is easy to see that , , and , when and or . By Hopf bifurcation theorem, along , Hopf bifurcation may occur. The bifurcation diagram is shown in Figure 2. The initial equilibrium solution is stable while belongs to area I. Along , Hopf bifurcation may occur. When crosses and goes into area II, the initial equilibrium solution becomes unstable.

Numercial results have been obtained by using a time integration scheme-fourth-order Runge-Kutta method. The numerical computation is performed on the base of the original differential equations (3a)–(3d). For the parameter , which is located in the area II (the unstable region of the initial equilibrium solution), numerical results show that any solution starting from an arbitrary initial point diverges to infinity, initiating that the initial equilibrium solution is unstable, as predicted by the analytic study. When the parameter value is chosen from the stable region of the E.S., say, , a numerical solution starting from an initial point is obtained, which converges to the origin, implying that the E.S. is stable. This is shown in Figure 3, where the phase trajectories are projected onto the and plane. It should be noted that since the study is focused on the local dynamic behaviors of the FGM in the vicinity of a critical point, so the parameter should be chosen near the critical point .

If one chooses , which is located in the curve , Hopf bifurcation may occur. The phase trajectories starting from are shown in Figure 4.

3.1.2. Double Zero and a Pair of Pure Imaginary Eigenvalues

Choosing the following parameter values: , , , , , , , , , , which imply that , , then the Jacobian matrix (4) has the eigenvalues .

Choosing and   as perturbation parameters, and using the parameter transformation , , the characteristic polynomial (5) of the Jacobian matrix (4) becomes where,

The stability conditions for the initial equilibrium solution are

From the four inequalities above, one may get the following three transition curves: Then, the transition curves are shown in Figure 5.

When ,   and are all satisfied, the initial equilibrium solution is stable. When belongs to region I, the initial equilibrium solution is stable, and when belongs to region II, the initial equilibrium solution of the system is unstable.

Similar to the case in Section 3.1, different parameters are chosen to conform the analytical results obtained in this Section. When the parameter is chosen as , numerical results show that a trajectory starting from converges to the origin, implying that the E.S. is stable. The phase trajectories are projected onto the and subspaces as shown in Figure 6.

3.2. Stability and Bifurcation Behaviors on the Detuning Parameters

In this section, the stability and bifurcation analysis on the parameters and are investigated, which can be divided into two cases due to degenerated equilibrium points.

3.2.1. The Case of Double Zero and Two Negative Eigenvalues

Corresponding to Section 3.1.1, take parameter values , , , , , , , , , .

Choosing and as perturbation parameters, and then using the parameter transformation   , and the state variable transformation (7), one can get where the nonlinear functions are exhibited in the appendix.

Similar to Section 3.1.1, we can get the equations in the critical variables and , up to leading orders, as follows:

To find the stability conditions of the initial equilibrium solution , one may evaluate the Jacobian matrix of (24a)–(24d) at and obtain where

The characteristic polynomial is

The stability conditions for the initial equilibrium solution are

It is easy to see that unless , so if , the initial equilibrium solution is stable. From this equality, one can denote

For or , Hopf bifurcation may occur. The discussion for Hopf bifurcation is similar to that of Section 3.1.1. The bifurcation diagram for this case is shown in Figure 7. It is observed from this figure that the critical curves , separate plane into two parts. While is in area I, the initial equilibrium solution is stable. In area II, the initial equilibrium solution is unstable. Hopf bifurcation may occur along and .

By choosing , which represents a point in the parameter space located in the stable region I of the E.S., one can use the time integration to find a numerical solution. The trajectory starting from the initial point   converges to the origin shown in Figure 8, implying that the E.S. is stable.

If one chooses , which is located in the line , Hopf bifurcation may occur. The phase trajectories starting from are shown in Figure 9. Similar to the case of , the trajectories of are omitted here.

3.2.2. Double Zero and a Pair of Pure Imaginary Eigenvalues

Corresponding to Section 3.1.2, take parameters , , , , , , , , , . It yields that ,   and the Jacobian matrix (4) has the eigenvalues , .

Choosing and as perturbation parameters, and then using the parameter transformation , , then characteristic polynomial (5) of the Jacobian matrix (4) becomes where

Denote . It is easy to see that , unless . One may get the following two facts:(i)there exist no zero eigenvalues when since ;(ii)when , except for or , there exist two pairs of pure imaginary eigenvalues; when or , there exist a double zero and a pair of pure imaginary eigenvalues.

In fact, supposing that there exist two pairs of pure imaginary eigenvalues, we have , which follows or . For , the characteristic (31) becomes

The solutions of (33) are   ,   . It is easy to show that if or , there exist a double zero and a pair of pure imaginary eigenvalues. For the other case, there always exist eigenvalues with positive real parts, so the initial equilibrium solution is unstable.

By the analytic study, we know that when , except for or , there exist two pairs of pure imaginary eigenvalues (including resonance and nonresonance cases). Choosing , one can use time integration to find a numerical solution shown in Figure 10. A 2D torus is obtained when the trajectory starts from the initial point , which implies that there is a 2D torus for two pairs of pure imaginary eigenvalues.

4. Conclusions

In this work, the stability of a simply supported functionally graded material rectangular plate subjected to the transversal and in-plane excitations has been studied in detail to show rich dynamical behaviors in the vicinity of critical points. Hao et al. [17] studied the bifurcation and chaos of FGM rectangular plate with numerical methods. We mainly investigate the local stability and bifurcation of the FGM plate in the vicinity of the equilibrium point using analytical methods. When the stable conditions for the initial equilibrium solution are not satisfied, bifurcations including Hopf bifurcation and 2D torus may occur. The stable conditions, stable regions, and critical bifurcation curves for the steady state solutions are presented explicitly in terms of the system parameters. With both analytical and numerical methods, bifurcation behaviors on damping parameters and detuning parameters are studied, respectively. Numerical computations have been performed and shown for each of the bifurcation cases. All numerical solutions agree with the analytical predictions.

Appendix

The nonlinear functions in (3a)–(3d) are

The nonlinear functions in (9a)–(9d) are

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (11172125) and the National Research Foundation for the Doctoral Program of Higher Education of China (20133218110025).