# Parametric and Internal Resonances of an Axially Moving Beam with Time-Dependent Velocity

Authors: Bamadev Sahoo; L. N. Panda; G. Pohit
Journal: Modelling and Simulation in Engineering
 (2013)

---

## Abstract

The nonlinear vibration of a travelling beam subjected to principal parametric resonance in presence of internal resonance is investigated. The beam velocity is assumed to be comprised of a constant mean value along with a harmonically varying component. The stretching of neutral axis introduces geometric cubic nonlinearity in the equation of motion of the beam. The natural frequency of second mode is approximately three times that of first mode; a three-to-one internal resonance is possible. The method of multiple scales (MMS) is directly applied to the governing nonlinear equations and the associated boundary conditions. The nonlinear steady state response along with the stability and bifurcation of the beam is investigated. The system exhibits pitchfork, Hopf, and saddle node bifurcations under different control parameters. The dynamic solutions in the periodic, quasiperiodic, and chaotic forms are captured with the help of time history, phase portraits, and Poincare maps showing the influence of internal resonance.


---

## Body

## 1. Introduction

Band saws, fibre textiles, magnetic tapes, paper sheets, aerial tramways, pipes transporting fluids, thread lines, and belts are some technological examples classified as axially moving continua. Analytical models for axially moving systems have been extensively used in the last few decades. The vast literature on axially moving continua vibration has been reviewed by Wickert and Mote Jr. [1] up to 1988. While a linear analysis provides natural frequencies, mode shapes, and critical speeds, its validity regarding the response of the system diminishes as the vibration amplitude becomes sufficiently large or as the critical speed is approached [2]. In these cases one must resort to a nonlinear analysis. Wickert and Mote Jr. [3, 4] studied the transverse vibration of axially moving strings and beams using an eigenfunction method. They also studied the dynamic response of an axially moving string loaded suspended mass. Wickert [5] presented a detailed study of the nonlinear vibrations and bifurcations of moving beams using the Krylov-Bogoliubov-Mitropolsky asymptotic method. Chakraborty et al. [6, 7] investigated both free and forced vibration of the nonlinear traveling beam using complex normal modes.There are papers devoted to the analysis of the dynamic behavior of traveling systems with time-dependent axial velocity or with time-dependent axial tension force. Öz and Pakdemirli [8] investigated principal parametric resonances and combination resonances of sum and difference types for any two modes for an axially accelerating beam using the method of multiple scales. They found that for combination resonances, instabilities occurred only for additive type but not for difference type. Öz et al. [9] extended the work to nonlinear transverse vibration and stability analysis. Comprehensive review of nonlinear modal interactions is there in [10–12]. Using method of multiple scales Riedel and Tan [13] studied the coupled and forced behavior of an axially moving strip with internal resonance. Özkaya et al. [14] investigated nonlinear transverse vibrations and 3 : 1 internal resonances of a beam with multiple supports and plotted frequency response curves for different support numbers. Bagdatli et al. [15] extended this work to find existence of internal resonance cases between different modes. Chin and Nayfeh [16] investigated three-to-one internal resonances in parametrically excited hinged-clamped beams. They studied principal parametric resonance of the first or the second mode and a combination parametric resonance of additive type of these modes. Panda and Kar [17, 18] investigated the principal parametric resonance of first mode, second mode and combination parametric resonance in presence of 3 : 1 internal resonance of a pipe conveying pulsating fluid with pinned-pinned end conditions. They studied the equilibria of these modulated equations and determined their stability, and then they did extensive dynamic study at typical initial conditions. Sze et al. [19] used incremental harmonic balance method for nonlinear vibration of axially moving beams. They investigated the fundamental, superharmonic, and subharmonic resonance in presence of internal resonance. Huang et al. [20] used the same method to analyse the stability and bifurcation analysis of a traveling beam tuned to 3 : 1 internal resonance with attention to fundamental and subharmonic resonances. Chen et al. [21] investigated dynamic stability of an axially accelerating viscoelastic beam undergoing parametric resonance. Ding et al. [22, 23] used Galerkin methods for finding natural frequencies of high-speed axially moving beams with hybrid boundary conditions. Marynowski and Kapitaniak [24–26] introduced several internal dissipation mechanisms like Kelvin-voigt and Zener in the modeling traveling continua. Pakdemirli and Öz [27] studied the transverse vibration of simply supported axially moving Euler-Bernoulli beam for infinite mode analysis and truncation to resonant modes. Ponomareva and van Horssen [28] investigated transversal vibrations of axially travelling continua based on a string model at the low frequencies and a tensioned beam model at the higher frequencies. Recently, a systematic research on travelling beam was pursued by Ghayesh et al. [29–33] involving nonlinear dynamic phenomenon of a variety of system models. The forced dynamics of an axially moving viscoelastic beam was investigated in [29]. The nonlinear vibrations and stability of an axially moving beam subjected to a distributed harmonic excitation load were investigated in [30]. In addition, the response of the system tuned to a three-to-one internal resonance was also examined. Subcritical dynamics of the system was examined [31] via the pseudoarclength continuation technique, while the global dynamics was investigated using direct time integration.The present work aims to investigate the problem of a simply supported traveling beam under parametric excitation due to the harmonic pulsation in the travelling velocity. For a certain range of mean velocity, the natural frequency of the second mode is approximately three times that of the first mode. This relationship between natural frequencies of first two modes of a system having cubic nonlinearity results in a condition of three-to-one autoparametric or internal resonance leading to energy exchange between the two involved modes through nonlinear modal interaction. In the present work, the principal parametric resonance of the first mode is considered. The system behavior shows pitchfork, Hopf, and saddle node bifurcations in steady state analysis. Decreasing internal frequency detuning parameter affects the amplitude of directly excited first mode and the number of Hopf bifurcation points and shifts the occurrence of jump phenomena. Due to the influence of internal resonance, the dynamic behavior of the system exhibits periodic, quasiperiodic, and chaotic responses. The results are presented in the form of time histories, phase plane portraits, fast Fourier transforms (FFTs), and Poincare maps.


