# Nonlinear Vibration Analysis for Stiffened Cylindrical Shells Subjected to Electromagnetic Environment

Authors: Xue-Qin Li; Guang-Chen Bai; Lu-Kai Song; Wei Zhang
Journal: Shock and Vibration
 (2021)

---

## Abstract

The nonlinear vibration behaviors of stiffened cylindrical shells under electromagnetic excitations, transverse excitations, and in-plane excitations are studied for the first time in this paper. Given the first-order shear deformation theory and Hamilton principle, the nonlinear partial differential governing equations of motion are derived with considering the von Karman geometric nonlinearity. By employing the Galerkin discretization procedure, the partial differential equations are diverted to a set of coupled nonlinear ordinary differential equations of motion. Based on the case of 1 : 2 internal resonance and principal resonance-1/2 subharmonic parametric resonance, the multiscale method of perturbation analysis is employed to precisely acquire the four-dimensional nonlinear averaged equations. From the resonant response analysis and nonlinear dynamic simulation, we discovered that the unstable regions of stiffened cylindrical shells can be narrowed by decreasing the external excitation or increasing the magnetic intensity, and their working frequency range can be expanded by reducing the in-plane excitation. Moreover, the different nonlinear dynamic responses of the stiffened cylindrical shell are acquired by controlling stiffener number, stiffener size, and aspect ratio. The presented approach in this paper can provide an eﬃcient analytical framework for nonlinear dynamics theories of stiffened cylindrical shells and will shed light on complex structure design in vibration test engineering.


---

## Body

## 1. Introduction

With the remarkable development of astronautic science and technology (i.e., space station, spacecraft, satellite, and rocket), the advanced and large-thrust electrodynamic vibration shakers are needed to simulate the actual working conditions and test the structural vibration performance [1–3]. As the critical component of electrodynamic vibration shakers, armature structures can generate the wanted vibration waveforms with given vibration signals, and its resonance frequency directly determines the upper limit of the working frequency of vibration shaker [4, 5]. During operations of vibration shaker, due to the large deformation and structural torsion induced by electromagnetic environment, unstable nonlinear dynamic behaviors such as bifurcation and chaos may appear in armature structure, which will produce vibration waveform deviations and considerably reduce the test accuracy. Therefore, effective nonlinear vibration analysis for armature structures to satisfy the test accuracy becomes increasingly important in space engineering. However, for armature structures operating in an electromagnetic environment (as shown in Figure 1), complex geometric configurations not only exist but also are subjected to transverse excitations, in-plane excitations, and electromagnetic excitations [6–8]. These factors significantly increase the analysis complexity of nonlinear vibrations for armature structures. Up to now, few studies have been performed on revealing the nonlinear vibrations of armature structures under electromagnetic environment, which greatly limits the development of high-accuracy vibration test equipment. The realistic demands on improving vibration tests motivate us to start this research on the nonlinear vibrations of armature structures. In this study, by considering the geometric configurations and stiffness distribution, we simplify the armature structure to the cylindrical shell structure with ring stiffeners, to facilitate reveal the dynamic characteristics by employing nonlinear shell theory, Hamilton principle, and perturbation analysis theory. The novelty of the study is to establish the dynamic model of stiffened cylindrical shells subjected to electromagnetic excitations, transverse excitations, and in-plane excitations for the first time, and the current studies can provide a heuristic way for space engineers to realize the nonlinear vibrations of armature structure in a macroscopic view.Figure 1
Armature structure under electromagnetic environment.Different from the regular cylindrical shells, the structural traits of stiffened cylindrical shell are closer to the real engineering structure (i.e., armature structure), and thereby the stiffened cylindrical shell can perform more accurately the nonlinear vibration analysis from geometric modeling perspective. Up till now, stiffened cylindrical shells have been widely used in the civil, mechanical, aerospace, and aeronautic engineering fields [9–13]. The nonlinear shell theories have been studied in the past several decades [14]. Just like Donnell’s shell theory [15, 16] and Sanders’ shell theory [17, 18], Novozhilov theory [19, 20], Koiter theory [21, 22], and Flügge-Lur’e-Byrne theory [23, 24], which are classic shell theories, assure the computing efficiency in dynamic equation establishment for thin cylindrical shells. To describe the nonlinear vibrations of thin shells accurately, the first-order shear deformation theory [25–27] and the third-order shear deformation theory [28–30] are usually applied. Amabili [31] analyzed the nonlinear vibration of cylindrical shells with simply supported boundary and retained all the nonlinear terms of displacement and rotation in plane by using first-order shear deformation theory. Yazdi [32] studied the large amplitude vibration of composite hyperbolic shells considering von Karman’s geometric nonlinearity and first-order shear deformation theory. Karimiasl [33] studied large amplitude vibration behaviors of multiscale doubly curved shells with piezoelectric layer by adopting Reddy’s third-order shear deformation theory. Gu et al. [34] explored the nonlinear dynamics of circular cylindrical shells considering small initial geometric imperfection. Amabili and Reddy [35] improved the shell deformation theory to consider the rotation inertia terms, the shear deformation terms, and the nonlinear terms of the in-plane and lateral nonlinear displacements. According to these investigations, it is found that different nonlinear shell theories have different application scenarios. Moreover, due to the existence of stiffeners, the nonlinear behavior will produce some different characteristics and result in serious consequences, such as the bifurcation point drifts backward, and if the current theories and analysis results of regular cylindrical shells are directly used for stiffened cylindrical shells, unexpected deviations will arise. Therefore, it is significant for us to choose the appropriate nonlinear shell theory according to the geometric configuration of stiffened cylindrical shells.Under the coupling effects of various complex excitations, the nonlinear vibrations and bifurcation phenomenon of stiffened cylindrical shells would be incurred and may lead to significant noises and resonance issues [36–38]. At present, there are more and more researches focusing on the frequency, chaotic motions, and instability for circular cylindrical shells under complex multiple excitations, and fruitful results have been achieved. Goncalves et al. [39] investigated the nonlinear vibration of cylindrical shells subjected to pulsating axial loads by using Donnell’s shadow shell equation. Liu et al. [40] studied the bifurcation, intermittent chaos, and nonlinear vibration of a large deployable space antenna under 1 : 3 internal resonant thermal load. Dogan [41] established nonlinear response model for double-wall Sandwich cylindrical shells under random excitation and studied the parameter sensitivity on the nonlinear responses. Kumar et al. [42] studied the nonlinear response of elliptical shell under transverse harmonic load. Duc et al. [43] investigated the nonlinear response of functionally graded cylindrical shell with a ceramic metal-ceramic layer under uniformly distributed radial loads. Breslavsky et al. [44] studied that the dynamic behavior of cylindrical shells subjected to multiharmonic excitation is particularly complex and exhibits many types of nonlinear behaviors: simple periodic vibration, quasiperiodic oscillations, subharmonic response, period-doubling bifurcations, and chaos.According to these researches on circular cylindrical shells, we find that the main objective of the current researches aims at investigating the regular cylindrical shells subjected to radial or/and in-plain excitations. However, multiscale composites, doubly curved shells and circular cylindrical shells with longitudinal or circumferential stiffeners, are not conventional in engineering practice. The existence of stiffeners will affect the dynamic characteristics of stiffened cylindrical shells and significantly increase the analysis complexity compared to that of ordinary cylindrical shells without stiffeners [45, 46]. Moreover, stiffened cylindrical shells operating in electromagnetic environment usually suffer the coupling effects of transverse excitations, in-plain excitations, and electromagnetic excitations. Without considering complex multiple excitations in one unified framework, it will bring about unacceptable analysis deviations when describing the complex amplitude-frequency characteristics and chaotic motion behaviors [47–49]. Especially subjected to the strong electromagnetic environment, it is evitable to bring in unacceptable analysis deviations if the electromagnetic excitation and various external excitations are not considered simultaneously.In this paper, the nonlinear vibrations of stiffened cylindrical shells subjected to electromagnetic environment are studied from a theoretical analysis view. A dynamic model of stiffened cylindrical shells subjected to transverse excitations, in-plane excitations, and electromagnetic excitations is proposed for the first time. To begin with, in the framework of the first-order shear deformation shell theory, the governing equations of motion for stiffened cylindrical shells are derived by using Hamilton’s principle. Then, the governing equations of motion are discretized into nonlinear ordinary differential governing equations by utilizing the Galerkin discretization procedure. The multiscale method of perturbation analysis is employed to obtain four-dimensional nonlinear averaged equations under the 1 : 2 internal resonance and principal resonance-1/2 subharmonic parametric resonance. Moreover, the nonlinear dynamic behaviors and the jump phenomena are revealed in the amplitude-frequency curves by considering the effects of geometric parameters and multiple excitations. Lastly, the periodic and chaotic motions of the stiffened cylindrical shell are investigated by Runge–Kutta approach, and the influence of stiffener number, stiffener size, and aspect ratio on the nonlinear dynamic responses of the stiffened cylindrical shell is acquired.


