#### Abstract

This work presents the nonlinear dynamical analysis of a multilayer piezoelectric macrofiber composite (MFC) laminated shell. The effects of transverse excitations and piezoelectric properties on the dynamic stability of the structure are studied. Firstly, the nonlinear dynamic models of the MFC laminated shell are established. Based on known selected geometrical and material properties of its constituents, the electric field of MFC is presented. The vibration mode-shape functions are obtained according to the boundary conditions, and then the Galerkin method is employed to transform partial differential equations into two nonlinear ordinary differential equations. Next, the effects of the transverse excitations on the nonlinear vibration of MFC laminated shells are analyzed in numerical simulation and moderating effects of piezoelectric coefficients on the stability of the system are also presented here. Bifurcation diagram, two-dimensional and three-dimensional phase portraits, waveforms phases, and Poincare diagrams are shown to find different kinds of periodic and chaotic motions of MFC shells. The results indicate that piezoelectric parameters have strong effects on the vibration control of the MFC laminated shell.

#### 1. Introduction

MFC is a piezoelectric fiber material, consisting of monotonic piezoelectric material, epoxy matrix, and electrodes with a specific arrangement, which can be considered as homogenized orthotropic materials with arbitrary piezoelectric fiber angles like composite structures.

There are generally three types of piezoelectric material which researches focus on mainly. The first type of the piezoelectric fiber composite material is produced by Smart Material Corp, treated as 1–3 composite. The second type is originally developed by MIT, named active fiber composite actuators. Macrofiber composite (MFC) is referred to as the third type presented by NASA Langley Research Center. Since several applications require conformable and packaged piezoelectric structures or actuators, MFC material become widely applied in both academic and industrial field, which would be used for an easier integration in smart, intelligent, or adaptive structures [1–3]. Great advantages of the MFC material, which is able to respond to changing environment and controlling structural deformation, have led to a new generation for aerospace structures especially for morphing aircrafts [4].

Theoretical analyses are the basis to obtain dynamical characteristics of MFC material subjected to different excitations. In some researches, piezoelectric characteristics of MFC are studied by nonlinear constitutive equations and numerical simulation, including higher-order terms and corresponding coefficients [5–7] or by using strain dependent or effective piezoelectric coefficients in constitutive equations [8–11].

Recently, different finite element models are developed to analyze the nonlinear dynamical behavior of MFC laminated materials [12–16]. Gohari et al. [17] present an analytical solution to obtain static deformation and optimal shape control of smart laminated cantilever piezoelectric composite hybrid plates and beams under several coupled loads using piezoelectric actuators. Park and Kim [18] show that the aerothermal large deflection of the composite panel can be suppressed using MFC actuators and the excessive actuation of the MFC can cause snap-through phenomena. Kim et al. [19] study the vibration suppression of an end-capped cylindrical shell structure with surface bonded macrofiber composite actuators and analyze the dynamic characteristics of the cylindrical shell structure.

Dano and Jullière [20] investigate the use of macrofiber composite (MFC) actuators to actively control thermally induced deformations in composite structures. Korayem and Homayooni [21] analyze the influence of the applied voltage on vibration frequency of a multilayer piezoelectric microplate as well as the size effects on the macroplate at different boundary conditions. Li et al. [22] present the active control of random vibration for piezoelectric fiber reinforced composites laminated plates and discuss the effect of piezoelectric fiber orientation in the PFRC layers. Suresh Kumar and Ray [23] analyze the geometrically nonlinear vibrations of doubly curved smart sandwich shells integrated with a patch of 1–3 piezoelectric composites active constrained layer damping (ACLD) treatment.

Experimental studies are effective in verifying theoretical results. The effective properties of piezoelectric composites on the interface material are reported by using theoretical and numerical methods in [24, 25]. Andrianov et al. [26] focus on the evaluation of the homogeneous properties of the active layers in both and MFCs. Biscani et al. [27] study the effective electromechanical properties of MFC transducers and compare the material properties of the system with experimental data. Prasath and Arockiarajan [28, 29] present total packing effects on the overall properties of MFC using experimental models. Based on the developed Kirchhoff plate theory, in some articles [30, 31], the actuation responses of and MFCs integrated smart structures are studied subjected to transverse excitations. Kashiwao et al. [32] propose an optimization method of a vibration energy harvesting system which is made of MFC piezoelectric elements.

#### 2. Mechanical Model