## 2. Formulation of the Problem

For the present work, a uniform horizontal beam simply supported at both ends and travelling with a harmonically variable velocity (Figure1) is considered. The assumptions taken here are (1) the motion of the beam is planar, (2) the uniform cross sections remain plane during the motion and the beam behaves like an Euler-Bernoulli beam in transverse vibration, and (3) the type of nonlinearity is geometric due to the midplane stretching effect of the beam. The nondimensional equation of transverse motion of the beam including the nonlinearity due to midplane stretching [5, 9] and viscous damping [34], along with viscoelastic damping [35], is given by
(1)w¨+2vw˙′+v˙w′+(v2-1)w′′+vf2w′′′′+2εαw˙′′′′+2εμw˙=12vl2w′′∫01w′2dx.
The nondimensional scheme used here is
(2)x=x*L,t=t*PρAL2,w=w*L,v=v*P/ρA,2εα=E*L2(ImE)1/2,2εμ=CL2mEI,vl=EAP,vf=EIPL2,
where the variables with asterisk denote dimensional ones. In the EOM, dot denotes derivatives with respect to time (t) and the prime denotes derivatives with respect to spatial derivative (x). m is mass per unit length, ρ is density, A is cross-sectional area, L is length, vf is nondimensional flexural stiffness, and vl is nondimensional longitudinal stiffness of beam. E* is the coefficient of internal dissipation of the beam material, which is assumed to be Kelvin-Voigt type viscoelastic, and C is the external damping factor. α is nondimensional material damping, and μ is nondimensional viscous damping. Reordering the transverse displacement with the relation, w=εw#, where ε<1, and putting it in the equation of motion (1), the system is converted into a weakly nonlinear one [9].Figure 1
Schematic diagram of an axially traveling simply supported beam with variable velocity.For convenience, the superscript “#” is removed, and the weakly nonlinear equation of motion becomes(3)w¨+2vw˙′+v˙w′+(v2-1)w′′+vf2w′′′′+2εαw˙′′′′+2εμw˙=12εvl2w′′∫01w′2dx,
The variable velocity of the beam is
(4)v=v0+εv1sinΩt,
where v0 is mean velocity, εv1 is the amplitude, and Ω is the frequency of the harmonically varying component. This harmonic component of velocity, which is the parametric excitation term in the equation of motion, gives various parametric responses of the system. Putting (4) in (3), we get the equation of transverse motion
(5)2εαw˙′′′′+vf2w′′′′+[v02+2εv0v1sinΩt-1]w′′+εv1Ωw′cosΩt+2(v0+εv1sinΩt)w˙′+2εμw˙+w¨=12εvl2w′′∫01w′2dx
with boundary conditions
(6)w(0,t)=w(1,t)=w′′(0,t)=w′′(1,t)=0,
The linear natural frequencies for various modes vary with variation in the mean velocity of beam. For a considerable range of mean velocity, the first two natural frequencies become commensurable, leading to modal interactions through internal resonance. The combination of internal resonance with different kinds of parametric resonances gives the system response in the form of directly excited mode and indirectly excited one because of the energy exchange between them [16–18].


## 3. Method of Analysis

