Abstract

In this paper, on the basis of taking the von Karman nonlinear factors into consideration, the constitutive equation of the antisymmetric cross-ply laminated composite was used to calculate the internal force and internal moment of the bistable structure, and the dynamic equilibrium equation and the compatible equation were constructed, respectively. The two equations were combined to establish a nonlinear dynamic model for the antisymmetric cross-ply laminated glass fiber resin bistable shell. Then, the finite element numerical simulation software ABAQUS was adopted to perform simulation modeling and numerical analysis on a series of bistable specimens, so as to study the impact of different geometric parameters on the frequencies, mode shapes, and other vibration characteristics of the antisymmetric laminated fiber resin bistable shell. Galerkin discretization was conducted on the vibration partial differential equation. Since there are only even-order partial differential terms of deflection with respect to x and y in the vibration partial differential equation at this time, the form of series obtained by each term is the same, which simplifies the discretization of the dynamic equilibrium equation and the compatible equation. Finally, the two equations after discretization were merged to obtain the three-degree-of-freedom nonlinear ordinary differential equation of the antisymmetric cross-ply laminated glass fiber resin bistable shell. The system averaged equation was acquired by perturbation analysis through a multiscale method, and the periodic solution of the antisymmetric laminated bistable system was studied. Moreover, the system’s nonlinear dynamic behavior characteristics such as bifurcation and chaos were explored when the main resonance is close to and , respectively, and the internal resonance is 1 : 2 : 3.

1. Introduction

A bistable structure is formed by laying multiple layers of fiber-reinforced materials, which is a composite laminated structure with two different stable states. Generally, under the driving of external loads such as mechanical force, smart materials, and temperature field, it can be changed from one stable state to the other, and its stable configuration can be maintained without the need for sustained energy input. Hyer [1] first discovered that asymmetric laminated cylindrical shells have bistable characteristics. In 1996, Daton-Lovett [2] found that antisymmetric laminated composite cylindrical shells also present the two stable states of a cylinder, which attracted considerable attention from many scholars, and relevant theoretical, numerical, and experimental studies are successively conducted. Innovatively, this paper found the law of deformation cloud diagrams of the antisymmetric laminated shallow shell and cylindrical shell made of T300 and studied the influence of geometric parameters (e.g., length, initial cross-section radius, initial central angle, number plies, and ply angle) on the modal frequency and mode shape of the antisymmetric laminated cylindrical shell.

The theoretical study of antisymmetric bistable laminated structures is mainly focused on statics and dynamics. In statics, Iqbal [3] proposed a simple linear elastic bistable model based on the classical laminate theory and predicted the bistable behavior of antisymmetric laminated shells according to the principle of minimum potential energy. Galletly and Guest [4,5] successively proposed the beam model and the shell model regarding the theoretical model developed by Iqbal. The theoretical results are in good agreement with the finite element simulation results, but greatly different from the experimental results. Pellegrino and Guest [6] proposed the double-parameter model by referring to previous studies. This model can be used to calculate and analyze the bistable characteristics of antisymmetric and symmetric laminated isotropic cylindrical shells as well as those with prestress. The College of Engineering, University of Cambridge successively studied the bistable characteristics of multistable composite shells with surface folds [710] and isotropic laminated shells with prestress [11]. Concerning the studies on the nonlinear dynamics of bistable structures, Zhang and Zheng [12] adopted the extended high-dimensional Melnikov method and numerical method and found that the multipulse double-parameter chaotic motion of asymmetric bistable laminated square plates will occur under a certain coupling effect of external excitation and parametric excitation. Zhang and Liu [13] studied the dynamic jump phenomenon of asymmetric bistable laminated square plates and its nonlinear vibration under base excitation from both theoretical and experimental aspects. Fei [14] analyzed the jump behavior of a two-layer orthogonal bistable laminated plate, while considering factors such as thickness, temperature, and external excitation of the laminated plate.