## 2. Dynamic Equations of Motion

The simplified geometric model of stiffened cylindrical shell in an external uniform electromagnetic field (0, 0,Bz) is presented, as shown in Figure 2. The basic parameters of stiffened cylindrical shell involve shell length L, shell middle surface radius R, shell thickness h, stiffener height hs, stiffener thickness ds, and stiffener spacing bs. To analyze the vibration characteristics of the stiffened cylindrical shell, the cylindrical coordinate system (x, θ, z) is built on the left-edge middle surface. Herein, x denotes the axis direction of stiffened cylindrical shell, and θ and z represent the circumferential and radial directions of the stiffened cylindrical shell.Figure 2
Simplified geometric model of stiffened cylindrical shell.To describe the vibration behaviors of stiffened cylindrical shells, by considering the transverse shear deformation, Reddy’s first-order shear deformation theory-based [25–27] dynamic model is constructed. Assuming that the transverse normal would not remain normal to the middle surface after shell deformation, the displacement field is established:(1)ux,θ,z,t=u0x,θ,t+zϕxx,θ,t,vx,θ,z,t=v0x,θ,t+zϕθx,θ,t,wx,θ,z,t=w0x,θ,t,where u0, v0, and w0 indicate middle surface displacements along the x-axis, θ-axis, and z-axis, respectively; ϕx and ϕθ are the rotation angle along the θ-axis and x-axis, respectively; t is the time variable. With considering the von Karman geometric nonlinearity rule [50], the geometric equation is built:(2)εxxεθθγxθγxzγθz=εxx0εθθ0γxθ0γxz0γθz0+zεxx1εθθ1γxθ1γxz1γθz1,εxx0εθθ0γxθ0γxz0γθz0=∂u0∂x+12∂w0∂x2∂v0R∂θ+12∂w0R∂θ2+w0R∂u0R∂θ+∂v0∂x+∂w0∂x∂w0R∂θϕx+∂w0∂xϕθ+∂w0R∂θ,εxx1εθθ1γxθ1γxz1γθz1=∂ϕx∂x∂ϕθR∂θ∂ϕxR∂θ+∂ϕθ∂x00,where εxx, εθθ, γxθ, γxz, and γθz denote the strain components; εxx0,εθθ0,γxθ0,γxz0,γθz0 are the membrane strains; εxx1,εθθ1,γxθ1,γxz1,γθz1 are the flexural strains.In light of the general Hook’s law, to further reveal the stress behavior in terms of displacement, the constitutive relation of stiffened cylindrical shells is indicated:(3)σxxσθθτxθτxzτθz=Q11Q12Q12Q22Q44Q55Q66εxxεθθγxθγxzγθz,where Qij (i = 1, 2 and j = 1, 2; i = 4–6 and j = 4–6) are the elastic constants of stiffened cylindrical shell:(4)Q11=Q22=E1−v2,Q12=vE1−v2,Q44=Q55=Q66=E21+v,σθθs=Esεθθ,in which σxx,σθθ,σθθs are the coefficients of the normal stress; τxθ are the coefficients of the shear stress; τxz,τθz are the coefficients of the bending stress; E is Young’s modulus of the cylindrical shell; Es is Young’s modulus of stiffeners.The force and bending moment distributions on middle surface of cylindrical shells are shown in Figure3, and they can be quantified in equation (5).Figure 3
(a) Force distribution and (b) bending moment distribution on middle surface of the cylindrical shell.
(a)(b)(5)Nxx=∫−h/2h/2σxxdz,Mxx=∫−h/2h/2σxxzdz,Nθθ=∫−h/2h/2σθθdz,Nθθs+∫h/2h/2+hsσθθdz,Nxθ=∫−h/2h/2τxθdz,Mθθ=∫−h/2h/2σθθzdz,Mθθs=∫h/2h/2+hsσθθzdz,Mxθ=∫−h/2h/2τxθzdz,Qx=∫−h/2h/2τxzdz,Qθ=∫−h/2h/2τθzdz.As operating in electromagnetic environment, stiffened cylindrical shells make deformation under the electromagnetic force, and the structural deformation further leads to the change of distribution condition of electromagnetic flux and electromagnetic force. These coupling effects between the electromagnetic field and the mechanical field bring in the high nonlinearity and large time-varying characteristics, which leads to great difficulty in performing the quantitative analysis of mechanical behavior for stiffened cylindrical shells. In this study, to clearly describe the mechanical behavior of stiffened cylindrical shells, considering the influence of electromagnetic force on virtual work under external excitations, the dynamic equations of motion under electromagnetic environment are established. By introducing the Lorentz law of force [51], the electromagnetic force acting on stiffened cylindrical shells is presented as follows.Based on the law of electromagnetic induction [52], the electric current induced by cylindrical shell motion in radial electromagnetic field can be expressed as(6)J=σV×B=σ∂u∂t×B,where J denotes the electric current vectors, σ denotes the conductivity, V denotes the speed, and B denotes the electromagnetic intensity, where B = (0, 0, Bz). Therein, i, j, and k are regarded as the vector in the three directions of x, θ, and z; then, the displacement vector of the stiffened cylindrical shell u can be expressed as(7)u=u0+zϕxi+v0+zϕθj+w0x,θ,tk.In the constant electromagnetic field, the Lorentz electromagnetic force per unit volume can be expressed as(8)ffu,fv,fw=J×B=σ∂u∂t×B×B=σ∂u0∂t+z∂ϕx∂ti,∂v0∂t+z∂ϕθ∂tj,∂w0∂tk×0,0,Bzk×0,0,Bzk=−σBz2∂u0∂t+z∂ϕx∂ti,−σBz2∂v0∂t+z∂ϕθ∂tj,0=fui+fvj.Furtherly, the Lorentz forcesFx and Fθ (in x and θ-axis) and electromagnetic moments Mx and Mθ (in θ and x-axis) can be derived as follows:(9)Fx=∫−h/2h/2fudz=−σBz2∂u0∂th,Fθ=∫−h/2h/2fvdz=−σBz2∂v0∂th,(10)Mx=∫−h/2h/2fuzdz=σh3Bz212∂ϕx∂t,Mθ=∫−h/2h/2fvzdz=σh3Bz212∂ϕθ∂t.Based on the nonlinear constitutive equation and the magnetic force described by the Lorentz law of force, the nonlinear dynamic equation of the system can be established. Given Hamilton’s principle [53], the nonlinear dynamic equation of stiffened cylindrical shell under electromagnetic environment is established as follows:(11)∫t2t1δT−U+δWdt=0.Therein, the kinetic energyδT, potential energy δU, and virtual work δW under external excitations can be signified as follows:(12a)δT=∫Vρu˙δu˙,v˙δv˙,w˙δw˙dV+∫Vsdsbsρu˙δu˙,v˙δv˙,w˙δw˙dVs,(12b)δU=∫Vσxxδεxx+σθθδεθθ+σxθδεxθ+σxzδεxz+σθzδεθzdV+∫VsdsbsσθθδεθθdVs,(12c)δW=∫VFxδu0+Fθδv0+Mxδϕx+MθδϕθdV+∫ANw0,xxδw0dA+∫AFδw0dA,where N = N1sin (ω1t) represents the in-plane excitation, N1 represents the in-plane excitation amplitude, F = F1sin (ω2t) represents the time-depend transverse excitation, F1 represents the transverse excitation amplitude, ω1 represents the in-plane excitation frequency, and ω2 represents the transverse excitation frequency.The nonlinear partial differential governing equations of motion for stiffened cylindrical shell are given by substituting the kinetic energy, potential energy, and virtual work into Hamilton’s equation:(13)δu0:Nxx,x+Nxθ,θR=I0u¨0+I1ϕ¨x+σBz2u˙0h,δv0:Nθθ,θR+Nxθ,x+Nθθ,θs=I1ϕ¨θ+I0v¨0+σBz2v˙0h,δw0:Nxx,x∂w0∂x+Nθθ,θs∂w0∂x+Nxx∂2w0∂x2+Nθθs∂2w0∂x2+Nθθ,θR∂w0R∂θ+Nθθ∂2w0R2∂θ2−NθθR+Nxθ,x∂w0R∂θ+Nxθ∂2w0R∂θ∂x+Nxθ,θR∂w0∂x+Nxθ∂2w0R∂θ∂x+Qx,x+Qθ,θR−N∂2w0∂x2+F=I0w¨0,δϕx:Mxθ,θR+Mxx,x−Qx=I2ϕ¨x+I1u¨0−σh3Bz212ϕ˙x,δϕθ:Mθθ,θR+Mxθ,x+Mθθ,θs−Qθ=I1v¨0+I2ϕ¨θ−σh3Bz212ϕ˙θ.By substituting the force and moment resultants (5) into equation (13), the governing equations of motion in the generalized displacements (u0, v0, and w0) can be rewritten as(14a)A11∂2u0∂x2+A11∂w0∂x∂2w0∂x2+A12∂2v0R∂θ∂x+A12∂w0R∂θ∂2w0R∂θ∂x+A12∂w0R∂x+B11∂2ϕx∂x2+Bs∂2ϕx∂x2+As∂2u0∂x2+As∂w0∂x∂2w0∂x2+B12∂2ϕθR∂θ∂x+A33∂2u0R2∂θ2+A33∂2v0R∂θ∂x+A33∂2w0R∂θ∂x∂w0R∂θ+A33∂w0∂x∂2w0R2∂θ2+B33∂2ϕxR2∂θ2+B33∂2ϕθR∂θ∂x=I0u¨0+I1ϕ¨x,(14b)A12∂2u0R∂x∂θ+A12∂w0∂x∂2w0R∂θ∂x+A22∂2v0R2∂θ2+A22∂w0R∂θ∂2w0R2∂θ2+A22∂w0R2∂θ+B12∂2ϕxR∂θ∂x+B22∂2ϕθR2∂θ2+A33∂2u0R∂θ∂x+A33∂2v0∂x2+A33∂2w0∂x2∂w0R∂θ+A33∂w0∂x∂2w0R∂θ∂x+B33∂2ϕxR∂θ∂x+B33∂2ϕθ∂x2=I1ϕ¨θ+I0v¨0,(14c)As∂2u0∂x2∂w0∂x+As∂w0∂x∂2w0∂x2∂w0∂x+Bs∂2ϕx∂x2∂w0∂x+As∂u0∂x∂2w0∂x2+As12∂w0∂x2∂2w0∂x2+Bs∂ϕx∂x∂2w0∂x2∂w0∂x+A11∂2u0∂x2+A11∂w0∂x∂2w0∂x2+A12∂2v0R∂θ∂x+A12∂w0R∂θ∂2w0R∂θ∂x+A12∂w0R∂x+B11∂2ϕx∂x2+B12∂2ϕθR∂θ∂x+A11∂u0∂x+A1112∂w0∂x2+A12∂v0R∂θ+A1212∂w0R∂θ2+A12w0R+B11∂ϕx∂x+B12∂ϕθR∂θ∂2w0∂x2+A12∂2u0R∂x∂θ+A12∂w0∂x∂2w0R∂θ∂x+A22∂2v0R2∂θ2+A22∂w0R∂θ∂2w0R2∂θ2+A22∂w0R2∂θ+B12∂2ϕxR∂θ∂x+B22∂2ϕθR2∂θ2∂w0R∂θ+A12∂u0∂x+A1212∂w0∂x2+A22∂v0R∂θ+A2212∂w0R∂θ2+A22w0R+B12∂ϕx∂x+B22∂ϕθR∂θ∂2w0R2∂θ2−1RA12∂u0∂x+A1212∂w0∂x2+A22∂v0R∂θ+A2212∂w0R∂θ2+A22w0R+B12∂ϕx∂x+B22∂ϕθR∂θ+A33∂2u0R∂θ∂x+A33∂2v0∂x2+A33∂2w0∂x2∂w0R∂θ+A33∂w0∂x∂2w0R∂θ∂x+B33∂2ϕxR∂θ∂x+B33∂2ϕθ∂x2∂w0R∂θ+A33∂u0R∂θ+A33∂v0∂x+A33∂w0∂x∂w0R∂θ+B33∂ϕxR∂θ+B33∂ϕθ∂x∂2w0R∂θ∂x+A33∂2u0R2∂θ2+A33∂2v0R∂x∂θ+A33∂2w0∂xR∂θ∂w0R∂θ+A33∂w0∂x∂2w0R2∂θ2+B33∂2ϕxR2∂θ2+B33∂2ϕθR∂x∂θ∂w0∂x+A33∂u0R∂θ+A33∂v0∂x+A33∂w0∂x∂w0R∂θ+B33∂ϕxR∂θ+B33∂ϕθ∂x∂2w0R∂θ∂x+A44∂ϕx∂x+A44∂2w0∂x2+A55∂ϕθR∂θ+A55∂2w0R2∂θ2+Bz2h2μ0∂w0∂x2+∂w0R∂θ2+ρgh∂w0∂x−ρghL−x∂2w0∂x2−N∂2w0∂x2+F=I0w¨0(14d)Bs∂2u0∂x2+Bs∂w0∂x∂2w0∂x2+Ds∂2ϕx∂x2+B33∂2u0R2∂θ2+B33∂2v0R∂θ∂x+B33∂2w0R∂θ∂x∂w0R∂θ+B33∂w0∂x∂2w0R2∂θ2+D33∂2ϕxR2∂θ2+D33∂2ϕθR∂θ∂x+B11∂2u0∂x2+B11∂w0∂x∂2w0∂x2+B12∂2v0R∂θ∂x+B12∂w0R∂θ∂2w0R∂θ∂x+B12∂w0R∂x+D11∂2ϕx∂x2+D12∂2ϕθR∂θ∂x−A44ϕx−A44∂w0∂x+σh3Bz212∂u0∂t+z∂ϕx∂t=I2ϕ¨x−I1u¨0,(14e)B12∂2u0R∂x∂θ+B12∂w0∂x∂2w0R∂θ∂x+B22∂2v0R2∂θ2+B22∂w0R∂θ∂2w0R2∂θ2+B22∂w0R2∂θ+D12∂2ϕxR∂x∂θ+D22∂2ϕθR2∂θ2+B33∂2u0R∂θ∂x+B33∂2v0∂x2+B33∂2w0∂x2∂w0R∂θ+B33∂w0∂x∂2w0R∂θ∂x+D33∂2ϕxR∂θ∂x+D33∂2ϕθ∂x2−A55ϕθ−A55∂w0R∂θ+σh3Bz212∂v0∂t+z∂ϕθ∂t=−I1v¨0+I2ϕ¨θ.Note that the parts which involve magnetic field strengthBZ are used to describe the contribution of the electromagnetic environment. Moreover, the coefficients Aij, Bij, Dij, As, Bs, and Ds can be further expressed as follows:(15a)Aij,Bij,Dij=∫−h/2h/2Qij1,z,z2dz,(15b)As,Bs,Ds=∫h/2h/2+hsEsbsds1,z,z2dz.