An approximate solution to this weakly nonlinear distributed parameter system in the form of a first order uniform expansion by using the direct perturbation technique of method of multiple scales (MMS) [14–18] is aimed. The time scale used here is Tn=εnt, n=0,1,2,3…, and the time derivatives are
(7)ddt=D0+εD1+⋯,d2dt2=D02+2εD0D1+⋯,Dn=∂∂Tn,n=0,1,2,3….
Assuming an expansion of the form
(8)w(x,t,ε)=w0(x,T0,T1)+εw1(x,T0,T1)+⋯.
Substituting (7) and (8) into (5) and (6) and equating coefficients of like powers of ε on both sides, we get
(9)O(ε0):D02w0+2v0D0w0′+(v02-1)w0′′+vf2w0′′′′=0,w0(0,t)=w0(1,t)=w0′′(0,t)=w0′′(1,t)=0,(10)O(ε1):D02w1+2v0D0w1′+vf2w1′′′′+(v02-1)w1′′=-2v0D1w0′-2v1sinΩtD0w0′-2D0D1w0-2αD0w0′′′′-2μD0w0-2v0v1sinΩtw0′′-v1ΩcosΩtw0′+12vl2w0′′∫01w0′2dx=0,w1(0,t)=w1(1,t)=w1′′(0,t)=w1′′(1,t)=0.
The solution of (9) may be written as
(11)w0(T0,T1,x)=∑n=1∞ϕn(x)An(T1)eiωnT0+cc,
where ϕn is the mode shapes, ωn is the natural frequencies, and cc is complex conjugate. The mode shapes are calculated previously [8] as
(12)ϕn(x)=C1n{eiβ1nx-(β4n2-β1n2)(eiβ3n-eiβ1n)(β4n2-β2n2)(eiβ3n-eiβ2n)eiβ2nx-(β4n2-β1n2)(eiβ2n-eiβ1n)(β4n2-β3n2)(eiβ2n-eiβ3n)eiβ3nx+[-1+(β4n2-β1n2)(eiβ3n-eiβ1n)(β4n2-β2n2)(eiβ3n-eiβ2n)+(β4n2-β1n2)(eiβ2n-eiβ1n)(β4n2-β3n2)(eiβ2n-eiβ3n)eiβ4nx]},
where βin are eigenvalues which satisfy the dispersive relation (13) and support condition (14) [8]:
(13)vf2βin4-(v02-1)βin2-2v0ωnβin-ωn2=0,i=1,2,3,4,(14)(ei(β1n+β2n)+ei(β3n+β4n))(β1n2-β2n2)(β3n2-β4n2)+(ei(β1n+β3n)+ei(β2n+β4n))(β2n2-β4n2)×(β3n2-β1n2)+(ei(β2n+β3n)+ei(β1n+β4n))×(β1n2-β4n2)(β2n2-β3n2)=0.
For a fixed velocity, the linear natural frequencies of the travelling beam vary with flexural stiffness (vf). For specific values of flexural stiffness, the lower natural frequencies can be commensurable for a range of mean travelling velocity of the beam. These phenomena result in internal resonance in the system and nonlinear interaction between the involved lower modes. In the present investigation, a three-to-one internal resonance (ω2≈3ω1) is considered for a range of mean velocity of the beam. Also it is assumed that there is no other commensurable frequency relationship with higher modes. The case of principal parametric resonance of the first mode (Ω≈2ω1) for subcritical flow velocities in presence of 3 : 1 internal resonance is analyzed in the present investigation. These first two modes are not in internal resonance with any higher modes, so the higher modes except the first two will decay with time due to the presence of damping and Coriolis terms present in the equation. Hence, the first two modes will contribute to the long term system response [10, 11]. Consequently, we replace (11) with
(15)w0(T0,T1,x)=A1(T1)ϕ1(x)eiω1T0+A2(T1)ϕ2(x)eiω2T0+cc.
Now, we write the frequency relations for the internal resonance and principal parametric resonance as
(16)ω2=3ω1+εσ1,Ω=2ω1+εσ2,
where σ1 and σ2 are detuning parameters. It is worthy to note that Ω=ω2-ω1+ε(σ2-σ1), a combination parametric resonance of the difference type, is also activated simultaneously. Substituting (15) and (16) into (10), we get
(17)D02w1+2v0D0w1′+(v02-1)w1′′+vf2w1′′′′=Γ1eiω1T0+Γ2ei(ω1T0+σ1T1)+Γ3eiω1T0+Γ4ei(ω1T0+σ1T1-σ2T1)+Γ5eiω2T0+Γ6ei(ω2T0-σ1T1)+Γ7ei(ω2T0+σ2T1-σ1T1)+cc+NST,
where the terms Γn are defined in the Appendix section. NST stands for terms that do not produce secular or small divisor terms. As the homogeneous part of (17) with its associated boundary conditions has a nontrivial solution, the corresponding nonhomogeneous problem has a solution only if a solvability condition is satisfied [36, 37]. This requires the right-hand side of (17) to be orthogonal to every solution of the adjoint homogeneous problem, which leads to the complex variable modulation equations for amplitude and phase:
(18)2A1′+8S1A12A-1+8S2A1A2A-2+8g1A-12A2eiσ1T1+2μC1A1+2αe1A1+2K1A-1eiσ2T1+2K2A2ei(σ1-σ2)T1=0,2A2′+8S4A22A-2+8S3A1A2A-1+8g2A13e-iσ1T1+2μC2A2+2αe2A2+2K3A1ei(σ2-σ1)T1=0,
where the prime denotes the differentiation with respect to slow time T1 and Si, gi, Ki, Ci, and ei are defined in the Appendix section. Overbar indicates complex conjugate. The terms in the previous equations involving the internal frequency detuning parameter σ1 are the contributions of the internal resonance in the system.


## 4. Stability and Bifurcations

The evolutions of the equilibrium solutions and their stability and bifurcation analysis for principal parametric resonance of first mode are carried out from the modulation equation (18). The Cartesian transformation is used for the complex amplitude as
(19)An=12[pn(T1)-iqn(T1)]eiλn(T1),n=1,2.
Putting this in (18), simplifying by trigonometric manipulations, and separating the real and imaginary parts, we get the normalized reduced equations or the Cartesian form of modulation equations:
(20)p1′=-ϑ1q1-S1R(p13+p1q12)-S1I(p12q1+q13)-S2R(p1p22+p1q22)-S2I(q1p22+q1q22)-g1R(p12p2-p2q12+2p1q1q2)+g1I(2p1q1p2-p12q2+q12q2)-μC1Rp1-μC1Iq1-αe1Rp1-αe1Iq1-K1Rp1+K1Iq1-K2Rp2-K2Iq2,q1′=ϑ1p1+S1I(p13+p1q12)-S1R(p12q1+q13)-S2R(q1p22+q1q22)+S2I(p1p22+p1q22)+g1R(2p1q1p2-p12q2+q12q2)+g1I(2p1q1q2+p12p2-p2q12)+K1Rq1+K1Ip1-K2Rq2+K2Ip2-μC1Rq1+μC1Ip1-αe1Rq1+αe1Ip1,p2′=-ϑ2q2-S4R(p23+p2q22)-S4I(q23+p22q2)-S3R(p12p2+p2q12)-S3I(p12q2+q12q2)-g2R(p13-3p1q12)+g2I(q13-3p12q1)-K3Rp1-K3Iq1-μC2Rp2-μC2Iq2-αe2Rp2-αe2Iq2,q2′=-ϑ2p2-S4R(q23+p22q2)+S4I(p23+p2q22)-S3R(p12q2+q12q2)+S3I(p12p2+p2q12)+g2R(q13-3p12q1)+g2I(p13-3p1q12)-K3Rq1+K3Ip1-μC2Rq2+μC2Ip2-αe2Rq2+αe2Ip2,
where
(21)ϑ1=0.5σ2,ϑ2=1.5σ2-σ1.
The previous equations are perturbed to evaluate the stability. The perturbed equation is
(22){Δp1′Δq1′Δp2′Δq2′}T=[Jc]{Δp1Δq1Δp2Δq2}T,
where T denotes transpose and [Jc] is the Jacobian matrix whose eigenvalues determine the stability and bifurcation of the system. The stability boundary for trivial state is obtained by setting p1=q1=p2=q2=0. The nonlinear steady state response behavior of the system is obtained from the normalized reduced equation (20) by setting p1′=q1′=p2′=q2′=0 and then solving the resulting set of nonlinear algebraic equations. The same set of equations is also used for the analysis of stability and bifurcation of trivial as well as nontrivial solutions. The analysis for dynamic solutions is carried out by numerically integrating (20) with different combinations of system parameters.