For the numerical analysis, the stable state transition of bistable laminated structures can be caught by the commercial finite element software, so the calculation result is more intuitive. Also, numerical analysis can be performed to perfect the theoretical model and predict the bistable characteristics of the structure, thus providing references for the experimental study. Iqbal and Pellegrinot [15] conducted numerical simulation on the bistable behavior of antisymmetric laminated shells using ABAQUS software and analyzed the stable state transition of antisymmetric laminated cylindrical shells and the stress distribution of the second stable state structure. Huang [16] introduced the bistable differential system into each cell of the grid and studied the effect of topology on the stable-state trend of the network.

In terms of the experimental study, Dano and Hyer [17] measured the jump loads of the bistable laminated plate by adjusting the water load, and the experimental measurements are in good agreement with the theoretical calculations. Potter et al. [18] adopted the loading method similar to “three-point bending” to induce the jump phenomenon of the bistable laminated plate and measured the load-displacement curve of the center of the form of the laminated plate. Daynes and Tawfik [19, 20] adopted similar experimental loading schemes and obtained the jump loads of orthogonal laminated plates with different geometric dimensions and the relations of curvature variation between two stable states. Etches [21] considered the effect of humidity on the mechanical properties of a bistable laminate plate and found that the jump load of the asymmetric orthogonal bistable laminated plate changes significantly at the initial stage of the experiment and then gradually levels off. Dano and Hyer [22] realized the intelligent control of asymmetric bistable laminated plates by controlling the temperature of the alloy wires to generate different driving forces. Schultz et al. [23] utilized a bistable laminated plate driven by piezoelectric materials to realize the stable state transition, but there is a big error between experimental and theoretical results. Tawfik et al. [24] realized the mutual transition between the two stable states by applying voltage, and the error between experimental and finite element simulation results is less than 10%. Giddings et al. [25] used piezoelectric patches to control the deformation of bistable laminated plates and studied the relationship between voltage and deformation, which played a guiding role in the design and control of the intelligent deformable structure. Kim et al. [26] realized the mutual conversion of bistable cantilever beams in two stable states by using MFC and SMA.

In the research of practical engineering application, the University of Cambridge [15] studied the bistable thin-walled tube similar to a tape measure and realized the stable state transition by controlling the energy accumulated in the folding process, which enables the bistable tap structure to be used for automatic hinge locks, detectors, expandable structures, and other aerospace devices. Schlecht [27] applied the bistable torsional assembly structure to the wing configuration of unmanned aerial vehicles, which in turn, can improve the sensitivity of aircraft during rotation because it requires only a small amount of energy to drive the wing to deformation. Diaconu and Weaver [28] used a bistable structure to design a wing model that can meet the actual conditions and analyzed the application prospect of such wing as a variant structure. Mattioni [29] made variable wing structures using bistable structures and used intelligent components to control their deformation to change the wing structure. The bistable structure can serve as both a deformable structure and as an energy harvesting structure. By utilizing the nonlinear vibration in the stable-state jump process, the bistable structure can be shaped into a broadband energy harvester [3032]. Three bistable structures have been applied to the piezoelectric energy capture, which is formed by using the principle of homopolar repulsion of magnetic force [3338], buckling of the beam, and thermal stress [39, 40], respectively, and the conversion of mechanical energy to electrical energy is achieved during vibration. Yao [4145] separated the end of the upper piezoelectric layer from the base layer in the traditional piezoelectric cantilever beam and conducted an experimental study on broadband energy collection and dynamic response of the L-shaped piezoelectric cantilever beam. Chen [46, 47] studied a micro/nanoscale bistable plate electrical-thermal-mechanical coupling system for energy harvesting and analyzed the voltage-modal frequency response of a nonlinear system.

2. Dynamic Modeling