## 3. Equations’ Discretization and Perturbation Analysis

In this section, a two-step solution procedure is considered to solve the nonlinear dynamic equations of the stiffened cylindrical shell: (i) the Galerkin discretization approach is applied to divert the partial differential equations into ordinary differential equations; (ii) the multiscale method of perturbation analysis is adopted to obtain the four-dimensional average equation of the system. The solution procedures of the nonlinear dynamic equations are summarized as follows.Considering the deformation features of armature structures, simply supported boundaries are employed to constrain both ends of the stiffened cylindrical shell, which can be expressed as(16)v0=w0=ϕθ=Mxx=0,∫02πRNxx|x=0,Ldθ=∫02πRNdθ,atx=0,x=L,where N0 represents the in-plane excitation.Regarding the fact that in-plane inertial terms are much smaller than the transverse inertial terms [54], the ordinary differential governing equations are derived by ignoring the in-plane inertial term in governing equations. Herein, the continuous system of the partial governing differential equation is truncated into a two-degree-of-freedom system of ordinary differential governing equations by using the Galerkin discretization approach, and the reasonable approximation functions are desired to expand the displacements in the middle surface. Since most of the vibration energy mainly concentrates in the first-order and second-order modes, accounting for 90% or more of the system energy [55], the first two modes are taken into account to describe the nonlinear dynamic behaviors of stiffened cylindrical shells. Based on [56, 57], the displacements u0, v0, w0, ϕx, and ϕθ which satisfy the boundary conditions are represented as(17)u0=U1tcosπxLcos3θ+U2tcos3πxLcosθ,v0=V1tsinπxLcos3θ+V2tsin3πxLcosθ,w0=W1tsinπxLsin3θ+W2tsin3πxLsinθ,ϕx=Ψx1tcosπxLsin3θ+Ψx2tcos3πxLsinθ,ϕθ=Ψθ1tsinπxLcos3θ+Ψθ2tsin3πxLcosθ,where W1 and W2 (t) denote the time-dependent total amplitude; U1 (t), U2 (t), V1 (t), V2 (t), Ψx1 (t), Ψx2 (t), Ψθ1 (t), and Ψθ2 (t) denote the time-dependent amplitude functions.Considering the evenly distributed periodic vibration load of the armature structure in engineering practice, to simulate the real loading conditions, the transverse uniformly distributed harmonic excitation is applied to the stiffened cylindrical shell. The uniformly distributed transverse harmonic excitation can be indicated as(18)F=F1sinπxLcos3θ+F2sin3πxLcosθ,where F1 and F2 represent the amplitudes of the transverse excitation.Considering that the transverse nonlinear vibrations are the primary motion type for stiffened cylindrical shells, the inertia terms in equation (14) are removed in partial differential governing equations. By substituting equations (17) and (18) into (14) and applying the Galerkin discretization procedure, the equations in u0, v0, ϕx, and ϕθ directions in displacements equation are established simultaneously, in which the generalized coordinates U1 (t), U2 (t), V1 (t), V2 (t), Ψx1 (t), Ψx2 (t), Ψθ1 (t), and Ψθ2 (t) can be expressed by the generalized coordinates W1 and W2 and brought into the w0 direction, and the two-degree-of-freedom discrete ordinary differential governing equations of transverse motion for the stiffened cylindrical shell can be driven as(19)W1..+ζ12W1.+ω12W1+ζ14W13+ζ15W1W22+ζ16N1W1cosQt=ζ17F1cosΩt,W2..+ζ22W2.+ω22W2+ζ24W12W2+ζ25W23+ζ26N1W2cosQt=ζ27F2cosΩt,where W = [W1W2]T represent the transverse displacements vectors; ζij (i = 1, 2, j = 2–7) represent the coefficients; Q and Ω represent the excitation frequencies; point above variables represents the time derivative; N1 represents the amplitude of in-plane excitation. Note that, even for the first two modes, the established governing equation also involves quadratic terms, cubic terms, and parametric and transverse excitations, which illustrates the high nonlinearity of transverse vibrations of stiffened cylindrical shells.For the high-dimensional nonlinear dynamic systems, when there exists a special internal resonant relationship between two linear natural frequencies, the large amplitude nonlinear responses may suddenly happen owing to the modal interactions [58]. In this section, to reveal the nonlinear vibration behaviors of stiffened cylindrical shells, the perturbation analysis for the system’s primary resonance is implemented with the multiscale method. By considering the case of 1 : 2 internal resonance and primary resonance-1/2 subharmonic parametric resonance, the resonant relations of stiffened cylindrical shells are given as follows:(20)2ω1=ω−εσ1,ω2=ω−εσ2,where ω = Q = Ω, ε is the small perturbation parameter, and σ1 and σ2 are the two detuning parameters.To obtain the dimensionless ordinary differential governing equations for the stiffened cylindrical shell, the transformations of the variables and the parameters are introduced as follows:(21)τ=tω1,W1=f1h,W2=f2h,Ω¯=Ωω1,Q¯=Qω1,ζ¯12=ζ12ω1,ζ¯14=ζ14h2ω12,ζ¯15=ζ15h2ω12,ζ¯16=ζ16ω12,ζ¯17=ζ17ω12h,ζ¯22=ζ22ω1,ζ¯24=ζ24h2ω12,ζ¯25=ζ25h2ω12,ζ¯26=ζ26ω12,ζ¯27=ζ27ω12h,ω¯1=1,ω¯2=ω2ω1,=2,Substituting equation (21) into equation (19), the dimensionless equations are acquired as follows:(22a)∂2f1∂τ2+ζ¯12∂f1∂τ+ω¯12f1+ζ¯14f13+ζ¯15f1f22+ζ¯16N1cosQ¯τf1=ζ¯17F1cosΩ¯τ,(22b)∂2f2∂τ2+ζ¯22∂f2∂τ+ω¯22f2+ζ¯24f12f2+ζ¯25f23+ζ¯26N1cosQ¯τf2=ζ¯27F2cosΩ¯τ.To precisely obtain the uniform solutions of equations (22a) and (22b), the multiscale method [54] is employed in the following form:(23)fnt,ε=un0T0,T1,T2+εun1T0,T1,T2+ε2un2T0,T1,T2+⋯,n=1,2,where time scale variables T0 = τ, T1 = ετ, and T2 = ε2τ.To express the equation conveniently, the time derivatives used in the multiscale method are given as follows:(24a)ddt=∂∂T0∂T0∂t+∂∂T1∂T1∂t+⋯=D0+εD1+⋯,(24b)d2dt2=D0+εD1+⋯2=D02+2εD0D1+⋯,where the partial derivative operators D0=∂/∂T0andD1=∂/∂T1.Substituting equations (23) and (24b, 24a) into (22a) and (22b) and balancing the coefficients of the like power of ε yield the following differential equations:(25a)ε0orderD02u10+ω12u10=0,(25b)D02u20+ω22u20=0,(25c)ε1orderD02u11+ω12u11=ζ¯12D0u10+2D0D1u10+ζ¯14u103+ζ¯15u10u202+ζ¯16N1u1012e−iQ¯T0+12eiQ¯T0−ζ¯17F112e−iT0Ω¯+12eiT0Ω¯,(25d)D02u21+ω22u21=ζ¯22D0u20+2D0D1u20+ζ¯24u102u20+ζ¯25u203+ζ¯26N1u2012e−iQ¯T0+12eiQ¯T0−ζ¯27F212e−iT0Ω¯+12eiT0Ω¯.Since equations (25a) and (25b) are linear ordinary differential equations, the complex form solutions can be expressed as follows:(26a)u10=A1T1eiω1T0+e−iω1T0A1¯T1,(26b)u20=A2T1eiω2T0+e−iω2T0A2¯T1,where A1¯ and A2¯ are the complex conjugate of A1 and A2. Substituting equations (26a) and (26b) into (25c) and (25d) yield(27a)D02u11+ω12u11=eiω1T02ζ15A1T1A2T1A2¯T1+12ζ16N1A1¯T1eiεσ1T0+iεσ2T0+3ζ14A1T12A1¯T1+iζ12ω1A1T1+2iD1A1T1ω1+cc+NST,(27b)D02u21+ω22u21=eiω2T02ζ24A1T1A2T1A1¯T1+3ζ25A2T12A2¯T1+iζ22ω2A2T1−12ζ27F2eiεσ2T0+2iD1A2T1ω2+cc+NST,where cc and NST denote the complex conjugate part on the right side of equation (27) and the terms that do not produce secular terms, respectively. Note that the variable symbol ζ¯ is rewritten as ζ to facilitate the following formula derivation.To write the modulation equations in polar coordinates, the amplitude functions are set as functionsA1 and A2, which can be denoted in the exponential form as follows:(28a)A1T1=12a1eiϕ1,(28b)A2T1=12a2eiϕ2.Substituting equations (28a) and (28b) into (27a) and (27b) and separating the real and imaginary parts, four-dimensional averaged equations in the polar form are obtained:(29a)a˙1=14ζ16a1N1sinϕ1+12ζ12a1ω1,(29b)a˙2=12ζ22a2ω2−12ζ27F2sinϕ2,(29c)ϕ˙1=14ζ15a1a22+14ζ16a1N1cosϕ1−12a1ω1σ1+σ2+3ζ14a138,(29d)ϕ˙2=14ζ24a12a2+3ζ25a238−12ζ27F2cosϕ2−a2T1σ2ω2.Ultimately, to obtain the numerical solutions effectively, the complex partial differential kinetic equations (i.e., equation (14)) are simplified into simple four-dimensional average equations (i.e., equation (29)).