## 5. Results and Discussions

The natural frequencies of the beam are numerically evaluated at different mean velocities (v0) with flexural stiffness vf=0.2 by simultaneous solution of dispersive relation (13) and support condition (14). The results are presented in Table 1. It is noticed that, at nondimensional mean velocity v0=0.513, the natural frequency of second mode is approximately equal to three times that of the first mode implying the existence of 3 : 1 internal resonance. It is also noticed that there are no other commensurable frequency relationships involving higher modes. Therefore, nonlinear interaction among higher modes is ruled out. The investigation is limited to the case of principal parametric resonance of first mode, that is, Ω≈2ω1, in presence of internal resonance in the subcritical mean velocity regime of a travelling beam.Table 1
Variation of natural frequencies for the first two modes with mean travelling velocityv0 showing 3 : 1 internal resonance for vf=0.2.
v
0
ω
1
ω
2
3
ω
1
ε
σ
1
=
ω
2
-
3
ω
1
(
ε
σ
1
/
ω
1
)
%
0.4500
3.3149
9.7032
9.9447
−0.2415
−7.2853
0.4600
3.2969
9.6857
9.8907
−0.2050
−6.2180
0.4700
3.2786
9.6677
9.8358
−0.1681
−5.1272
0.4800
3.2599
9.6493
9.7797
−0.1304
−4.0001
0.4900
3.2407
9.6305
9.7221
−0.0916
−2.8265
0.5000
3.2211
9.6113
9.6633
−0.0520
−1.6144
0.5100
3.2011
9.5917
9.6033
−0.0116
−0.3624
0.5110
3.1991
9.5897
9.5973
−0.0076
−0.2376
0.5120
3.1970
9.5878
9.5910
−0.0032
−0.1001
0.5130
3.1950
9.5858
9.5850
0.0008
0.0250
0.5140
3.1930
9.5838
9.5790
0.0048
0.1503
0.5600
3.0948
9.4876
9.2844
0.2032
6.5659
0.5800
3.0493
9.4431
9.1479
0.2952
9.6809
0.6000
3.0020
9.3968
9.0060
0.3908
13.0180The trivial state stability boundary shown in Figure2 is plotted in terms of principal parametric frequency detuning (σ2) and amplitude of fluctuating velocity component (v1) for system parameters vf=0.2, vl=40, v0=0.7, ω1=2.7388, and ω2=9.1403 and for different damping values. The book keeping parameter is taken as ε=0.01, and the corresponding internal frequency detuning parameter is assumed to be σ1=92.39. The region inside the boundary denotes instability. Higher values of damping have the effect of raising and narrowing the instability zones. It is revealed that the effect of material damping (α) raises and narrows the instability zone more compared to that of viscous damping (μ) on the trivial state stability boundary. However, this trivial state stability plot may not reveal completely the real system behavior, because in the unstable zone of the trivial stability plot, the system may have stable nontrivial equilibrium solution or stable dynamic solution like periodic or quasiperiodic solution. In addition, there is a possibility of chaotic solutions or multiple stable solutions. For this reason, it is required to carry out the dynamic analysis as well as the stability and bifurcation study of equilibrium solutions of the system.Figure 2
Trivial state boundary for different damping parameters. Values of the nondimensional damping parameters(μ,α), indicated on the curves.


### 5.1. Stability and Bifurcations of Equilibrium Solutions