The antisymmetric cross-ply laminated structure is formed by alternately laying the orthotropic single-layer structure at an angle of and between the main direction of the material and the coordinate axis. The stiffness coefficient of the antisymmetric cross-ply laminated structure has the following relationship: , , and . Thus, , , and , and it can be proved that and . In this study, four layers antisymmetric cross-ply laminated bistable thin shell was selected as the model, as shown in Figure 1, and the specific laying mode is shown in Figure 2. Glass fiber and epoxy resin substrate were used as the materials for the shell. The thickness of a single layer is 0.185 mm. Other material parameters are displayed in Table 1.

Based on the classical nonlinear shallow shell theory and the von Karman large deformation theory, the nonlinear response of the antisymmetric cross-ply laminated bistable shallow shell was studied. The nonlinear geometric equation of the middle plane of the Donnell cylindrical shell was applied:where , , and represent the displacement of the point on the middle plane of the cylindrical shell in the directions of , , and and and represent the initial radius of curvature of the cross-section in the and directions of the cylindrical shell.

Since the main direction of vibration of the shallow shell is transverse deflection, the inertia in the plane is ignored. Because of the coupling between the in-plane displacement and the transverse displacement in the geometric equation, it is necessary to introduce a compatible equation to form a complete set of equations. The compatible equation of the double-curved shallow shell under plane stress is expressed as

There are six independent components of internal force in the theory of shallow shell vibration: , , , , , and . After the shear force and was eliminated, three dynamic equilibrium equations can be obtained:where represents the external excitation received by the shell and represents the damping coefficient.

The three components of internal force , , and can be connected by the Airy stress function according to equation (3) in the following way:

By substituting equation (1) into equation (2), the first two dynamic equilibrium equations (3a) and (3b) can be automatically satisfied. The expressions of the internal force and internal moment of the antisymmetric cross-ply laminated composite material can be written as

It can be seen from the above equation that the antisymmetric cross-ply laminated bistable structure has the coupling between stretching and bending. Combining equation (4a4b4c) with equation (5), the expression of the middle plane strain expressed by can be obtained:where , , , and .

By substituting simplified equation (5) of the internal force and internal moment on the cross section of the antisymmetric cross-ply thin shell and midsurface strain expression (6a)–(6c) into dynamic equilibrium equation (3c), we have

Then, by substituting simplified equation (5) of internal force and internal moment and the expression of midsurface strain (6a6b6c) into compatible equation (2), we have

In the discussion of the mode superposition method, the solution of a continuous system can be written as a linear combination of all mode functions. Since the boundary condition of the model is that the four sides are free, the center is fixed, the selected shape function does not need to satisfy any geometric boundary conditions, and the displacement shape function in the direction of is given bywhere represents generalized coordinates and is the mode function of the system, where the subscript and , and all the shape function terms are the combined operations of nonzero trigonometric functions, where , and is the number of degrees of freedom selected in the approximation. Since the model built is a shallow shell structure and the area defined by the transverse displacement and stress function is the same, it is assumed that they can be expanded with the same shape function . Take the mode shape function of the transverse displacement and stress function , as shown below:where and are, respectively, the generalized coordinates of transverse displacement and stress functions . First, substituting equation (10a10b) into dynamic equilibrium equation (7) leads to

Then, equation (10a10b) was substituted into compatible equation (8), and we have

3. Galerkin Discretization

Then, the two equations were discretized by Galerkin so that the infinite degree of the freedom system can be discretized into a finite degree of the freedom system. Take , . Both sides of equilibrium equation (11) were multiplied by , and they were integrated over the area to obtain

Both sides were divided by , and the equilibrium equation obtained is as follows:where , , and . Multiply both sides of compatible equation (12) by and integrate over the area, and we have

Let the coefficient ahead of be 1, and the compatible equation obtained iswhere , , , , and

Due to the decoupling between the generalized coordinates of the stress function and the generalized coordinates of the transverse displacement, the equation set can be decoupled and the governing motion equation of the transverse displacement can be given. By substituting equation (16) into equation (14), the final vibration differential equation can be obtained:where .