## 4. Analysis of Resonant Responses

In this section, based on the established four-dimensional averaged equations (equations (29a)–(29d)) in the polar form, by considering the electromagnetic excitations, transverse excitations, and in-plane excitations, the amplitude-frequency characteristics of stiffened cylindrical shells are investigated in this section. It should be noted that the established four-dimensional averaged equations are derived and solved on MATHEMATICA and MATLAB, respectively. The basic parameters of the stiffened cylindrical shell are set as follows: elastic modulus E = 2.1 × 1011 Pa, density ρ = 7800 kg/m3, Poisson’s ratio ν = 0.3, conductivity σ = 2.3 × 106 (Ω·m)−1, L = 0.485 m, R = 0.1665 m, h = 0.007 m, and in-plane and transverse excitation frequencies are Q = Ω = 150 Hz [26]. Moreover, by setting the left parts of equations (equations (29a)–(29d)) to equal zero, the effects of different parameters on the resonant responses of the stiffened cylindrical shell are studied, and the nonlinear dynamic responses are described in amplitude-frequency curves under the cases of the primary resonance-1/2 subharmonic parametric resonance and 1 : 2 internal resonance, as shown in Figures 4–7. Therein, it is worth noting that the dotted line represents the unstable region, the solid line indicates the stable area, and the red dot represents the limit point in these amplitude-frequency curves.Figure 4
Amplitude-frequency curves with different transverse excitations.Figure 5
Amplitude-frequency curves with different in-plane excitations.Figure 6
Amplitude-frequency curves with different electromagnetic intensities.Figure 7
Amplitude-frequency curves with different small tuning parameters.Figure4 shows the influence of different transverse excitation amplitude on the amplitude-frequency curve. The simulation results reveal that the amplitude of the first-order mode is larger than that of the second-order mode. Both the first two modes make left bending formants, so it can be judged that the system makes nonlinear softening characteristics. Moreover, with the increase of transverse excitation amplitude F2, the amplitudes and resonant interval become larger and the unstable region increases gradually. These changes indicate that the stronger resonance responses are excited since more energy is absorbed with the increase of the transverse excitation. Furthermore, the resonance regions of the first-order modes are consistent with that of the second-order modes, which illustrates that the resonance of the first-order modes is caused by the main resonance phenomenon of the second-order modes. Besides, the bending of the amplitude-frequency curves also leads to the multiple amplitudes and the jump phenomena. Herein, the multiple amplitudes mean that there are many possible responses in stiffened cylindrical shells, and the actual responses are determined by the initial conditions; the jump phenomenon in stiffened cylindrical shells means that the amplitude changes abruptly.To investigate the effects of in-plane excitation amplitudeN1 on the amplitude-frequency response, the amplitude variation of the modes is analyzed, and the results are shown in Figure 5. In Figure 5, the upper and lower four curves represent the amplitude-frequency curve of the first two modes when N1 = 3.2 × 107 N/m2, N1 = 3.5 × 107 N/m2, N1 = 3.8 × 107 N/m2, and N1 = 4.1 × 107 N/m2, respectively. It can be seen from the figure that, with the increase of tuning parameter σ1, there will be a jump of amplitude in the first two modal amplitudes in an unstable region. Moreover, we also find that the first-order and the second-order modal limit points are σ1 = 0.2021 and σ1 = 0.2818 when N1 = 3.2 × 107 N/m2; σ1 = 0.2402 and σ1 = 0.3194 when N1 = 3.5 × 107 N/m2; σ1 = 0.2769 and σ1 = 0.3570 when N1 = 3.8 × 107 N/m2; σ1 = 0.3145 and σ1 = 0.3945 when N1 = 4.1 × 107 N/m2. The results show that, with the increase of in-plane excitation amplitude N1, the unstable region gradually moves to the right, which means that the natural frequency of the system is increasing. Besides, with the increase of N1, the ordinates of limit points of the first-order modes move upward, while the ordinates of limit points of the second-order modes remain unchanged, indicating that the increase of N1 increases the energy of the first-order modes but has little effect on the amplitude of the second-order modes. These results illustrate that the subharmonic parametric resonance which is forced by axial in-plane excitation occurs in the first modes.To verify the effects of electromagnetic intensitiesB on the amplitude-frequency response of the system, the amplitude changes of the modes under different electromagnetic intensities are analyzed. By setting transverse excitation amplitude F2 = 4.5 × 105 N/m2, in-plane excitation amplitude N1 = 3.2 × 107 N/m2, and the amplitude-frequency curves are drawn in Figure 6. Note that the upper and lower four curves represent the amplitude-frequency response curve of the first two modes when B = (18–21) T, respectively. From the figure, we find that limit points of first two modal tuning parameters are σ1 = 0.2021 and σ1 = 0.2818 when B = 18 T; σ1 = 0.2402 and σ1 = 0.3194 when B = 19 T; σ1 = 0.2769 and σ1 = 0.3570 when B = 20 T; σ1 = 0.3145 and σ1 = 0.3945 when B = 21 T. Besides, with the increase of the electromagnetic intensity, the energy of the first two modes is decreasing gradually. Moreover, when passing through the resonance frequency, a certain degree of the mutation of the first two modes is discovered.Figure7 shows the effect of the small tuning parameter σ2 on the amplitude-frequency responses of the system, in which σ2 is used to adjust the second-order frequency. The upper and lower four curves represent the amplitude-frequency curve of the first two modes when σ2 = 0; σ2 = 0.05 and σ2 = 0.10; σ2 = 0.15, respectively. The two limit points of the first two modal tuning parameters are σ1 = 0.2021 and σ1 = 0.2818 when σ2 = 0; σ1 = 0.2658 and σ1 = 0.3328 when σ2 = 0.05; σ1 = 0.3220 and σ1 = 0.3784 when σ2 = 0.10; σ1 = 0.3721 and σ1 = 0.4193 when σ2 = 0.15. It can be found that when increasing the small tuning parameter σ2, the unstable region gradually shifts to the right. Meanwhile, with the increase of σ2, the amplitudes and resonant intervals become greater in amplitude-frequency curves. Besides, the amplitude of the first-order modes uniformly moves up and the ordinates of limit points move down, while the ordinates of limit points of the second-order modes remain unchanged. It is indicated that, with the increase of small tuning parameter σ2, the energy of the first-order modes will increase obviously, while the amplitude of the second-order modes will not be affected basically.