Continuation algorithm is used to determine the nonlinear steady state response by solving the set of algebraic equations generated after settingpi′=qi′=0 in the normalized reduced equation (20). The stability and bifurcation of the equilibrium solutions are obtained from the eigenvalues of the Jacobean matrix at each point of the solution. In order to validate the results obtained by the present analysis, the frequency response and amplitude response curves of Chin et al. [16] are generated once again using continuation algorithm. They are shown in Figure 3, and the results are found to be in good agreement. Since the frequency and amplitude response curves are symmetrical about σ2 and v1 axes, respectively, only positive sides of the response curves are shown.Frequency response curves as obtained by continuation algorithm for the first and second modes when the first mode is parametrically excited for the system parameters, Chin and Nayfeh [16].
(a)
(b)Frequency response curves are obtained against variation in frequency detuning parameterσ2 for first and second mode for μ=0.1, α=0, v1=10, vl=40, and σ1=92.39 and are shown in Figure 4. The normal continuous lines in the figure represent stable equilibrium solutions, the bold lines represent unstable foci, and the dotted lines denote saddles. Different parameter values for characteristic points on different branches of the plot are indicated in Table 2. The response curves exhibit a hardening-spring type of nonlinearity. With increase in σ2 from a small value, the trivial stable solution loses stability at σ2=-22.071, through supercritical pitchfork bifurcation, and results in a two-mode nontrivial stable equilibrium solution. It is observed that the amplitude of the first mode increases initially then decreases, but the amplitude of second mode increases monotonically. When the value of σ2 increases, the equilibrium solution becomes unstable at H1(σ2=206.588) through Hopf bifurcation, and out of two pairs of complex conjugate eigenvalues, one pair crosses the imaginary axis from the left half of the complex plane to the right half. With further increase of σ2, the same state continues until a saddle node bifurcation occurs at SN(σ2=301.873), where the system response jumps to one of the two stable equilibrium branches, one trivial and the other nontrivial, depending on the initial conditions as the solution converges to the closer equilibrium state as per the concept of region of attraction. With further increase in frequency detuning parameter, amplitude of the first mode increases monotonically along the nontrivial branch, whereas the amplitude of second mode decreases continuously. Thus the amplitude of the indirectly excited second mode is limited to a fixed higher magnitude, and for high values of σ2, it becomes stagnant at fixed low amplitude while there is no such limitation for the directly excited first mode.Table 2
Typical points on the different branches of the frequency response forμ=0.1, α=0, v1=10, v1=40, and σ1=92.39.
p
1
q
1
p
2
q
2
σ
2
Bifurcation point
0.00500579
−0.01569813
0.00004925
0.01964672
100
−0.02181265
−0.00742674
0.01533298
0.00020989
100
0.010164208
0.00356793
0.01833426
0.00011845
100
−0.004819409
0.01479703
0.00009386
0.01071881
100
−0.007939
−0.003674
0.042869
0.006271
206.588482
H
1
−0.002453
−0. 004177
0. 039646
0. 039917
301.873000
SN
0.032196
0.010862
0.005854
0.000026
190.033742
H
2
0.026628
0.008991
0.010665
0.000036
144.502843
H
3
0.022130
0.007495
0.015112
0.000039
127.980336
H
4
0.004908
0.001822
0.013291
0.000245
80.408029
H
5
0.000000
0.000000
0.000000
0.000000
49.474028
H
6
0.000000
0.000000
0.000000
0.000000
44.608834
H
7Frequency response curves for the (a) first mode and (b) second mode when the first mode is parametrically excited for the system parametersμ=0.1, α=0, v1=10, vl=40, and σ1=92.39.
(a)
(b)When the detuning parameter decreases from a high value, the system follows either trivial or nontrivial stable equilibrium path, depending on the initial conditions. If the solution is nontrivial, with decrease ofσ2 value, the nontrivial stable branch loses stability via Hopf bifurcation at H2(σ2=190.033) and regains stability via a reverse Hopf bifurcation at H3(σ2=144.502). When the detuning parameter is further decreased, again the system loses stability via Hopf bifurcation at H4(σ2=127.98) and regains stability via a reverse Hopf bifurcation at H5(σ2=80.408) on the same path. With further decrease in frequency detuning parameter, the nontrivial stable equilibrium branch merges with stable trivial equilibrium solution, the system losing and regaining the stability at H6(σ2=49.474) and H7(σ2=44.608), respectively. At σ2=14.190, the trivial equilibrium solution loses stability via subcritical/reverse pitchfork bifurcation and results in a jump of the response to the stable nontrivial branch of the solution. Again with further decrease of frequency detuning, the nontrivial stable solution branch loses stability through pitchfork bifurcation at σ2=-22.071, giving the way to trivial solution. The directly excited first mode dominates the indirectly excited second mode.Figure5 shows the frequency response curves for first and second mode of the system for higher amplitude of the fluctuating velocity component (v1). The system parameters considered are μ=0.1, α=0, v1=15, vl=40, and σ1=92.39. Even though the solution curves are similar in shape to the curves obtained in case of lower amplitude of excitation (v1=10) as shown in Figure 4, the jump phenomena at the saddle node bifurcation (SN) occur at higher value of the detuning parameter (σ2=457.7476) compared to the previous case (σ2=301.8730 (Figure 4)). In addition, unstable zone in trivial solution gets broadened which is commensurate with the trivial state stability plot.Frequency response curves for the (a) first mode and (b) second mode when the first mode is parametrically excited for the system parametersμ=0.1, α=0, v1=15, vl=40, and σ1=92.39.
(a)
(b)Figure6 shows typical frequency response curves for two modes considering the effect of the internal damping for system parameters μ=0, α=0.001, v1=10, vl=40, and σ1=92.39. The strength of nonlinear interaction due to internal resonance gets weakened due to internal damping (Figure 6) compared to the case of external damping (Figure 4). The influence of internal detuning parameter (σ1) on the frequency response is shown in Figure 7. It is evident that the decrease in internal detuning parameter (σ1) to 20.320 (Figures 7(a) and 7(b)) and −27.680 (Figures 7(c) and 7(d)) from 92.390 (Figure 4) weakens the strength of nonlinear interaction due to three-to-one internal resonance. The amplitude of the directly excited first mode decreases more pronouncedly than the indirectly excited second mode. Beside this, the number of Hopf bifurcation points on the upper nontrivial curve decreases from four for σ1=92.39 to two for σ1=20.320 and totally vanishes for σ1=-27.680, and also there is decreasing trend in σ2 value at which saddle node bifurcation occurs. Figure 8 shows the effect of decrease of external damping (μ=0.05) on the frequency response of the system. The nature of the shape of the curves is similar to that of Figure 4, but the jump phenomena occur at a higher value of parametric excitation frequency detuning parameter (σ2=382.9953), and the amplitudes of both directly and indirectly excited modes are amplified.Frequency response curves for the (a) first mode and (b) second mode when the first mode is parametrically excited for the system parametersα=0.001, μ=0, v1=10, vl=40, andσ1=92.39.
(a)
(b)Effect of internal detuning parameter on the frequency response of (a) first mode and (b) second mode for system parametersμ=0.1, α=0, v1=10, vl=40, (a,b)σ1=20.320, and (c,d)σ1=-27.680.
(a)
(b)
(c)
(d)Effect of external damping parameter on the frequency response of first mode (a) and second mode (b) forμ=0.05, α=0, v1=10, vl=40, and σ1=92.39.
(a)
(b)


### 5.2. Dynamic Solutions