Since low-order mode shapes play a major role in vibration, the higher the order, the smaller the role. To facilitate writing, the displacement in the deflection direction of the plate was rewritten as , and the load excitation was replaced by . Therefore, the first three modes were selected to obtain the three-degree-of-freedom ordinary differential governing equations:where is the damping parameter of the system, is the one of the natural frequencies of the system, and is the coefficient before each quadratic nonlinear term.

4. Linear System

4.1. Numerical Analysis of Vibration Characteristics

After obtaining the three-degree-of-freedom ordinary differential equation of the antisymmetric cross-ply bistable thin shell, Figure 3 shows the mode shapes of the first 9 antisymmetric cross-ply bistable shallow shell models obtained by the finite element numerical simulation software ABAQUS. It can be seen from the figure that each is symmetric or antisymmetric with respect to the axis and the axis. The first mode shape is the biaxially antisymmetrical torsion, the second mode shape is the biaxially symmetrical bending vibration, and the third mode shape is the torsional vibration symmetric with respect to -axis and antisymmetric with respect to -axis. Starting from the fourth order, the mode shapes are all bending-torsional coupling vibrations. The fourth mode shape is antisymmetric with respect to -axis and symmetric with respect to -axis, the fifth mode shape is symmetric with respect to the two axes, and the sixth mode shape is antisymmetric with respect to the axis and symmetric with respect to axis. The seventh mode shape is antisymmetric with respect to the two axes. The eighth mode shape is symmetric with respect to the axis and antisymmetric with respect to the axis, and the ninth mode shape is antisymmetric with respect to the two axes. Higher-order mode shapes also satisfy this law.

4.2. The Impact of Geometric Parameters on the Frequency and Mode Shapes of the Bistable Cylindrical Shell

Since the laying angle of antisymmetrical cross-ply was selected as the alternately laying mode at and , simulation modeling and analysis were performed on the first 9 groups of the test specimen in Table 2. By changing the shell length , the initial crosssection radius , the initial cross-section center angle , and the number of layers , the impact of different geometric parameters on the frequency and mode shapes of the antisymmetric laminated bistable shell was studied.

Figure 4 shows the relationship between different shell length and the frequency value of the first 14 antisymmetric cross-ply laminated cylindrical shell. The black, red, and blue lines represent the three situations where , , and , respectively. It can be seen from the figure that the frequency value increases with the increase of the order of mode shape and, at the same time, increase with the decrease of the shell length. The frequency values before the 7th mode shape are close, and the values of the 10th and 11th mode shapes are close to each other, while the values of other mode shapes are obviously different. The frequency values of the 8th and 9th mode shapes of each model are slightly different, and the curves in the figure tend to be flat. The frequency values of the 12th and 13th mode shapes also have the same law. The small figure gives the details of the frequency values of the first 5 models with three shell lengths.

Figures 5 to 7 present the stress cloud of the first 9 mode shapes of the models with three shell lengths obtained by simulation. It can be seen that the main mode shape at each order of the models with three shell lengths has changed. In the figure, the distribution law of the stress on the shell in each mode shape can be observed. The magnitude of the stress increases with the changes of blue, green, yellow, and red. When there are two straight opposite sides of the mode shape, the maximum stress is distributed on the two straight sides of the shell, and the distribution of minimum values on the shell surface is close to a regular rectangle; when the four sides of the mode shape are all curved sides, the maximum values of stress are distributed on the four corner points, and the distribution of the minimum values is irregular. Only the second, fourth, and fifth mode shapes remain unchanged, while other mode shapes are different. In the first 6 mode shapes, bending vibration and torsional vibration appear alternately. Starting from the 7th order, all mode shapes are bending-torsion coupling vibrations.