## 5. Numerical Simulation of Periodic and Chaotic Dynamics

In this section, with consideration of the influence of stiffener number, stiffener size, and aspect ratio, the periodic and the chaotic motions of stiffened cylindrical shells in the case of 1 : 2 internal resonance and primary resonance-1/2 subharmonic parametric resonance are studied in the following four cases. Therein, Case 1 is used to verify the feasibility of the proposed method, and its solutions are regarded as the benchmark results of the following cases; Cases 2–4 are applied to reveal the influence of stiffeners on the nonlinear vibration behaviors by adjusting the parameters (i.e., stiffener number, stiffener height, and aspect ratio). The change of parameters in Cases 1–4 is illustrated in Table1.Table 1
The change of parameters in Cases 1–4.
ParametersStiffener numberStiffener height (m)Aspect ratioOthersCase 150.0032SameCase 2100.0032SameCase 350.0062SameCase 450.0033.5Same


### 5.1. Case 1

To describe the nonlinear vibration responses of stiffened cylindrical shells, the basic parameters of the stiffened cylindrical shell are set as follows: elastic modulusE = 2.1 × 1011 Pa, density ρ = 7800 kg/m3, Poisson’s ratio ν = 0.3, conductivity σ = 2.3 × 106 (Ω·m)−1, shell length L = 0.485 m, shell middle surface radius R = 0.1665 m, shell thickness h = 0.007 m, stiffener height hs = 0.003 m, stiffener spacing bs = 0.097 m, stiffener thickness ds = 0.015 m, and electromagnetic intensity B = 15 T. Considering the transverse excitation amplitude F1 = F2 = 1 × 104 N and in-plane and transverse excitation frequencies Q = Ω = 150 Hz, the periodic and chaotic motions of the stiffened cylindrical shells are analyzed by adopting the fourth-order Runge–Kutta approach. Note that the in-plane excitation amplitude N1 is assumed as a variable, while the basic parameters remained as the fixed values in numerical analysis.To globally view the nonlinear vibration behaviors of the stiffened cylindrical shell, the bifurcation diagrams of the first two modes are depicted in Figure8. Herein, the horizontal axis represents the in-plane excitation amplitude N1 and the vertical axis represents the value of the first two modal displacements W1 and W2. Figure 8 shows that, with the increase of in-plane excitation amplitude, there exist periodic and chaotic motions in two main different regions. As N1 increased past N1 = 5.92 × 107, the system undergoes changes from periodic motions to chaotic motions. The variation results are displayed in Figures 9 and 10. In these figures, (a) and (c) represent the phase portraits on the planesW1,W˙1 and W2,W˙2, respectively; (b) and (d) depict the waveforms on the planes t,W1andt,W2, respectively; (e) and (f) represent three-dimensional phase portrait in the space of W1,W˙1,W2 and the Poincare map on the plane W1,W˙1, respectively; Figure 9(f) shows one point in Poincare maps, which indicates that the system makes the periodic motions when N1 = 4 × 107 N/m2. Under the circumstances of N1 = 7 × 107 N/m2, we observe that each excitation cycle produces different motion states in Figures 10(a)–10(e) and the Poincare map shows a lot of irregular points in Figure 10(f); the dynamic phenomenon illustrates that the system makes the chaotic motion behavior. Note that N1 indicates the in-plane excitation amplitude, while W1 and W2 represent the first and second modal deformation, respectively.Figure 8
Bifurcation diagrams of stiffened cylindrical shell. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 9
Periodic motion of the stiffened cylindrical shell whenN1 = 4 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 10
Chaotic motion of the stiffened cylindrical shell whenN1 = 7 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)