Frequency response and amplitude response plots reveal different stability and bifurcations of the equilibrium solutions with variation of control parameters. Dynamic analysis of the system which is dependent on initial conditions is studied in the form of periodic, quasiperiodic, and chaotic responses, and some selected results are presented.Figures9(a)–9(d) show typical system response in terms of phase portraits (a, b) and time traces (c, d) at σ2=68.2799 corresponding to the upper nontrivial stable branch of the frequency response curve (Figure 8) for μ=0.05, α=0, v1=10, and σ1=92.39. The response is periodic about the nontrivial equilibrium solution when the time integration is started with the initial values p1=0.002, q1=0.0007, p2=0.0067, and q2=0.0001. Further along the same branch at σ2=75.2799, the response is quasiperiodic in both modes being more prominent in the second mode as shown in the two-dimensional projections of the phase portraits onto the p-q planes in Figures 12(a) and 12(b) and the time traces in Figures 10(c) and 10(d). The response remains quasiperiodic in both modes for higher frequency detuning parameter values, typically at σ2=76.2799 as shown in the closed loop Poincare maps and FFT power spectra in Figures 11(c), 11(d), and 11(b), respectively. With further increase in the value of detuning parameter, typically at σ2=78.2799, we find the closed loops of Poincare map get merged and give a way to chaotic response in second mode in Figure 12(d). However, in first mode the system response is still quasiperiodic as seen from phase portrait (Figure 12(a)) and Poincare map (Figure 12(c)), respectively.Phase portraits (a, b) and time histories (c, d) in the upper nontrivial stable branch of the frequency response curve of Figure8 for σ2=68.2799, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Phase portraits (a, b) and time histories (c, d) forσ2=75.2799, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Phase portrait (a), FFT power spectra (b), and Poincare maps (c, d) forσ2=76.2799, μ=0.05, v1=10, and α=0, σ1=92.39.
(a)
(b)
(c)
(d)Phase portraits (a, b) and Poincare maps (c, d) forσ2=78.2799, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Figures13(a)–13(d) show the typical system behavior at σ2=115.7201, μ=0.05, α=0, v1=10, and σ1=92.39 in terms of phase portraits and time traces in Figures 13(a), 13(b), 13(c), and 13(d), respectively. The response is initially chaotic and jumps to the nearby stable trivial attractor. With further increase in detuning parameter, the system behavior changes drastically as shown in Figures 14(a)–14(d), at σ2=116.4201 illustrating quasiperiodic motion in the first mode and chaotic motion in second mode. The dynamic response in both modes becomes chaotic at detuning parameter σ2=117.0201 as shown in phase portrait, FFT power spectra, and Poincare maps in Figures 15(a), 15(b), 15(c), and 15(d), respectively. The changes in system response from periodic in both modes to mixed mode, that is, quasiperiodic in first and chaotic in second mode to chaotic in both modes, as explained above, happen in the zone of frequency response plot where all three kinds of curves stable, saddle, and unstable are in very close proximity and crossing each other. Existence of multiple branches is possible due to the presence of internal resonance in the system. The nonlinear modal interaction influences simultaneously both the stable and unstable attractors, which finally results in such varied system responses.Phase portraits (a, b) and time histories (c, d) forσ2=115.7201, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Phase portraits (a, b), time trace (c), and Poincare map (d) forσ2=116.4201, μ=0.05, v1=10, and α=0, σ1=92.39.
(a)
(b)
(c)
(d)Phase portrait (a), FFT power spectra (b), and Poincare maps (c, d) forσ2=117.4201, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Similar investigation is carried in the lower nontrivial stable branch of the frequency response plot of Figure8 as well. For a point on the same branch at σ2=280.5201, the system response exhibits one periodic and one quasiperiodic system behavior as shown in Figure 16 in terms of phase portraits, time trace, and FFT power spectra. The same behavior is noticed for another point on the same branch at σ2=310.5201 though the figures are not presented to avoid repetition. Thus, due to the presence of internal resonance, a wide range of dynamic behavior can be observed with variation of control parameters.Phase portraits (a, b), time trace (c), and FFT power spectra (d) in the lower nontrivial stable branch of the frequency response plot of Figure8 for σ2=280.5201, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)


## 5.1. Stability and Bifurcations of Equilibrium Solutions