Figure 8 illustrates the impact of different initial cross section radius on the frequency of the antisymmetric cross-ply laminated cylindrical shell. The black, red, and blue lines, respectively, represent the curves of the frequency values at the first 14 orders in the three cases where , , and . It can be seen from the figure that the frequency value increases with the increase of the order of mode shape and, at the same time, increases with the decrease of the initial cross-section radius . The overall trend in the three cases is similar to that in Figure 3, except that the frequency value increases slowly at the 8th, 9th, 12th, and 13th mode shapes, and the frequency values at other orders show a significant increase. The small figure displays the detailed diagram of the frequency values of the first 5 mode shapes. Figures 9 and 10 present the stress cloud of the first 9 mode shapes when and obtained by simulation, and they were compared with the situation where in Figure 4. In the figure, the stress distribution on the shell in each mode shape can be found, and the stress distribution law is the same as that shown in Figures 57. Only the 4th, 5th, and 8th mode shapes remain unchanged, while other mode shapes are different. In the first 6 mode shapes, bending vibration and torsional vibration appear alternately. Starting from the 7th mode shape, all modes are bending-torsion coupling vibrations.

Figure 11 exhibits the relationship between different initial cross-section central angle and the frequency value . The black, red, and blue lines represent the curves of frequency at the first 14 orders when the central angle , , and , respectively. The small figure is the detailed diagram of the frequency values of the 6th to 9th mode shapes. It can be seen that the frequency values of the 6th to 9th mode shapes when are relatively close, and the curves tend to be flat. In other cases, the frequency value increases with the decrease of the central angle. Figures 12 and 13 show the stress cloud of the first 9 mode shapes under the initial cross-section center angle and obtained by simulation, respectively. The distribution rules of the maximum and minimum stress on the shell surface are the same as those described above. Only the first, second, fourth, and eighth mode shapes remain unchanged, while other mode shapes are quite different. In the first 6 mode shapes, bending vibration and torsional vibration appear alternately. Starting from the 7th mode shape, all mode shapes are bending and torsion coupling vibrations.

Figure 14 describes the impact of different layer thicknesses on the frequency of antisymmetric cross-ply laminated cylindrical shell. The black, red, and blue lines represent the curves of frequency at the first 14 orders when the layer thickness , , and . It can be found from the figure that the frequency value increases with the increase of the order of mode and, at the same time, increases with the increase of the layer thickness. The overall trend in these three cases is similar to that in the previous cases, except that the frequency values of the 3rd, 5th, and 9th mode shapes are significantly increased, and the values at other orders all increase slowly. The small figure gives the detailed diagram of the frequency values of first 5 mode shapes. Figures 15 and 16 are the stress cloud of the first 9 mode shapes when and obtained by the simulation and were compared with the situation in Figure 5 where . In the figure, the stress distribution on the shell in each mode shape can be seen, and the stress distribution law is the same as that shown in Figures 57. Only the first, second, fourth, and eighth mode shapes remain unchanged, while other mode shapes are different. Bending vibration and torsional vibration appear alternately in the first 6 mode shapes. Starting from the 7th mode shape, the mode shapes are all bending-torsion coupling vibrations.

5. Nonlinear System

5.1. Perturbation Analysis When the Main Resonance
5.1.1. Noninternal Resonance Situation

According to the frequency of the linear system, there is a 1 : 2 : 3 internal resonance relationship between the natural frequencies of the first three orders in the bistable system. On the premise that the main resonance is close to , internal resonance situation and noninternal resonance situation were discussed separately.

First, regarding noninternal resonance situation, the solvability condition (the condition for eliminating the long-term term) iswhere is the derivative of about and is the tuning parameter.

The steady-state response that can be obtained iswhere . Therefore, in the case where there is no internal resonance, any nonlinear term will not produce a long-term term, and its first approximation is not affected by the nonlinear term; thus, it is actually the solution to the corresponding linear problem.

5.1.2. Internal Resonance Situation

When there is a 1 : 2 : 3 internal resonance, let , , and , where , , and are the tuning parameters, and the solvability condition is

At this time, the symbol was introduced, let , and the result was divided into a real part and an imaginary part, and we havewhere