### 5.2. Case 2

The flexural, compressive, or bearing capacity of stiffened cylindrical shells is seriously affected by the geometric configurations of stiffeners. To improve the mechanical performance (such as strength and stiffness) and increase the structural reliability of stiffened cylindrical shells, the effects of stiffener number on nonlinear responses under electromagnetic environment are investigated in this subsection. The geometric parameters of the stiffened cylindrical shell are set as follows: shell lengthL = 0.485 m, shell middle surface radius R = 0.1665 m, shell thickness h = 0.007 m, stiffener number n = 10, stiffener height hs = 0.003 m, stiffener spacing bs = 0.0485 m, and stiffener thickness ds = 0.015 m.The bifurcation diagrams of the first two modes are described in Figure11. From Figure 11, it can be seen that the shell’s motions change from periodic motions to complex chaotic motions when N1 = 8.155 × 107 N/m2. Figure 12 describes the nonlinear dynamic responses when in-plane excitation amplitude N1 = 6.5 × 107 N/m2. At this time, there is only one fixed point in the Poincare map, which indicates that the stiffened cylindrical shell is in single period motion. When N1 = 9.5 × 107 N/m2, the nonlinear dynamic responses of the stiffened cylindrical shell are drawn in Figure 13: there are infinite irregular points in the Poincare map and the motion in the time history diagram is also irregular; this phenomenon indicates that the stiffened cylindrical shell makes chaotic motion behavior. Compared with Figures 8 and 11, we observed that when increasing stiffener number n from 5 to 10, the critical value of N1 will increase from N1 = 5.92 × 107 N/m2 to 8.155 × 107 N/m2, which illustrates that the chaotic motion of stiffened cylindrical shells can be effectively restrained by increasing the stiffener number.Figure 11
Bifurcation diagrams of stiffened cylindrical shell with stiffener numbern = 10. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 12
Periodic motion of the stiffened cylindrical shell whenN1 = 6.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 13
Chaotic motion of the stiffened cylindrical shell whenN1 = 9.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)