Continuation algorithm is used to determine the nonlinear steady state response by solving the set of algebraic equations generated after settingpi′=qi′=0 in the normalized reduced equation (20). The stability and bifurcation of the equilibrium solutions are obtained from the eigenvalues of the Jacobean matrix at each point of the solution. In order to validate the results obtained by the present analysis, the frequency response and amplitude response curves of Chin et al. [16] are generated once again using continuation algorithm. They are shown in Figure 3, and the results are found to be in good agreement. Since the frequency and amplitude response curves are symmetrical about σ2 and v1 axes, respectively, only positive sides of the response curves are shown.Frequency response curves as obtained by continuation algorithm for the first and second modes when the first mode is parametrically excited for the system parameters, Chin and Nayfeh [16].
(a)
(b)Frequency response curves are obtained against variation in frequency detuning parameterσ2 for first and second mode for μ=0.1, α=0, v1=10, vl=40, and σ1=92.39 and are shown in Figure 4. The normal continuous lines in the figure represent stable equilibrium solutions, the bold lines represent unstable foci, and the dotted lines denote saddles. Different parameter values for characteristic points on different branches of the plot are indicated in Table 2. The response curves exhibit a hardening-spring type of nonlinearity. With increase in σ2 from a small value, the trivial stable solution loses stability at σ2=-22.071, through supercritical pitchfork bifurcation, and results in a two-mode nontrivial stable equilibrium solution. It is observed that the amplitude of the first mode increases initially then decreases, but the amplitude of second mode increases monotonically. When the value of σ2 increases, the equilibrium solution becomes unstable at H1(σ2=206.588) through Hopf bifurcation, and out of two pairs of complex conjugate eigenvalues, one pair crosses the imaginary axis from the left half of the complex plane to the right half. With further increase of σ2, the same state continues until a saddle node bifurcation occurs at SN(σ2=301.873), where the system response jumps to one of the two stable equilibrium branches, one trivial and the other nontrivial, depending on the initial conditions as the solution converges to the closer equilibrium state as per the concept of region of attraction. With further increase in frequency detuning parameter, amplitude of the first mode increases monotonically along the nontrivial branch, whereas the amplitude of second mode decreases continuously. Thus the amplitude of the indirectly excited second mode is limited to a fixed higher magnitude, and for high values of σ2, it becomes stagnant at fixed low amplitude while there is no such limitation for the directly excited first mode.Table 2
Typical points on the different branches of the frequency response forμ=0.1, α=0, v1=10, v1=40, and σ1=92.39.
p
1
q
1
p
2
q
2
σ
2
Bifurcation point
0.00500579
−0.01569813
0.00004925
0.01964672
100
−0.02181265
−0.00742674
0.01533298
0.00020989
100
0.010164208
0.00356793
0.01833426
0.00011845
100
−0.004819409
0.01479703
0.00009386
0.01071881
100
−0.007939
−0.003674
0.042869
0.006271
206.588482
H
1
−0.002453
−0. 004177
0. 039646
0. 039917
301.873000
SN
0.032196
0.010862
0.005854
0.000026
190.033742
H
2
0.026628
0.008991
0.010665
0.000036
144.502843
H
3
0.022130
0.007495
0.015112
0.000039
127.980336
H
4
0.004908
0.001822
0.013291
0.000245
80.408029
H
5
0.000000
0.000000
0.000000
0.000000
49.474028
H
6
0.000000
0.000000
0.000000
0.000000
44.608834
H
7Frequency response curves for the (a) first mode and (b) second mode when the first mode is parametrically excited for the system parametersμ=0.1, α=0, v1=10, vl=40, and σ1=92.39.
(a)
(b)When the detuning parameter decreases from a high value, the system follows either trivial or nontrivial stable equilibrium path, depending on the initial conditions. If the solution is nontrivial, with decrease ofσ2 value, the nontrivial stable branch loses stability via Hopf bifurcation at H2(σ2=190.033) and regains stability via a reverse Hopf bifurcation at H3(σ2=144.502). When the detuning parameter is further decreased, again the system loses stability via Hopf bifurcation at H4(σ2=127.98) and regains stability via a reverse Hopf bifurcation at H5(σ2=80.408) on the same path. With further decrease in frequency detuning parameter, the nontrivial stable equilibrium branch merges with stable trivial equilibrium solution, the system losing and regaining the stability at H6(σ2=49.474) and H7(σ2=44.608), respectively. At σ2=14.190, the trivial equilibrium solution loses stability via subcritical/reverse pitchfork bifurcation and results in a jump of the response to the stable nontrivial branch of the solution. Again with further decrease of frequency detuning, the nontrivial stable solution branch loses stability through pitchfork bifurcation at σ2=-22.071, giving the way to trivial solution. The directly excited first mode dominates the indirectly excited second mode.Figure5 shows the frequency response curves for first and second mode of the system for higher amplitude of the fluctuating velocity component (v1). The system parameters considered are μ=0.1, α=0, v1=15, vl=40, and σ1=92.39. Even though the solution curves are similar in shape to the curves obtained in case of lower amplitude of excitation (v1=10) as shown in Figure 4, the jump phenomena at the saddle node bifurcation (SN) occur at higher value of the detuning parameter (σ2=457.7476) compared to the previous case (σ2=301.8730 (Figure 4)). In addition, unstable zone in trivial solution gets broadened which is commensurate with the trivial state stability plot.Frequency response curves for the (a) first mode and (b) second mode when the first mode is parametrically excited for the system parametersμ=0.1, α=0, v1=15, vl=40, and σ1=92.39.
(a)
(b)Figure6 shows typical frequency response curves for two modes considering the effect of the internal damping for system parameters μ=0, α=0.001, v1=10, vl=40, and σ1=92.39. The strength of nonlinear interaction due to internal resonance gets weakened due to internal damping (Figure 6) compared to the case of external damping (Figure 4). The influence of internal detuning parameter (σ1) on the frequency response is shown in Figure 7. It is evident that the decrease in internal detuning parameter (σ1) to 20.320 (Figures 7(a) and 7(b)) and −27.680 (Figures 7(c) and 7(d)) from 92.390 (Figure 4) weakens the strength of nonlinear interaction due to three-to-one internal resonance. The amplitude of the directly excited first mode decreases more pronouncedly than the indirectly excited second mode. Beside this, the number of Hopf bifurcation points on the upper nontrivial curve decreases from four for σ1=92.39 to two for σ1=20.320 and totally vanishes for σ1=-27.680, and also there is decreasing trend in σ2 value at which saddle node bifurcation occurs. Figure 8 shows the effect of decrease of external damping (μ=0.05) on the frequency response of the system. The nature of the shape of the curves is similar to that of Figure 4, but the jump phenomena occur at a higher value of parametric excitation frequency detuning parameter (σ2=382.9953), and the amplitudes of both directly and indirectly excited modes are amplified.Frequency response curves for the (a) first mode and (b) second mode when the first mode is parametrically excited for the system parametersα=0.001, μ=0, v1=10, vl=40, andσ1=92.39.
(a)
(b)Effect of internal detuning parameter on the frequency response of (a) first mode and (b) second mode for system parametersμ=0.1, α=0, v1=10, vl=40, (a,b)σ1=20.320, and (c,d)σ1=-27.680.
(a)
(b)
(c)
(d)Effect of external damping parameter on the frequency response of first mode (a) and second mode (b) forμ=0.05, α=0, v1=10, vl=40, and σ1=92.39.
(a)
(b)


## 5.2. Dynamic Solutions