MFC mainly consist of piezoceramic fibers, epoxy matrix, and electrodes, which have two different types of structures, named and modes. Here, a MFC middling thick cross-ply laminated shell is considered, where the piezomacrofiber is polarized in the thickness direction and fully embedded in the epoxy matrix as shown in Figure 1. The shell is a cantilever doubly curved cylindrical shell with a rectangular base. The orthogonal curvilinear coordinate is shown in Figure 2, and *α* and *β* curves are along the lines of curvatures on the middle surface and is perpendicular to the middle surface of the shell. Geometric dimensions of the shell in the curvilinear coordinate are the lengths and , and the thickness , and the principal radii of the curvatures and . The displacements of an arbitrary point within the shell in the coordinate () are denoted with , , and , respectively. is taken as positive vector going outward from the center of the smallest radius of the curvature. The shell is subjected to a transverse excitation given as , as shown in Figure 2.

According to Reddy’s third-order theory, the displacement fields at an arbitrary point in the composite shell are given in the following form:where , , and are the original displacements at the mid-plane of the MFC shell in the directions; and represent the rotations of transverse normal at the mid-plane about the and axes.

Using von Karman’s geometric relationship, the strain and displacements of the MFC shell can be written as follows: Here, , are the Lame coefficients of the shell and can be expressed as , .

The stress-strain relationships for the MFC material are given bywhere is the electric field intensity and represents the piezoelectric constant.

All piezoelectric fibers are considered to be poled in (through-thickness) direction. Therefore, it can be assumed that in-plane electric fields vanish (i.e., ), and the following reduced constitutive equations hold:where , , and , , are the stiffness modulus of the MFC shell, denote Poisson’s ratio, and represent the shear modulus.

The stiffness elements of the symmetric cross-ply composite laminated shell are expressed in terms of the stiffness coefficients as follows:

The relationship between curvilinear coordinate system and rectangular coordinate system is , , . Substitute these transformations into (1a), (1b), and (1c) and then apply Hamilton’s principle; the nonlinear governing equations of motion in terms of generalized displacements for the MFC shell can be obtained as follows: where , , and in (6c) is the damping coefficient, and all the parameters of the specific expression can be found in the Appendix.

The boundary conditions of the cantilever shell are expressed as

#### 3. Perturbation Analysis

Since vibration amplitudes of the lower frequencies are much larger than that of the higher frequencies for the shell, the mainly dynamical damage or instability of structures is caused by resonances in lower frequencies. Here, the first two modes of the MFC laminated shell are considered. Thus, displacements , , , , and , which satisfy the boundary conditions for the shell, are represented aswhere

To obtain the dimensionless equations, the transformation of variables and parameters are introduced as

Then, taking all these derived expressions in (8a), (8b), (8c), (8d), and (8e)–(9a), (9b), (9c), (9d), (9e), and (9f) into (6a), (6b), (6c), (6d), and (6e) and applying the Galerkin procedure, a two-degree-of-freedom nonlinear ordinary differential equation of the MFC laminated shell with dimensionless is obtained as follows:Here, , , , and present the piezoelectric coefficients of the shell.

Considering the case of primary parametric and 1 : 1 internal resonance, the following relations can be established accordingly:where and are the detuning parameters.

To study the oscillations and bifurcations of the established nonlinear system, the method of multiscale is a powerful tool to determine the solutions for conservative and nonconservative systems. Hence, this method is adopted to investigate the nonlinear vibration responses of the MFC shell.

The uniformly approximate solutions of (11a) and (11b) are obtained aswhere and .

Then, the derivatives with respect to become

Substituting (13a), (13b), (14a), and (14b) into (11a) and (11b) and equating the coefficients of to zero yield the following differential equations:

The averaged equations in the Polar Coordinates form are obtained for the MFC laminated shell as follows:

#### 4. Numerical Simulation

In this section, a series of numerical experiments are conducted for the nonlinear dynamical behavior of the MFC laminated shell. The nonlinear governing equations (17a), (17b), (17c), and (17d) are used to perform the numerical simulation through the Runge-Kutta algorithm. Whereas the transverse excitation is one of the most important factors on the nonlinear vibration of the system, the excitation is chosen as a controlling parameter in the subsequent studies.