Equations (22a)–(22f) shows that the phase modulation between response and excitation begins to change, which is not available in linear systems. The aforementioned phase modulation change is important since it limits the amplitude of the response to a finite value. and were eliminated to obtain the average equation under polar coordinates:

Let , , and . The average equation under rectangular coordinates is as follows:

5.2. Dynamic Analysis of Bifurcation and Chaos When the Main Resonance

Based on rectangular coordinate average equations (25a)–(25f) of the bistable system under 1 : 2 : 3 internal resonance and main resonance, the nonlinear dynamic behavior characteristics of the bistable system with the change of external excitation were analyzed by numerical methods. First, the initial conditions were defined as , , , , , and , and other parameters were set to be , , , , , , , and . In rectangular coordinate average equation (25a), f1 is the parameter related to the external excitation. Here, f1 was selected as the control parameter to study the nonlinear dynamic characteristics of the system affected by the external excitation.

Figure 17 is the bifurcation diagram of the second and third mode shapes of the antisymmetric cross-ply bistable system. The abscissa axis represents the changes of related to the external excitation, and the ordinate axis represents the changes of and related to the second-order and third-order amplitudes of the system. As the external excitation expands, the bifurcation diagram of the system shows obvious up-and-down vibrations, indicating that the nonlinear dynamic behavior of the antisymmetric cross-ply laminated bistable cylindrical shell has the basic characteristics of the bistable state. In the interval , the bistable structure is in a period-doubling state. In the interval , the bistable structure is in chaotic motion, with the amplitude constantly expanding.

Figure 18 is the maximum Lyapunov exponent diagram of the antisymmetric cross-ply laminated bistable system. The abscissa axis represents the changes of related to external excitation, and the ordinate axis represents the maximum Lyapunov exponent. It can be seen from the figure that, in the interval , the maximum Lyapunov exponent is less than zero. In the interval , the maximum Lyapunov exponent is always greater than zero, indicating that the bistable structure has been in a state of chaos. Based on the selected value of , the waveform diagram, phase diagram, Poincaré cross-section diagram, and frequency spectrum diagram of the system under different external excitation values are displayed.

Figures 19 to 22 are the phase diagram, time history diagram, Poincaré cross-section diagram, and frequency spectrum diagram of the corresponding system motion when the parameter values of the external excitation are different. In each figure, (a) and (b), respectively, show the three-dimensional phase diagram in the space and the time history diagram in . (c) and (d), respectively, show the three-dimensional phase diagram in space and the second mode shape phase diagram on the two-dimensional plane . (e) and (f), respectively, show the Poincaré cross-section and frequency spectrum on the two-dimensional plane .

When , it can be seen from Figure 19 that a periodic phenomenon occurs in the system. It can be found from the figure that a single closed loop appears in the two-dimensional and three-dimensional phase diagrams of the system, the waveform diagram shows obvious periodic characteristics, a mapping point appears in the Poincaré mapping, and there is a main frequency component in the spectrogram. Moreover, the largest Lyapunov exponent is close to zero. These characteristics indicate that the system has periodic movement at this time.

When , it can be seen from Figure 20 that an almost periodic phenomenon occurs in the system. It can be found from the figure that the two-dimensional and three-dimensional phase diagrams of the system show dense circular orbits, dense mapping points appear in the Poincaré map, the spectrogram is a discrete spectrum, and the maximum Lyapunov exponent is close to zero. These characteristics reveal that the system has almost periodic motion at this time.

Figure 21 shows the motion characteristics of the system when . The phase diagram of the system is messy, dense, and clearly bounded. There are continuous spectra in the spectrogram, countless mapping points appear in the Poincaré map, and the maximum Lyapunov exponent is greater than zero. These features demonstrate that chaotic movement occurs in the system at this time.

Figure 22 shows the motion characteristics of the system when . The phase diagram trajectory of the system is messy, dense, and clearly bounded. The waveform diagram shows that the system motion has random characteristics. Continuous spectra appear in the spectrogram, and there are countless mapping points in the Poincaré map. Meanwhile, the maximum Lyapunov exponent is greater than zero. These features indicate that chaotic motion occurs in the system at this time.