Frequency response and amplitude response plots reveal different stability and bifurcations of the equilibrium solutions with variation of control parameters. Dynamic analysis of the system which is dependent on initial conditions is studied in the form of periodic, quasiperiodic, and chaotic responses, and some selected results are presented.Figures9(a)–9(d) show typical system response in terms of phase portraits (a, b) and time traces (c, d) at σ2=68.2799 corresponding to the upper nontrivial stable branch of the frequency response curve (Figure 8) for μ=0.05, α=0, v1=10, and σ1=92.39. The response is periodic about the nontrivial equilibrium solution when the time integration is started with the initial values p1=0.002, q1=0.0007, p2=0.0067, and q2=0.0001. Further along the same branch at σ2=75.2799, the response is quasiperiodic in both modes being more prominent in the second mode as shown in the two-dimensional projections of the phase portraits onto the p-q planes in Figures 12(a) and 12(b) and the time traces in Figures 10(c) and 10(d). The response remains quasiperiodic in both modes for higher frequency detuning parameter values, typically at σ2=76.2799 as shown in the closed loop Poincare maps and FFT power spectra in Figures 11(c), 11(d), and 11(b), respectively. With further increase in the value of detuning parameter, typically at σ2=78.2799, we find the closed loops of Poincare map get merged and give a way to chaotic response in second mode in Figure 12(d). However, in first mode the system response is still quasiperiodic as seen from phase portrait (Figure 12(a)) and Poincare map (Figure 12(c)), respectively.Phase portraits (a, b) and time histories (c, d) in the upper nontrivial stable branch of the frequency response curve of Figure8 for σ2=68.2799, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Phase portraits (a, b) and time histories (c, d) forσ2=75.2799, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Phase portrait (a), FFT power spectra (b), and Poincare maps (c, d) forσ2=76.2799, μ=0.05, v1=10, and α=0, σ1=92.39.
(a)
(b)
(c)
(d)Phase portraits (a, b) and Poincare maps (c, d) forσ2=78.2799, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Figures13(a)–13(d) show the typical system behavior at σ2=115.7201, μ=0.05, α=0, v1=10, and σ1=92.39 in terms of phase portraits and time traces in Figures 13(a), 13(b), 13(c), and 13(d), respectively. The response is initially chaotic and jumps to the nearby stable trivial attractor. With further increase in detuning parameter, the system behavior changes drastically as shown in Figures 14(a)–14(d), at σ2=116.4201 illustrating quasiperiodic motion in the first mode and chaotic motion in second mode. The dynamic response in both modes becomes chaotic at detuning parameter σ2=117.0201 as shown in phase portrait, FFT power spectra, and Poincare maps in Figures 15(a), 15(b), 15(c), and 15(d), respectively. The changes in system response from periodic in both modes to mixed mode, that is, quasiperiodic in first and chaotic in second mode to chaotic in both modes, as explained above, happen in the zone of frequency response plot where all three kinds of curves stable, saddle, and unstable are in very close proximity and crossing each other. Existence of multiple branches is possible due to the presence of internal resonance in the system. The nonlinear modal interaction influences simultaneously both the stable and unstable attractors, which finally results in such varied system responses.Phase portraits (a, b) and time histories (c, d) forσ2=115.7201, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Phase portraits (a, b), time trace (c), and Poincare map (d) forσ2=116.4201, μ=0.05, v1=10, and α=0, σ1=92.39.
(a)
(b)
(c)
(d)Phase portrait (a), FFT power spectra (b), and Poincare maps (c, d) forσ2=117.4201, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)Similar investigation is carried in the lower nontrivial stable branch of the frequency response plot of Figure8 as well. For a point on the same branch at σ2=280.5201, the system response exhibits one periodic and one quasiperiodic system behavior as shown in Figure 16 in terms of phase portraits, time trace, and FFT power spectra. The same behavior is noticed for another point on the same branch at σ2=310.5201 though the figures are not presented to avoid repetition. Thus, due to the presence of internal resonance, a wide range of dynamic behavior can be observed with variation of control parameters.Phase portraits (a, b), time trace (c), and FFT power spectra (d) in the lower nontrivial stable branch of the frequency response plot of Figure8 for σ2=280.5201, μ=0.05, v1=10, α=0, and σ1=92.39.
(a)
(b)
(c)
(d)


## 6. Conclusions

In the present investigation, principal parametric resonance of first mode in presence of 3 : 1 internal resonance of a beam moving with variable velocity is considered. Stability boundaries of trivial state are obtained for different values of internal and external dissipations. It has been observed that higher values of damping have the effect of raising and narrowing the instability zones. Bifurcations of equilibrium solutions are analyzed in the form response plots. It has been shown in frequency response plot that the nontrivial steady state solutions bifurcate from trivial solutions through supercritical pitchfork bifurcations.Besides the pitchfork bifurcations, the system also experiences Hopf bifurcation and saddle node bifurcation due to variation of different system parameters. Damping decreases the strength of nonlinear interaction due to internal resonance. Increasing amplitude of fluctuating velocity component broadens the range of trivial state instability and increases the value of parametric frequency detuning at which jump phenomena occur. Decreasing internal frequency detuning parameter affects the amplitude of directly excited first mode and number of Hopf bifurcation points. It also shifts the occurrence of jump phenomena.A detailed study is carried out to determine the influence of different control parameters on dynamic behavior of the system. The dynamic solutions in the periodic, quasiperiodic, and chaotic forms, are captured with the help of time history, phase portraits, and Poincare maps. A wide array of dynamic behavior is noticed when nontrivial stable and saddle branches are formed due to internal resonance and also in the zone where the three branches are very close and crossing each other.In case of conventional nontravelling beams with simply supported boundary conditions, occurrence of internal resonance is not possible due to vanishing of the nonlinear interaction coefficients [10]. However, a varied system response is possible in case of travelling beams due to nonlinear modal interaction leading to simultaneous influence of both stable and unstable attractors.



---

*Source: 919517-2013-11-10.xml*