### 5.3. Case 3

To verify the effects of various stiffening configurations on the nonlinear dynamic responses, the periodic and chaotic motions of the modes for stiffened cylindrical shells are analyzed. The geometries of the stiffened cylindrical shell are chosen as follows: shell lengthL = 0.485 m, shell middle surface radius R = 0.1665 m, shell thickness h = 0.007 m, stiffener height hs = 0.006 m, stiffener spacing bs = 0.097 m, stiffener thickness ds = 0.015 m, and electromagnetic intensity B = 15 T. The bifurcation diagrams when stiffener height hs = 0.006 m are drawn in Figure 14. From Figure 14, it can be seen that the shell’s motions change from periodicities to complex chaotic motions with the increase of in-plane excitation amplitude N1. And the chaotic motion of the system would emerge when N1 = 9.45 × 107 N/m2. Moreover, the nonlinear dynamic responses of N1 = 6.5 × 107 N/m2 and 9.5 × 107 N/m2 are shown in Figures 15 and 16, respectively.Figure 14
Bifurcation diagrams of stiffened cylindrical shell with stiffener heighths = 0.006 m. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 15
Periodic motion of the stiffened cylindrical shell whenN1 = 6.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 16
Chaotic motion of the stiffened cylindrical shell whenN1 = 9.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)From the comparison of Figures8 and 14, we find that the critical value of N1 increases from 5.92 × 107 N/m2 to 9.45 × 107 N/m2 as the stiffener height increases from 0.003 m to 0.006 m. It can be seen that the higher stiffener height can effectively suppress the chaotic bifurcation behavior of the stiffened cylindrical shell.


### 5.4. Case 4

Considering that different aspect ratios of stiffened cylindrical shells in different operating environments will show vastly different response characteristics, we further explore the influence of aspect ratio on dynamic responses of stiffened cylindrical shells in this subsection. By setting the aspect ratio asL/R = 3.5, the bifurcation diagrams of the first two modes are described in Figure 17. As shown in Figure 17, the shell’s motions translate periodic motions to chaotic motions when the critical value N1 reaches 8.91 × 107 N/m2. Moreover, the nonlinear dynamic responses of N1 = 7 × 107 N/m2 and 10.19 × 107 N/m2 are depicted in Figures 18 and 19, respectively. By comparing Figures 8 and 17, we find that, with the increase of aspect ratio from 2 to 3.5, the critical value of in-plane excitation N1 increases from 5.92 × 107 N/m2 to 9.37 × 107 N/m2, which indicates that the larger aspect ratio can well prevent the chaotic phenomenon in the system.Figure 17
Bifurcation diagrams of stiffened cylindrical shell with length diameter ratioL/R = 3.5. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 18
Periodic motion of the stiffened cylindrical shell whenN1 = 7 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 19
Chaotic motion of the stiffened cylindrical shell whenN1 = 10.19 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Note that, in order to clearly explore the nonlinear behavior under various parameters, the representative parameters (i.e., stiffener number, stiffener height, and aspect ratio) are selected. Although the other stiffened cylindrical shells with different parameter combinations can also be calculated by adopting the proposed method, it should be explored in further investigations.


## 5.1. Case 1

To describe the nonlinear vibration responses of stiffened cylindrical shells, the basic parameters of the stiffened cylindrical shell are set as follows: elastic modulusE = 2.1 × 1011 Pa, density ρ = 7800 kg/m3, Poisson’s ratio ν = 0.3, conductivity σ = 2.3 × 106 (Ω·m)−1, shell length L = 0.485 m, shell middle surface radius R = 0.1665 m, shell thickness h = 0.007 m, stiffener height hs = 0.003 m, stiffener spacing bs = 0.097 m, stiffener thickness ds = 0.015 m, and electromagnetic intensity B = 15 T. Considering the transverse excitation amplitude F1 = F2 = 1 × 104 N and in-plane and transverse excitation frequencies Q = Ω = 150 Hz, the periodic and chaotic motions of the stiffened cylindrical shells are analyzed by adopting the fourth-order Runge–Kutta approach. Note that the in-plane excitation amplitude N1 is assumed as a variable, while the basic parameters remained as the fixed values in numerical analysis.To globally view the nonlinear vibration behaviors of the stiffened cylindrical shell, the bifurcation diagrams of the first two modes are depicted in Figure8. Herein, the horizontal axis represents the in-plane excitation amplitude N1 and the vertical axis represents the value of the first two modal displacements W1 and W2. Figure 8 shows that, with the increase of in-plane excitation amplitude, there exist periodic and chaotic motions in two main different regions. As N1 increased past N1 = 5.92 × 107, the system undergoes changes from periodic motions to chaotic motions. The variation results are displayed in Figures 9 and 10. In these figures, (a) and (c) represent the phase portraits on the planesW1,W˙1 and W2,W˙2, respectively; (b) and (d) depict the waveforms on the planes t,W1andt,W2, respectively; (e) and (f) represent three-dimensional phase portrait in the space of W1,W˙1,W2 and the Poincare map on the plane W1,W˙1, respectively; Figure 9(f) shows one point in Poincare maps, which indicates that the system makes the periodic motions when N1 = 4 × 107 N/m2. Under the circumstances of N1 = 7 × 107 N/m2, we observe that each excitation cycle produces different motion states in Figures 10(a)–10(e) and the Poincare map shows a lot of irregular points in Figure 10(f); the dynamic phenomenon illustrates that the system makes the chaotic motion behavior. Note that N1 indicates the in-plane excitation amplitude, while W1 and W2 represent the first and second modal deformation, respectively.Figure 8
Bifurcation diagrams of stiffened cylindrical shell. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 9
Periodic motion of the stiffened cylindrical shell whenN1 = 4 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 10
Chaotic motion of the stiffened cylindrical shell whenN1 = 7 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)


## 5.2. Case 2