5.3. Perturbation Analysis When the Main Resonance
5.3.1. Noninternal Resonance Situation

First, regarding the situation of noninternal resonance, the solvability condition (the condition for eliminating the long-term term) is

The resulting solution expression is

When and , we have

Therefore, in the case where there is no internal resonance, any nonlinear term will not produce a long-term term, and its first-order approximation is not affected by the nonlinear term; thus, it is actually the solution of the corresponding linear problem.

5.3.2. Internal Resonance Situation

When there is a 1 : 2 : 3 internal resonance, let , , and , where , , and are the tuning parameters, and the solvability condition is

At this time, the symbol was introduced, let , and we have

The result was divided into a real part and an imaginary part:where

The phase modulation between response and excitation begins to change, which is not available in linear systems. The aforementioned phase modulation change is important since it limits the amplitude of the response to a finite value. and were eliminated to obtain the average equation under polar coordinates:

Let , , and . The average equation under rectangular coordinates was obtained, as follows:

5.4. Dynamic Analysis of Bifurcation and Chaos When the Main Resonance

Next, based on rectangular coordinate average equation (34a34b34c34d34e34f) of the bistable system under 1 : 2 : 3 internal resonance and main resonance, numerical methods were used to analyze the nonlinear dynamic behavior characteristics of the bistable system with the change of external excitation. First, the initial conditions were defined as , , , , , and , and the values of other parameters are , , , , , , , and . In rectangular coordinate average equation (34a34b34c34d34e34f), is the parameter related to the external excitation. Here, was selected as the control parameter to study the nonlinear dynamic characteristics of the system affected by the external excitation.

Figure 23(a) is the bifurcation diagram of the first mode shape of the antisymmetric cross-ply bistable system. The abscissa axis represents the changes of related to external excitation, and the ordinate axis represents the changes of related to the first-order amplitude of the system. As the external excitation increases, the bifurcation diagram of the system shows obvious upper and lower areas, which indicates that the nonlinear dynamic behavior of the antisymmetric laminated bistable cylindrical shell has the basic characteristics of the bistable state. In the interval , the bistable structure has been in a monostable state and transforms between two steady states. In the interval , the bistable structure continues to transform between two steady states, with the amplitude continuously expanding. Figure 23(b) is a bifurcation diagram of the second mode shape of the antisymmetric cross-ply laminated bistable system. The abscissa axis represents the changes of related to external excitation, and the ordinate axis represents the changes of related to the second-order amplitude of the system. With the increase of the external excitation, the bifurcation diagram of this system does not have obvious upper and lower areas. In the interval , the bistable structure is always in the chaotic state of the two steady states up and down near the equilibrium point.

Figure 24 is the maximum Lyapunov exponent diagram of the antisymmetric cross-ply laminated bistable system. The abscissa axis represents the changes of related to external excitation, and the ordinate axis represents the maximum Lyapunov exponent. It can be seen from the figure that, in the interval , the maximum Lyapunov exponent is always greater than zero, indicating that the bistable structure has been in a state of chaos. Based on the selected value of , the waveform diagram, phase diagram, Poincaré cross-section diagram, and frequency spectrum diagram of the system are given below.

Figures 25 to 28 are the phase diagram, time history diagram, Poincaré cross-section diagram, and frequency spectrum diagram of the corresponding system motion when the parameter values of external excitation are different. In each figure, (a) and (b), respectively, show the three-dimensional phase diagram in the space and the time history diagram in . (c) and (d), respectively, show the three-dimensional phase diagram in space and the second mode shape phase diagram on the two-dimensional plane . (e) and (f), respectively, show the Poincaré cross-section and spectrogram on a two-dimensional plane .