Choose the following parameters: , , , , , , , , , , , , , , ,and . The frequency-response curves of the first-order and second-order are obtained through changing the initial condition of the MFC shell. The behavior shown in Figure 3 features the characteristics of soften-type systems with the piezoelectric coefficient increased from 0.5 to 1.5. It is also shown that the stiffness would be decreased with the increase in .

**(a)**

**(b)**

The bifurcation diagrams of Poincare sections for the displacements of the middle surface of the shell are shown in Figure 4 with the transverse excitation increasing from 700 to 2000. The complex nonlinear vibration responses of the system are discussed corresponding to the above parameters. The system keeps the chaotic motion until increased to 800, at which point period motions of the system occur. After short windows of period- motions occur in 950–1000, the motions of the system dominantly enter into chaotic motions thereafter. Then the chaotic motion is observed again until the transverse excitation equals 1500. The motion status of the MFC laminated shell alternately appears from unstable to stable and then again to unstable, which depends on the amplitude of the transverse excitation .

**(a)**

**(b)**

To reveal the specific form of different sections in the bifurcation diagram, the response-frequency curves, phase portraits, power spectrums, and waveforms of the shell are depicted as shown in Figures 5(a)–5(f). Diagrams (a) and (c) are, respectively, the wave forms on the planes and , and Diagrams (b) and (d) are the two-dimensional phase portraits on the planes and . Diagrams (e) show the three-dimensional phase portraits in the space while Diagrams (f) represent the Poincare diagram. The typical chaotic motion of the system is given as in Figure 5 when the amplitude of the transverse excitation arrives at 760.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

From the frequency-response curves of the system analyzed, it is known that the structure stiffness would alter when the piezoelectric coefficients of the system changed, so it would be effective to adjust the motion status of the MFC shell from unstable to stable through modulating the piezoelectric coefficients of the system. Fix the system parameters, based on which Figure 5 is depicted, and only adjust the piezoelectric parameters to perform numerical simulations and the results are shown in Figure 6. The bifurcation diagram reveals that the system would enter into periodic motions when is increased from 0 to 1.6. The specific periodic motion form of the system is also shown in Figure 7 when is equal to 2.

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

In order to describe the influence of the piezoelectric parameters on the nonlinear vibrations of the system profoundly, another set of parameters is selected in the following: , , , , , , , , , , , , , , , , and . It is found from the bifurcation diagram in Figure 8 that the response of the system varies from periodic motions to chaotic motions when the transversal excitation increases from 3 to 18. The different periodic and chaotic motions of the system are also given as in Figures 9-10, respectively, when equals 5 and 12.

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

Based on the motion of the system in Figure 8, we only increase the piezoelectric parameter to describe the vibration responses of the system. The bifurcation diagram of the system also shows the evolution from chaos to period- and then back into chaos again as shown in Figure 11. It is also found that the system enters into periodic motion from chaotic motion when the piezoelectric parameter increases from 1.0 to 2.0. Figure 12 shows the specific periodic motion of the system as equals 1.5.

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

From the results of the numerical simulation, it is discovered that the piezoelectric parameters could adjust the vibration responses of the structure effectively. Since the piezoelectric parameters could vary the stiffness of the shell, it is generalized: once the inherent frequency of the system changes, the resonance would be restrained for piezoelectric structures.

#### 5. Conclusion

In this paper, nonlinear dynamical behaviors of a simply supported cantilever MFC shell are presented, which is subjected to transverse loads. Based on known geometrical and material properties of its constituents, their electric field dependence is presented. The vibration mode-shape functions are obtained according to the boundary conditions, and then Galerkin method is employed to transform the partial differential equations into two nonlinear ordinary differential equations. The illustrative case of 1 : 1 internal resonances is considered. The externally excited system is transformed into a set of averaged equations by using the method of multiple scales. Next, the effect of the transverse excitations and the piezoelectric coefficients on the MFC laminated shell is described in numerical simulation.

The results of the numerical simulation demonstrate the complex nonlinear vibration responses of the MFC shell that occurred under the transverse excitation, including the periodic and chaotic motions. The energy transformation existed between two resonances modes. The appropriate control technique of the forcing excitations contributes significantly to the responses of autonomous nonlinear systems. It is also revealed that the piezoelectric parameters of the structure could adjust the dynamic stability of the structure from unstable to stable which would be a good way to control vibration responses for MFC structures.

#### Appendix

#### 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 (NNSFC) through Grant nos. 11572006, 11202009, 11072008, and 10732020 and the Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (PHRIHLB).