The flexural, compressive, or bearing capacity of stiffened cylindrical shells is seriously affected by the geometric configurations of stiffeners. To improve the mechanical performance (such as strength and stiffness) and increase the structural reliability of stiffened cylindrical shells, the effects of stiffener number on nonlinear responses under electromagnetic environment are investigated in this subsection. The geometric parameters of the stiffened cylindrical shell are set as follows: shell lengthL = 0.485 m, shell middle surface radius R = 0.1665 m, shell thickness h = 0.007 m, stiffener number n = 10, stiffener height hs = 0.003 m, stiffener spacing bs = 0.0485 m, and stiffener thickness ds = 0.015 m.The bifurcation diagrams of the first two modes are described in Figure11. From Figure 11, it can be seen that the shell’s motions change from periodic motions to complex chaotic motions when N1 = 8.155 × 107 N/m2. Figure 12 describes the nonlinear dynamic responses when in-plane excitation amplitude N1 = 6.5 × 107 N/m2. At this time, there is only one fixed point in the Poincare map, which indicates that the stiffened cylindrical shell is in single period motion. When N1 = 9.5 × 107 N/m2, the nonlinear dynamic responses of the stiffened cylindrical shell are drawn in Figure 13: there are infinite irregular points in the Poincare map and the motion in the time history diagram is also irregular; this phenomenon indicates that the stiffened cylindrical shell makes chaotic motion behavior. Compared with Figures 8 and 11, we observed that when increasing stiffener number n from 5 to 10, the critical value of N1 will increase from N1 = 5.92 × 107 N/m2 to 8.155 × 107 N/m2, which illustrates that the chaotic motion of stiffened cylindrical shells can be effectively restrained by increasing the stiffener number.Figure 11
Bifurcation diagrams of stiffened cylindrical shell with stiffener numbern = 10. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 12
Periodic motion of the stiffened cylindrical shell whenN1 = 6.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 13
Chaotic motion of the stiffened cylindrical shell whenN1 = 9.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)


## 5.3. Case 3

To verify the effects of various stiffening configurations on the nonlinear dynamic responses, the periodic and chaotic motions of the modes for stiffened cylindrical shells are analyzed. The geometries of the stiffened cylindrical shell are chosen as follows: shell lengthL = 0.485 m, shell middle surface radius R = 0.1665 m, shell thickness h = 0.007 m, stiffener height hs = 0.006 m, stiffener spacing bs = 0.097 m, stiffener thickness ds = 0.015 m, and electromagnetic intensity B = 15 T. The bifurcation diagrams when stiffener height hs = 0.006 m are drawn in Figure 14. From Figure 14, it can be seen that the shell’s motions change from periodicities to complex chaotic motions with the increase of in-plane excitation amplitude N1. And the chaotic motion of the system would emerge when N1 = 9.45 × 107 N/m2. Moreover, the nonlinear dynamic responses of N1 = 6.5 × 107 N/m2 and 9.5 × 107 N/m2 are shown in Figures 15 and 16, respectively.Figure 14
Bifurcation diagrams of stiffened cylindrical shell with stiffener heighths = 0.006 m. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 15
Periodic motion of the stiffened cylindrical shell whenN1 = 6.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 16
Chaotic motion of the stiffened cylindrical shell whenN1 = 9.5 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)From the comparison of Figures8 and 14, we find that the critical value of N1 increases from 5.92 × 107 N/m2 to 9.45 × 107 N/m2 as the stiffener height increases from 0.003 m to 0.006 m. It can be seen that the higher stiffener height can effectively suppress the chaotic bifurcation behavior of the stiffened cylindrical shell.


## 5.4. Case 4

Considering that different aspect ratios of stiffened cylindrical shells in different operating environments will show vastly different response characteristics, we further explore the influence of aspect ratio on dynamic responses of stiffened cylindrical shells in this subsection. By setting the aspect ratio asL/R = 3.5, the bifurcation diagrams of the first two modes are described in Figure 17. As shown in Figure 17, the shell’s motions translate periodic motions to chaotic motions when the critical value N1 reaches 8.91 × 107 N/m2. Moreover, the nonlinear dynamic responses of N1 = 7 × 107 N/m2 and 10.19 × 107 N/m2 are depicted in Figures 18 and 19, respectively. By comparing Figures 8 and 17, we find that, with the increase of aspect ratio from 2 to 3.5, the critical value of in-plane excitation N1 increases from 5.92 × 107 N/m2 to 9.37 × 107 N/m2, which indicates that the larger aspect ratio can well prevent the chaotic phenomenon in the system.Figure 17
Bifurcation diagrams of stiffened cylindrical shell with length diameter ratioL/R = 3.5. (a) First-order bifurcation diagram. (b) Second-order bifurcation diagram.
(a)(b)Figure 18
Periodic motion of the stiffened cylindrical shell whenN1 = 7 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Figure 19
Chaotic motion of the stiffened cylindrical shell whenN1 = 10.19 × 107 N/m2. (a) Phase portraits on the plane W1,W˙1. (b) Waveform on the plane t,W1. (c) Phase portraits on the plane W2,W˙2. (d) Waveform on the plane t,W2. (e) Three-dimensional phase portrait W1,W˙1,W2. (f) Poincare map on plane W1,W˙1.
(a)(b)(c)(d)(e)(f)Note that, in order to clearly explore the nonlinear behavior under various parameters, the representative parameters (i.e., stiffener number, stiffener height, and aspect ratio) are selected. Although the other stiffened cylindrical shells with different parameter combinations can also be calculated by adopting the proposed method, it should be explored in further investigations.


## 6. Conclusions

In this study, the nonlinear vibration analysis for stiffened cylindrical shells subjected to electromagnetic excitations, transverse excitations, and in-plane excitations is studied for the first time. Firstly, in the framework of the first-order shear deformation theory, the nonlinear partial differential governing equations of motion for stiffened cylindrical shells are derived by using Hamilton’s principle. Then, the nonlinear partial differential governing equations of motion are diverted into a set of coupled nonlinear ordinary differential equations by the Galerkin discretization procedure. By using the multiscale method of perturbation analysis, we obtain four-dimensional nonlinear averaged equations in the polar forms under the case of 1 : 2 internal resonance and principal resonance-1/2 subharmonic parametric resonance. At last, we investigate the effects of different parameters on the resonant responses and the chaotic dynamics of the stiffened cylindrical shell. During the aforementioned analysis, the nonlinear vibrations of stiffened cylindrical shells can be concluded as follows:(1)
From the analysis of resonant responses, we find that multiple amplitudes and jump phenomena have appeared in amplitude-frequency curves. By decreasing the external excitation or increasing the magnetic intensity, the unstable regions of stiffened cylindrical shells can be narrowed and the test accuracy of armature structures can be guaranteed effectively.(2)
With the numerical simulation of periodic and chaotic motions, the bifurcation diagrams, phase portraits, waveforms, and Poincare maps are acquired. From the simulation results, we discover that the periodic motions are translated to chaotic motions as the increase of in-plane excitation. Therefore, the working frequency range of armature structures can be expanded by reducing the in-plane excitation.(3)
With the influence investigations of geometric configurations on nonlinear dynamics, we notice that, with the increase of stiffener number, stiffener size, and aspect ratio, the chaotic phenomenon can be suppressed effectively and the high-accuracy vibration waveform can be produced successfully. Thus, the efforts are of great significance.Since unstable nonlinear dynamic behaviors impose huge waveform deviations of armature structures, we can control multiple excitations (i.e., external excitation and magnetic intensity) and geometric configurations (i.e., stiffener number, stiffener size, and aspect ratio) in order to enhance the test accuracy of armature structures. The presented approach in this paper can provide an eﬃcient analytical framework for nonlinear dynamics theories of stiffened cylindrical shells and will supply important technique guidance for armature structure design in vibration test engineering. Note that the analysis results of this paper are calculated by strict mathematical reasoning, which can effectively ensure its correctness. Other methods will be used to solve partial differential equations of dynamics equation (16), to further verify the proposed numerical model and its solutions.



---

*Source: 9983459-2021-07-20.xml*