When , it can be seen from Figure 25 that a periodic phenomenon occurs in the system. There are multiple closed loops in the two-dimensional and three-dimensional phase diagrams of the system, the waveform diagram has obvious periodic characteristics, multiple mapping points appear in the Poincaré map, and there are multiple main frequency components in the spectrogram. In addition, the maximum Lyapunov exponent is close to zero. These characteristics reveal that the system has periodic movement at this time.

When , it can be seen from Figure 26 that an almost periodic phenomenon occurs in the system. The two-dimensional and three-dimensional phase diagrams of the system show dense circular orbits, dense mapping points appear in the Poincaré map, the spectrogram is a discrete spectrum, and the maximum Lyapunov exponent is close to zero. Based on the characteristics shown in the figure, it can be confirmed that the system has almost periodic motion at this time.

Figure 27 shows the motion characteristics of the system when . The phase diagram of the system is messy, dense, and clearly bounded. There are continuous spectra in the spectrogram, countless mapping points appear in the Poincaré map, and the maximum Lyapunov exponent is greater than zero. These features indicate that chaotic movement occurs in the system.

Figure 28 shows the motion characteristics of the system when . The phase diagram of the system is messy, dense, and clearly bounded. The waveform diagram shows that the system motion has random characteristics. Continuous spectra appear in the spectrogram, and there are countless mapping points in the Poincaré map. The maximum Lyapunov exponent exceeds zero. These characteristics indicate that chaotic motion occurs in the system at this time.

6. Conclusion

First, the relevant vibration characteristics of the antisymmetric cross-ply laminated bistable shallow shell were studied. The finite element software ABAQUS was used to perform numerical simulation on the linear model to explore the influence of different geometric parameters on the vibration characteristics such as frequency and main mode shapes. It can be found that the frequency value increases with the decrease of the model size and increases with the increase of the structure thickness and the laying angle, exhibiting stepped growth overall. The main mode shape of models in each group is similar overall, with little change. The distribution of stress on the shell structure is regular. When there are two straight opposite sides of the mode shape, the maximum values of stress are distributed on the two straight sides of the shell, and the distribution of the minimum values on the shell surface is close to a regular rectangle; when all the four sides of the mode shape are curved, the maximum stress is distributed at the four corner points, and the distribution of the minimum values is irregular. Besides, bending vibration and torsional vibration appear alternately in the first 6 mode shapes, and from the 7th order, the mode shapes are all coupled bending and torsion vibrations.

Secondly, the dynamic equilibrium equation and the compatible equation were simultaneously established to build a low-order nonlinear dynamic model for the antisymmetric cross-ply laminated bistable shallow shell. Perturbation analysis was conducted on the three-degree-of-freedom nonlinear ordinary differential equation of the antisymmetric laminated bistable cylindrical shallow shell under steady-state excitation through the multiscale method, and the average equations in the six-dimensional polar coordinate form and the rectangular coordinate form under 1 : 2 : 3 internal resonance and the combination of the two different main resonances were obtained. The impact of excitation amplitude on the bifurcation behavior of the system was studied, and the following conclusions were drawn: first, from the bifurcation diagram of the system, it can be seen that there are typical nonlinear phenomena such as period doubling and chaos in the bistable system. Based on the fixed parameter value, two-dimensional and three-dimensional phase diagrams, time history diagram, frequency spectrum, Poincaré mapping, and maximum Lyapunov exponent were used to analyze the specific motion state of the system and to explain that the external excitation has a significant impact on the dynamic characteristics of the bistable model. In the case where two main resonances exist and have similar parameters, there are huge differences in the variation law of the nonlinear dynamic behavior characteristics of the antisymmetric cross-ply bistable system. With the increase of damping coefficient, tuning parameters, and excitation amplitude, the system presents complex nonlinear dynamic behavior, and several typical phenomena were found, such as the alternate occurrence of period-doubling bifurcation and chaos.

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 regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the support of the National Natural Science Foundation of China (NNSFC) through Grant no. 11832002 and Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (PHRIHLB).