The nonlinear resonant responses, mode interactions, and multitime periodic and chaotic oscillations of the cantilevered pipe conveying pulsating fluid are studied under the harmonic external force in this research. According to the nonlinear dynamic model of the cantilevered beam derived using Hamilton’s principle under the uniformly distributed external harmonic excitation, we combine Galerkin technique and the method of multiple scales together to obtain the average equation of the cantilevered pipe conveying pulsating fluid under 1 : 3 internal resonance and principal parametric resonance. Based on the average equation in the polar form, several amplitude-frequency response curves are obtained corresponding to the certain parameters. It is found that there exist the hardening-spring type behaviors and jumping phenomena in the cantilevered pipe conveying pulsating fluid. The nonlinear oscillations of the cantilevered pipe conveying pulsating fluid can be excited more easily with the increase of the flow velocity, external excitation, and coupling degree of two order modes. Numerical simulations are performed to study the chaos of the cantilevered pipe conveying pulsating fluid with the external harmonic excitation. The simulation results exhibit the existence of the period, multiperiod, and chaotic responses with the variations of the fluid velocity or excitation. It is found that, in the cantilevered pipe conveying pulsating fluid, there are the multitime nonlinear vibrations around the left-mode and the right-mode positions, respectively. We also observe that there exist alternately the periodic and chaotic vibrations of the cantilevered pipe conveying pulsating fluid in the certain range.
National Natural Science Foundation of China116721881183200211427801Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality1. Introduction
Pipes conveying fluid are widely utilized in many engineering fields, such as aeronautic, astronautic, and mechanical engineering systems. It is extremely important for us to ensure the efficient utilization and safe operation of the pipe conveying fluid system, and its stable and safe operations are closely related to all aspects of the personal life and industrial production. However, the applications of the pipes conveying pulsating fluid are particularly challenging because they undergo the large deformations and significant stresses. The large deformations often lead to the nonlinear vibrations of the pipes conveying pulsating fluid. One of the main reasons for the nonlinear vibrations of the pipes conveying pulsating fluid is the time-varying flow speed and external harmonic excitation. Pulsating flow due to the pump operation can cause a parametric excitation loading in the pipes conveying fluid. The nonlinear oscillations of the pipes conveying pulsating fluid will lead to the structure damages. As we all know, there are three typical types of nonlinear oscillations in the structures and systems, namely, the periodic, quasi-periodic, and chaotic oscillations. In fact, the chaotic oscillations of the pipes conveying pulsating fluid are dangerous because the amplitudes of the chaotic oscillations are larger than those of the periodic oscillations, which have been the object of increasing attention in engineering applications. However, there is less research on the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid with 1 : 3 internal resonance when the fluid is transported at a critical speed through the pipe. Therefore, it is of great significance for us to study the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid under the case of 1 : 3 internal resonance.
The pipe conveying fluid mainly consists of three important elements: pipeline, fluid, and external environment. It is necessary to establish a mathematical model for obtaining a reasonable description of the pipes conveying fluid. The beam model is usually used for analysis of the vibration when the pipe diameter is much smaller than the length. The nonlinear dynamic study of the pipe conveying fluid system began in 1980s. Researches for the pipes conveying pulsating fluid have become a hot field of engineering and science [1–4]. Holmes and Marsden [5, 6] established the first dynamic model of motion for the pipe conveying fluid through considering the nonlinear factors. They studied the bifurcation phenomena caused by the velocity and the axial force and summed up the bifurcation motion characteristics of the pipe conveying fluid system. According to Hamilton’s principle and Euler-Bernoulli beam theory, Huo and Wang [7] derived the differential governing equation of motion for a vertical cantilevered pipe conveying fluid when the pipe exhibits the deploying or retracting motions and discussed the influence of the deploying or retracting speed, mass ratio, and fluid velocity on the dynamic responses and stability. Based on the modified strain gradient theory in conjunction with Euler-Bernoulli beam model, Hosseini and Bahaadini [8] studied the size-dependent stability of the cantilever micropipes conveying fluid and examined the influences of the geometric parameters on the natural frequencies and the flutter critical speeds. Askarian et al. [9] researched the dynamic stability of a vertical clamped-free pipe conveying pulsatile flow by using Euler-Bernoulli beam theory.
In addition, some scholars also provided several different mathematical models to investigate the vibrations of the pipes conveying fluid. According to the nonlinear Novozhilov shell theory for the isotropic materials, Tubaldi et al. [10] established a fully coupled fluid-structure interaction model and studied the nonlinear vibrations of the circular cylindrical shells conveying pulsatile flow with the flexible boundary conditions subjected to the pulsatile pressure. Bai et al. [11] simulated the varying density fluid of the vertical cantilevered pipe conveying fluid by a new mathematical model.
After the establishment of the rational mathematical model, the further researches on the vibrations of the pipes conveying fluid are mainly focused on three significant aspects. Firstly, the vibration characteristics and instability conditions of the pipes conveying fluid are studied, which have different boundary conditions, material properties, and functions. Secondly, the problems of the nonlinear dynamics are studied to understand the nonlinear vibration characteristics of the pipes conveying fluid because the pipes can be regarded as the complex nonlinear dynamical systems. Thirdly, the transfer mechanism of the energy between two modes and the internal resonance are studied to avoid the chaotic vibrations of the pipes conveying fluid.
More and more researches about model and analysis have been published with the continuous development of modern computing and analytical techniques. Bajaj et al. [12] researched the nonlinear vibration responses of the pipes conveying fluid and analyzed the influences of flow rate, mass ratio, and pressure on the nonlinear vibrations and found that Hopf bifurcation occurs under the certain conditions. Sri Namchchivaya and Tien [13, 14] examined the nonlinear vibrations of the supported pipes conveying pulsating fluid and found that the trivial solution of the averaged equation loses its stability through the simple or double zero bifurcations in the vicinity of the subharmonic resonance. Based on a semianalytical approach, Sarkar and Paidoussis [15] obtained the proper orthogonal modes to describe the nonlinear oscillations of a cantilevered pipe conveying fluid and to explore the nonlinear dynamics of the pipe by means of the low-dimensional model. McDonald and Sri Namachchivaya [16] studied the local bifurcation behaviors of parametrically excited simply supported pipes conveying fluid and the stability of solution where the energy transfer may happen from the high-frequency to low-frequency vibration modes.
Yoshizawa et al. [17] established the dynamic model of the pipes conveying fluid under the fixed-hinge boundary conditions and studied the vibration responses under the effect of pulse flow. For the pipe conveying fluid system, Hou and Zeng [18] obtained numerical solutions of the transverse vibration equation through using the finite element method. Setoodeh and Afrahim [19] investigated the nonlinear vibrations of the functionally graded materials micropipes conveying fluid. An explicit expression of the nonlinear fundamental frequency was given by using the homotopy analysis method. Kheiri and Païdoussis [20] derived the equation of motion for a typical flexible pipe conveying fluid by using the generalized Hamilton’s principle. Gan et al. [21] simplified the equation of motion to the random variable and researched the vibration characteristics of the pipes conveying fluid clamped at both ends. Zhang et al. [22] established three-dimensional nonlinear equations of motion for the pipes conveying fluid under the general boundary conditions. The natural frequencies of the pipes under different boundary conditions were calculated and the nonlinear dynamic characteristics were analyzed. Wang and Liu [23] studied the transverse vibration and stability of the functionally graded material pipe conveying fluid by utilizing the symplectic method.
Based on the fluid-structure interaction, Liang et al. [24] studied the free vibration of the pipes conveying fluid through using the linear and nonlinear complex mode approach. Liang et al. [25] gave the analysis of the nonlinear free vibration for the spinning viscoelastic pipes conveying fluid. Liang et al. [26] studied the transverse free vibration of the spinning pipes conveying fluid and found that the qualitative stability of the pipes mainly relies on the fluid-structure interaction and mass ratio. Liang et al. [27] established a dynamic model of simply supported spinning pipes conveying fluid with the axial deployment and studied the transverse free vibration involving the time-dependent parameters. Liang et al. [28] researched the coupled flexural-torsional vibrations of the pipes conveying fluid spinning on an eccentric axis. Recently, Liang et al. [29] investigated the parametric vibrations of the pipes conveying fluid by the nonlinear normal modes and numerical iterative approach.
Several papers researched the nonlinear dynamic characteristics of the pipes conveying fluid. Li and Paidoussis [30] studied the nonlinear oscillations of a standing cantilevered pipe conveying fluid and showed that the chaotic oscillations of the cantilevered pipe exist when the gravity parameter was sufficiently perturbed off the doubly degenerate point. Based on the theoretical and experimental methods, Semler and Paidoussis [31, 32] analyzed the nonlinear oscillations of the cantilevered pipe conveying fluid with a sinusoidally perturbed flow velocity and studied the nonlinear dynamic responses of the cantilevered pipe conveying fluid with a small mass attached at the free end. They found that there exist the jumping phenomena and quasi-periodic and chaotic oscillations. Jin and Zou [33] studied the stability and nonlinear dynamics of the cantilevered pipe conveying fluid with the motion limiting constraints and a linear spring support and analyzed the local behaviors in the neighborhood of a double degenerate point. Ghayesh et al. [34] researched the nonlinear dynamic characteristics of the pipes conveying fluid through considering the lateral and longitudinal displacements and found that the cantilevered pipe will generate the flutter through Hopf bifurcation at the critical velocity. Wang et al. [35] studied the nonlinear oscillations of the pipes conveying fluid under the loose constraint boundary conditions. Askarian et al. [36] investigated the nonlinear oscillations of an extensible cantilevered pipe conveying pulsating flow with a nozzle attached to the end of the pipe.
For high-dimensional nonlinear dynamical systems, due to the existence of the modal interactions, there exist the relationships among several types of internal resonances, which can lead to different nonlinear oscillations [37–40]. Panda and Kar [41] established the nonlinear dynamic model of the pipes conveying fluid under the hinged boundaries at both ends and studied the nonlinear dynamic characteristics of the pipes. Ni et al. [42] simplified the pipes conveying fluid to the constant coefficients gyroscopic system and studied the nonlinear oscillations of the pipes on the nonlinear elastic foundation under 3 : 1 internal resonance. Zhang and Chen [43] used the multiple-scale method to determine the steady-state solutions of a pipe conveying fluid under 2 : 1 internal resonance. Mao et al. [44] investigated the forced oscillations of the pipe conveying fluid with 3 : 1 internal resonance around the bending configuration. Zhang et al. [45] established the nonlinear dynamic model of the cantilevered pipe conveying pulsating fluid under the external harmonic excitation and analyzed the multipulse orbits and chaotic oscillations with 1 : 2 internal resonance through utilizing the energy phase method. Recently, Ding et al. [46] investigated the nonlinear vibration isolation of the pipes conveying fluid through using the quasi-zero stiffness characteristics.
To avoid the damage of the pipes conveying pulsating fluid caused by chaotic oscillations, we study the nonlinear oscillations of the pipes conveying pulsating fluid under 1 : 3 internal resonance. Based on the nonlinear dynamic model of the cantilevered beam, the nonlinear dynamic equations of motion for the cantilevered pipe conveying pulsating fluid under the uniformly distributed harmonic excitation and average equations are obtained through using the combination of the multiple-scale method and Galerkin technique under 1 : 3 internal resonance and principal parametric resonance. We analyze the nonlinear resonant responses and the mode interactions of the cantilevered pipe conveying pulsating fluid. Moreover, numerical simulations are performed to study the multitime periodic and chaotic oscillations of the cantilevered pipe conveying pulsating fluid under the external harmonic excitation. It is found that, in the cantilevered pipe conveying pulsating fluid, there are the multitime nonlinear and chaotic oscillations around the left-mode and the right-mode positions, respectively.
2. Equation of Motion and Perturbing Analysis
We consider a cantilevered pipe conveying pulsating fluid, where A is the internal cross-sectional area, L is the length of a tubular cantilevered beam, m is the mass per unit length of the pipe, EI is the flexural rigidity, M is the conveying fluid mass per unit length with an axial velocity U which may vary with respect to time, and f is an external harmonic force, as shown in Figure 1. It is assumed that initially vertical position along the X-axis is located in the direction of gravity and the nonlinear vibrations occur in the X,Y plane for the cantilevered pipe conveying pulsating fluid.
The dynamic model of a cantilevered pipe conveying pulsating fluid is given.
The basic assumptions of the cantilevered pipe and the fluid are made as follows:
The fluid is incompressible.
The diameter of the pipe is small compared to its length. Therefore, the pipe behaves like an Euler-Bernoulli beam.
The vibration of the pipe is planar, and the deflections of the pipe are large.
The rotatory inertia and shear deformation are neglected.
The pipe centerline has inextensible property in the case of a cantilevered pipe.
Based on researches given by Semler et al. [3], Païdoussis [4], and Zhang et al. [45], the nonlinear partial differential governing equation of motion for the cantilevered pipe conveying pulsating fluid is derived through utilizing Hamilton’s principle under the external harmonic excitation:(1)m+My¨+2MUy˙′1+y′2+m+Mgy′1+12y′2+y″MU21+y′2+MU˙−m+MgL−s1+12y′2+EIy″″+4y′y″y″′+y″3+y″″y′2−y″∫sL∫0sm+My˙′2+y′y¨′dsds+∫sL12MU˙y′2+2MUy′y˙′+MU2y′y″ds+y′∫0sm+My˙′2+y′y¨′ds+f=0.
In this paper, we assume that the cantilevered pipe conveying pulsating fluid is made of Kelvin-Voigt type viscoelastic material. Therefore, we have(2)σ=Eε+Eε˙,E=E1+a∂∂t.
In order to obtain the dimensionless governing equation of motion for the cantilevered pipe conveying pulsating fluid, the transformations of the variables and the parameters are introduced as(3)ξ=sL,η=yL,τ=EIm+M1/2tL2,α=EIm+M1/2aL2,U∗=MEI1/2UL,γ=m+MEIL3g,β=mm+M,f∗=L3EIf.
For simplicity, the notation ∗ is omitted in the following analysis. We obtain the dimensionless partial differential governing equation of motion for the cantilevered pipe conveying pulsating fluid:(4)αη˙″″+η″″+η¨+2Uβη˙′1+η′2+η″U21+η′2+U˙β−γ1−ξ1+32η′2+γη′1+12η′2+1+α∂∂τη″″η′2+4η′η″η″′+η″3−η″∫ξ1∫0ξη′2+η′η¨′dξdξ+∫ξ112U˙βη′2+2Uβη′η˙′2+U2η′η″dξ+η′∫0ξη˙′2+η′η¨′dξ+f=0.
It is assumed that the flow velocity U is represented as a periodic perturbation U=U0+U1cosΩ1t. In the meantime, the external excitation f is described as a periodic perturbation f=F+f1cosΩ2t.
We find that there exist the gyroscopic terms in equation (4), which means that the damping of the cantilevered pipe conveying pulsating fluid is too large to ignore. The traditional methods dealing with the nonlinear oscillations are first discretized by Galerkin method and then are analyzed by the perturbation method. However, these methods usually cancel the coefficients related to the gyroscopic terms and lose a lot of important information. In order to retain the gyroscopic terms as much as possible, we use a combining method of both the multiple-scale method and Galerkin method to obtain the average equation of the cantilevered pipe conveying pulsating fluid.
For the perturbation analysis, equation (4) is simplified and the following expression is given for the cantilevered pipe conveying pulsating fluid:(5)Mη¨+Gη˙+Kη=εFη−εNη−2εU1βη˙′cosΩ1t−2εη″U0U1cosΩ1t−εη″Ω1U1β1−ξsinΩ1t,where ε represents the small-scale transformation and other parameters are given as follows:(6a)M=I,G=2U0β∂∂x+α∂∂x4,K=U02−γ1−ξ∂∂x2+γ∂∂x+∂∂x4,(6b)Fη=−F−f1cosΩ2t,(6c)Nη=2U0βη˙′η′2+η″η′2U02−32γ1−ξ+12γη′3+1+a∂∂τη″″η′2+4η′η″η″′+η″3−η″∫ξ1∫0ξη′2+η′η¨′dξdξ−η″∫ξ12U0βη′η˙′2+U02η′η″dξ+η′∫0ξη˙′2+η′η¨′dξ.
The method of multiple scales [47] is applied to the partial differential equation (5) to obtain the uniform solutions in the following form:(7)ηt,ε=η0T0,T1+εη1T0,T1,where T0=t and T1=εt.
Then, the derivatives with respect to t become(8a)ddt=∂∂T0∂T0∂t+∂∂T1∂T1∂t+⋯=∂∂T0+ε∂∂T1+⋯=D0+εD1+⋯,(8b)d2dt2=D0+εD1+⋯2=D02+2εD0D1+ε2D12+⋯,where D0=∂/∂T0 and D1=∂/∂T1.
Substituting equations (7), (8a), and (8b) into equations (6a)–(6c) and eliminating the secular terms, we obtain the following:
The boundary conditions are given as(11)y0=y′0=0,y″L=y″′L=0.
In order to simplify the cantilevered pipe conveying pulsating fluid to the finite dimension by using Galerkin discretization, the modal function is selected as a beam function:(12)φjξ=coshβjx−cosβjx−λjsinhβjx−sinβjx.
The general solution of equation (9) is expressed as follows:(13)ηT0,T1,ξ=∑j=0NφjξAjT1eiωjT0+cc,where we take the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid into account in the first two oscillation modes:(14)ηT0,T1,x=φ1xA1T1eiω1T0+φ3xA2T1eiω3T0+cc.
Substituting equation (14) into equation (5) and after a series of calculations and simplification, the approximate frequency relationships are obtained as follows:(15a)ω12=α12−α11L˜2γ+α11UU+L˜β2,(15b)ω32=α22−α21L˜2γ+α21UU+L˜β2,where L˜, β, and U, respectively, represent the length, mass ratio, and flow velocity, and α11, α12, α21, α22, and γ are the system parameters and dimensionless parameters.
According to the geometries and the material properties of the cantilevered pipe conveying pulsating fluid, the parameters are chosen as β=0.9, U=3, α11=−1.9, α12=11, α21=−4.8, α22=30, and γ=2. The length L˜ is selected as the variable. The first-order and third-order natural frequencies are solved, as shown in the Campbell diagram in Figure 2. We find that the natural frequencies of the first-order and third-order oscillation modes for the cantilevered pipe conveying fluid decrease as the length increases. It is observed that 1 : 3 internal resonance happens in the cantilevered pipe conveying pulsating fluid.
The Campbell diagram of the cantilevered pipe conveying pulsating fluid is obtained in the case of 1 : 3 internal resonance.
We only consider the case of 1 : 3 internal resonance, principle parameter resonance, and 1/2 subharmonic resonance for equation (5). In this resonant case, there exist the following relations:(16)ω1=12Ω1+εσ1,ω3=32Ω1+εσ2,Ω1=12Ω2,where σ1 and σ2 are two detuning parameters.
Substituting equations (12), (14), and (16) into equation (10) yields(17)MD02η1+GD0η1+Kη1=Γ1eiT01/2Ω1+εσ1+Γ2eiT03/2Ω1+εσ2−q1e1/2iT0Ω1+Γ3eiT01/2Ω1−εσ1+Γ4eiT01/2Ω1+εσ2+Γ5eiT03/2Ω1+εσ1+cc+NST,where(18a)Γ1=iΩA112A12U0βm1+32αm2+iΩA11A21A222U0βm3+32αm4+A112A123U02m5−γm6−3m7+34Ω2m8+A11A21A22U02m9−γm10+m4+Ω2m11−iΩ1A˙11Mφ1−A˙11Gφ1,(18b)Γ2=iΩA212A223U0βn1+αn2+iΩA11A12A216U0βn3+αn4+A212A223U02n5−γn6+n7+Ω2n8+A11A12A21U02n9−γn10+n4+Ω2n11−3iΩ1A˙21Mφ2−A˙21Gφ2,(18c)Γ3=A11A122−12Ω12m12−iΩ1U0βm13+U02m14+iΩA11A12A216U0βn3+αn4+A11A21A22−Ω12m15−iΩ1U0βm16+U02m17+12iU1βA12Ω1m18−U0U1A12m19,(18d)Γ4=12iA21Ω1U1βm20−A21U0U1m21,(18e)Γ5=12iA11Ω1U1βn12−A11U0U1n13.
The solution exists when the nonhomogeneous equation corresponding to equation (17) satisfies that the right side of the equation is orthogonal to the solution of its homogeneous accompanying equation. Thus, we have(19a)∫01Γ1−f1e−iσ1T1+Γ3e−2iσ1T1+Γ4e−iT1σ1−σ2φ¯1dx=0,(19b)∫01Γ2+Γ5e−iT1σ2−σ1φ¯2dx=0.
Let A1 and A2 be of the following forms:(20)A1T1=12a1eiφ1⋅e−iσ1T1−2mπ,A2T1=12a2eiφ2⋅e−iσ2T1−2nπ.
Substituting equation (20) into equations (19a) and (19b) and separating the real and imaginary parts, the averaged equations in the polar form are obtained as follows:(21a)a˙1=18β14+β111−β10cos2φ1+β16sin2φ1a1a22+12β17sin2φ1−β18cos2φ1a1+18β13+β110a13+12β19sinφ1−φ2−β20cosφ1−φ2a2−Fsinφ1−μa1,(21b)a1φ˙1=18β12+β10sin2φ1+β16cos2φ1a1a22+12β17cos2φ1+β18sin2φ1a1+18β11−β15a13+12β19cosφ1−φ2+β20sinφ1−φ2a2−Fcosφ1−σ1a1,(21c)a˙2=18β23a23+18β24a2a12+12β25sinφ1−φ2+β26cosφ1−φ2a1−μa2,(21d)a2φ˙2=18β21a23+18β22a2a12+12β25cosφ1−φ2+β26sinφ1−φ2a1−σ2a2.
In order to obtain the averaged equations in the Cartesian form, we express A1 and A2 in the following forms:(22)A1T1=x1T1−ix2T1e−iο1T1−2mπ,A2T1=x3T1−ix4T1e−iο2T1−2nπ.
Based on the same way as the aforementioned analysis, the averaged equations in the Cartesian form are obtained for the cantilevered pipe conveying pulsating fluid:(23a)x˙1=β17−μx1−β18+σ1x1+β19x3+β20x4+β11−β15x1x12+x22+β12−β16x1x32+x42−β13+β110−p15x2x12+x22−β14+β111−2p17x2x32+x42,(23b)x˙2=β17−μx2−β18−σ1x1−β19x4+β20x3+β11+β15x2x12+x22+β12+β16x2x32+x42−β13−β110−p15x1x12+x22−β14−β111−2p17x1x32+x42−Fp0,(23c)x˙3=−μx3−σ2x4+β25x1+β26x2+β21x3x32+x42+β22x3x12+x22−β23x4x32+x42−β24x4x12+x22,(23d)x˙4=−μx4+σ2x3−β25x2+β26x1+β21x4x32+x42+β22x4x12+x22+β23x3x32+x42+β24x3x12+x22,where the coefficients are given, respectively, as(24)β10=Ω1U0βp16,β11=U0βp1+32αp2,β12=2U0βp3+32αp4,β13=3U02p5−γp6+3p7+34Ω2p8,β14=U02p9−γp10+p4+Ω2p11,β15=Ω1U0βp13,β16=U02p17−Ω2p15,β17=12U1Ω1βp18,β18=U0U1p19,β19=12Ω1U1βp20,β110=U02p4−12Ω2p12,β111=U02p17−Ω2p15,β20=U0U1p21,β21=Ω3U0βq1+αq2,β22=Ω6U0βq3+αq4,β23=3U02q5−γq6+q7+Ω2q8,β24=U02q9−γq10+q4+Ω2q11,β25=12Ω1U1βq12,β26=U0U1q13.
Equations (21a)–(21d) and (23a)–(23d) describe the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid under the external harmonic excitation. It is known that the constant solutions of the averaged equation correspond to the periodic solutions of the original system, the periodic solutions of the averaged equation correspond to the quasi-periodic solutions of the original system, and the chaotic solutions of the averaged equation correspond to the chaotic solutions of the original system.
3. Analysis of Resonant Responses
The resonant response curves are the important basis for judging and studying the nonlinear oscillations. They include many complex nonlinear dynamic phenomena. The practical problems can be solved better by analyzing these phenomena. Based on the averaged equations (21a)–(21d) in the polar form, we analyze the amplitude-frequency response curves to reveal the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid. Since φ1 and φ2 are all constant, we make φ1=φ2=0.25π. Therefore, we obtain the amplitude-frequency response functions of the cantilevered pipe conveying pulsating fluid:(25a)18β13+β110a13+18β14+β111a1a22−12β20a2−22F−μa12+18β11+β15a13+18β12a1a22+12β19−β20a2−22F−σ1a12=18β110a1a22−12β18a12+18β16a1a22+12β17a12,(25b)18β24a2a12+18β23a232+18β22a2a12+18β21a232=12β26a1+μa22+12β25a1+σ2a22.
Two cases of the amplitude-frequency response curves [48, 49] are considered for the cantilevered pipe conveying pulsating fluid:
There is no coupling effect between the first-order and third-order oscillation modes. We only consider the decoupled case and set a2=0 in equation (25a) and a1=0 in equation (25b).
We let an amplitude change when another is fixed because the amplitude is much smaller than that in the case of weak coupling; namely, set a2=1 in equation (25a) and a1=1 in equation (25b).
We obtain the amplitude-frequency response curves in the decoupled and coupled cases based on equations (25a) and (25b). Based on the geometries and the material properties of the cantilevered pipe conveying pulsating fluid, the basic parameters are chosen as(26)U0=1.5,F=2.5,β=0.9,γ=2,μ=4,U1=3.5,Ω=1.2,Ω1=10,p1=2×103,p2=4,p3=2,p4=8,p5=10,p6=2,p7=−20,p8=1,p9=100,p10=−1,p11=−3,p12=−14,p13=−6,p14=−10,p15=14,p16=10,p17=4,p18=−1,p19=20,p20=30,p21=12,p22=1,p23=1,p24=6,q1=15,q2=16,q3=−4.7,q4=18,q5=18,q6=7,q7=4,q8=14,q9=11,q10=11,q11=260,q12=3,q13=0.25.
For both the decoupled and coupled cases, when the velocity parametric excitation is U0=1.5 and external excitation is F=2.5, it is found that the cantilevered pipe conveying pulsating fluid has the hardening spring characteristics because the amplitude-frequency response curves bend to right, as shown in Figures 3(a) and 3(b). The blue line represents the coupled case and the red line gives the decoupled case. The bending of the amplitude-frequency response curves also leads to the multiple amplitudes and bifurcations. The multiple amplitudes mean that there are several possible responses in the cantilevered pipe conveying pulsating fluid. The actual response among several possible responses depends on the initial conditions. Hence, the typical jumping phenomena of the nonlinear oscillations can be observed for the cantilevered pipe conveying pulsating fluid. Comparing the amplitude-frequency response curves between two cases given in Figure 3(a), we find that the amplitude of the coupled case is larger and the resonant response interval is wider than the decoupled case. The coupling effect of two oscillation modes will not change the system’s stiffness but will cause a stronger nonlinear resonant response. Furthermore, it can be also observed that the amplitude of the first-order mode is larger than that of the third-order mode in the coupled case, as shown in Figure 3(b). In the coupled case, there exists the mode interaction between two oscillation modes. Therefore, the amplitude-frequency response curves of the coupled case are considered for the cantilevered pipe conveying pulsating fluid in the following analysis.
The amplitude-frequency response curves of the cantilevered pipe conveying pulsating fluid are obtained when flow velocity is U0=1.5 and external force is F=2.5. (a) Coupling influence on amplitude-frequency response curves; (b) amplitude-frequency response curves of the first-order and third-order modes.
In the weak coupling cases, the influences of the flow velocity U0, external excitation F, and mass parameter β on the amplitude-frequency response curves are analyzed in Figures 4(a)–4(c), where the mass parameter β denotes the ratio of the pipe mass and total mass. Figure 4(a) illustrates the amplitude-frequency response curves of the cantilevered pipe conveying pulsating fluid when the flow velocities are U0=1.5 and U0=2.0, respectively. With the increases of the flow velocity U0, it is obviously found that the amplitude is larger and the resonant response interval is wider for the cantilevered pipe conveying pulsating fluid, as shown in Figure 4(a). The increase of flow velocity will not change the system’s stiffness but will cause a stronger nonlinear resonant response. These characteristics indicate that the cantilevered pipe conveying pulsating fluid obtains the more external energy. Figure 4(b) gives the amplitude-frequency response curves of the cantilevered pipe conveying pulsating fluid when the external excitations are F=2.5 and F=7.0, respectively. In Figure 4(b), we find that the increase of the external excitation F mainly results in the larger amplitude, which has little effect on the resonant response interval. Figure 4(c) shows the amplitude-frequency response curves of the cantilevered pipe conveying pulsating fluid when the mass parameters are β=0.9 and β=0.7, respectively. In Figure 4(c), as the mass parameter β decreases, which means that the pipe mass is lighter, we find that the nonlinear stiffness of the cantilevered pipe conveying pulsating fluid does not change. However, it is observed that the resonance of the cantilevered pipe conveying pulsating fluid is more easily excited and the amplitudes of the nonlinear oscillations are greater when mass parameters are smaller.
The amplitude-frequency response curves of the cantilevered pipe conveying pulsating fluid are obtained in the coupling case with different parameters. (a) Amplitude-frequency response curves when flow velocities U0=1.5 and U0=2.0, (b) amplitude-frequency response curves when external forces F=2.5 and F=7.0, and (c) amplitude-frequency response curves when mass parameters β=0.9 and β=0.7.
The effects of the flow velocity and external excitation on the amplitude-frequency response curves are investigated for the cantilevered pipe conveying pulsating fluid in Figures 5(a) and 5(b) under the coupling cases. We set σ=10 and other parameters are the same as the parameters given in Figures 3 and 4. The amplitude-frequency response curves of the flow velocity are given in Figure 5(a). It can be shown that the amplitudes of the cantilevered pipe conveying pulsating fluid increase with the increase of the flow velocity. Thus, the energy of the fluid motion is transferred into the energy of the cantilevered pipe. The fluid-structure interaction occurs in the cantilevered pipe conveying pulsating fluid. There are two jumping phenomena on the amplitude-frequency response curves with the increase of the flow velocity: one is jumping up and one is jumping down. The oscillations of the cantilevered pipe conveying pulsating fluid are stable at a lower amplitude. The increase and decrease of the amplitudes alternatively change with the increase of the flow velocity. It is indicated that the flow velocity will affect the nonlinear oscillation behaviors of the cantilevered pipe conveying pulsating fluid in a certain range. The amplitude-frequency response curves of the external excitation are given in Figure 5(b). The amplitudes of the cantilevered pipe conveying pulsating fluid become larger with the gradual increase of the external excitation. The sudden increase in the amplitude of the cantilevered pipe conveying pulsating fluid appears when a jumping phenomenon occurs. It is illustrated that the cantilevered pipe conveying pulsating fluid will obtain more energy to create the larger amplitude with the increase of the external excitation.
The amplitude-frequency response curves of the cantilevered pipe conveying pulsating fluid are obtained. (a) Amplitude-response curves on flow velocity; (b) amplitude-response curves on external excitation.
We explore the influences of different parameters on the amplitude-frequency and force-amplitude response curves in the cantilevered pipe conveying pulsating fluid, respectively. In Figures 6(a)–6(d), the changes of the amplitude-frequency and force-frequency response curves in each figure are caused by changing one parameter while other parameters remain unchanged. The influences of different parameters, such as the parametric excitation and mass parameter, on the amplitude-frequency response curves are studied, as shown in Figures 6(a) and 6(b). Figure 6(a) gives the relations between the amplitudes and flow velocity U0 for the cantilevered pipe conveying pulsating fluid when the parametric excitations induced to the flow velocity are F=5.0 and F=7.0, respectively. Figure 6(b) indicates the relations on the amplitudes versus the flow velocity U0 in the cantilevered pipe conveying pulsating fluid when the mass parameters are β=0.9 and β=0.7, respectively. We find that the amplitudes of the nonlinear oscillations for the cantilevered pipe conveying pulsating fluid become larger as the parametric excitation increases and the mass parameter decreases. This indicates that when the pipe mass is lighter, the amplitudes of the nonlinear oscillations are greater. In Figures 6(c) and 6(d), we study the influences of the external excitations on the force-amplitude response curves. Figure 6(c) demonstrates the relations between the amplitudes versus the external excitations F in the cantilevered pipe conveying pulsating fluid when the flow velocities are U0=15.8 and U0=16.0, respectively. We find that the amplitudes of the nonlinear oscillations for the cantilevered pipe conveying pulsating fluid become larger when the flow velocity is larger. Figure 6(d) shows the relations on the amplitudes versus the external excitations F in the cantilevered pipe conveying pulsating fluid when the mass parameters are β=0.9 and β=0.7, respectively. It is found that the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid more easily happen with the decrease of the mass parameter β.
The influences of different parameters on the amplitude-frequency response curves and force-amplitude response curves are given for the cantilevered pipe conveying pulsating fluid. (a) Amplitude-frequency response curves on flow velocity when F=5 and F=7, (b) amplitude-frequency response curves on flow velocity when mass parameters β=0.9 and β=0.5, (c) force-amplitude response curves on external force when flow velocities U0=15.8 and U0=16, and (d) force-amplitude response curves on external force when mass parameters β=0.9 and β=0.7.
Based on the analyses of the amplitude-frequency response curves, it can be found that the hardening-spring type behaviors and jumping phenomena are exhibited for the cantilevered pipe conveying pulsating fluid. The jumping phenomena also occur in the force-amplitude response curves versus the flow velocity and external excitation. Moreover, it is known that the flow velocity, external excitation, and coupling degree of two oscillation modes can affect the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid under the external harmonic force. The nonlinearity bends the amplitude-frequency response curves to the right when the cantilevered pipe conveying pulsating fluid has the hardening-spring type behaviors. The bending of the amplitude-frequency response curves leads to the occurrence of the jumping phenomena. The hardening-spring type behavior and jumping phenomena mean that the oscillation amplitudes of the cantilevered pipe conveying pulsating fluid increase firstly and change abruptly. The energy of the fluid motion is transferred into the energy of the cantilevered pipe when the flow velocity increases. Thus, the fluid-structure interaction happens for the cantilevered pipe conveying pulsating fluid. The nonlinear resonance of the large amplitude for the cantilevered pipe conveying pulsating fluid can be stimulated with the increase of these factors.
4. Numerical Simulations of Periodic and Chaotic Oscillations
In order to study the nonlinear dynamic properties of the cantilevered pipe conveying pulsating fluid, the influences of the velocity parametric excitation and external excitation on the nonlinear oscillations of the pipe are investigated. In this section, the average equations (23a)–(23d) in the Cartesian coordinate system are numerically simulated to find the multitime periodic and the chaotic oscillations of the cantilevered pipe conveying pulsating fluid under the external harmonic excitation based on the fourth-order Runge-Kutta algorithm [50–52]. We choose the velocity parametric and external excitations as the controlling parameters to study the complicated nonlinear dynamics of the cantilevered pipe conveying pulsating fluid, respectively. The initial conditions are chosen as x1=0, x2=0.1, x3=0, and x4=0.1. Other parameters are chosen as(27)β=0.9,γ=7,μ=15,σ=7.7,Ω=15,Ω1=8.6,α=0.7,p0=−3,p1=0.5,p2=−11,p3=2,p4=8,p5=−10,p6=−13,p7=16,p8=−10,p9=−1,p10=−10,p11=−3,p12=−14,p13=−6,p14=−1,p15=14,p16=10,p17=4,p18=−9,p19=2,p20=1,p21=4,p22=1,p23=1,p24=6,q1=−5,q2=−6,q3=−4.7,q4=18,q5=18,q6=7,q7=4,q8=14,q9=−1,q10=−11,q11=9,q12=−10,q13=3.
To reveal the nonlinear dynamic behaviors of the cantilevered pipe conveying pulsating fluid, the bifurcation diagrams, maximum Lyapunov exponents, phase portraits, waveforms, and Poincare map are depicted. Figures 7 and 8 present the bifurcation diagrams and maximum Lyapunov exponents to demonstrate the multitime nonlinear oscillation of the cantilevered pipe conveying pulsating fluid by varying the velocity parametric excitation and external excitation, respectively.
The bifurcation diagrams and maximum Lyapunov exponents of the cantilevered pipe conveying pulsating fluid are given for the parametric excitations of the flow velocities when the external excitation is F=7. (a) x1 versus flow velocity, (b) x3 versus flow velocity, and (c) maximum Lyapunov exponents versus flow velocity.
The bifurcation diagrams and maximum Lyapunov exponents of the cantilevered pipe conveying pulsating fluid are obtained for the external excitations when the flow velocity is U0=12. (a) x1 versus external excitation, (b) x3 versus external excitation, and (c) maximum Lyapunov exponents versus flow velocity.
In order to study the influences of the flow velocities on the multitime nonlinear oscillations of the cantilevered pipe conveying pulsating fluid, we set the external excitation while other parameters remain unchanged. The bifurcation diagrams and maximum Lyapunov exponents of the cantilevered pipe conveying pulsating fluid with the change of the flow velocity are shown in Figure 7 when the external excitation F=7. According to changing characteristics of the bifurcation diagram, with the increase of the flow velocities, the period oscillations of the cantilevered pipe conveying pulsating fluid appear at the beginning, the multitime quasi-period oscillations appear, and the multitime chaotic oscillations happen. With the increase of the flow rate, the multiple periodic windows appear for the cantilevered pipe conveying pulsating fluid. This illustrates the alternating change of the periodic and chaotic oscillations. When the flow rate increases to a particular value, the motion of the cantilevered pipe conveying pulsating fluid evolves to a stable period-1 oscillation. The maximum Lyapunov exponents given in Figure 7(c) describe the intervals of the flow velocity for the chaotic vibrations.
Based on numerical simulations corresponding to the aforementioned analyses, we further investigate the effects of the external excitations on the nonlinear dynamics of the cantilevered pipe conveying pulsating fluid when the flow velocity is U0=12 and other parameters remain unchanged, as shown in Figure 8. It is found that the changing characteristics of the oscillations with the increase of the external excitations are similar to those with the increase of the flow velocities for the cantilevered pipe conveying pulsating fluid. The periodic oscillation appears at the beginning, and then the multitime periodic oscillation appears, and, finally, the multitime chaotic oscillations occur for the cantilevered pipe conveying pulsating fluid under the external harmonic excitation. With the further increase of the external excitation, several periodic oscillation windows appear for the cantilevered pipe conveying pulsating fluid. When the external excitation increases to a particular value, the motions of the cantilevered pipe conveying pulsating fluid eventually become the chaotic oscillations from the intermittent chaos. The intermittent chaos is verified in Figure 8(c) by whether the maximum Lyapunov exponent is positive.
We present a variety of figures to confirm the vibrations of the cantilevered pipe conveying pulsating fluid corresponding to different flow rates and different external excitations. In Figures 9–17, Figures(a) and (c) represent the phase portraits on the planes x1,x2 and x3,x4, Figures (b) and (d) depict the waveforms on the planes t,x1 and t,x3, Figure(e) denotes three-dimensional phase portrait in the space x1,x2,x3, and Figure(g) depicts Poincare map on the plane x1,x2, respectively. The aforementioned figures are the effective tools to judge the chaotic oscillations of the cantilevered pipe conveying pulsating fluid. Whether the waveforms are regular or irregular, these indicate that the motions of the cantilevered pipe conveying pulsating fluid are the periodic or chaotic oscillations. The phase diagram is the projection curve of the solution in the phase plane and phase space. The closed curve represents the periodic oscillation, and the phase portraits of the chaotic oscillations are the repeating unclosed curve confined within a bounded region. According to Poincare mapping, the periodic and quasi-periodic oscillations are the isolated points and curves on Poincare cross section. The nonlinear oscillations of the cantilevered pipe conveying pulsating fluid are the multitime chaotic oscillations when the discrete point sets appear with the infinite self-similar structures on Poincare cross section.
The periodic oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the parametric excitation of the flow velocity is U0=11.0. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The period-2 oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the parametric excitation of the flow velocity is U0=19.5. (a) phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The multiperiodic oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the parametric excitation of the flow velocity is U0=23.2. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The chaotic oscillations of the cantilevered pipe conveying pulsating fluid are obtained when the parametric excitation of the flow velocity is U0=31.5. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The period-1 oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the parametric excitation of the flow velocity is U0=42.5. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The multiperiodic oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the external excitation is F=2.9. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The multitime periodic oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the external excitation is F=4.5. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The quasi-periodic oscillation of the cantilevered pipe conveying pulsating fluid is obtained when the external excitation is F=6.3. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
The chaotic oscillations of the cantilevered pipe conveying pulsating fluid are obtained when the external excitation is F=6.5. (a) Phase portrait on plane x1,x2, (b) waveform on plane t,x1, (c) phase portrait on plane x3,x4, (d) waveform on plane t,x3, (e) three-dimensional phase portrait in space x1,x2,x3, and (g) Poincare map on plane x1,x2.
According to the characteristics of each figure for different motions, there exist the periodic, multitime quasi-periodic, and multitime chaotic oscillations of the cantilevered pipe conveying pulsating fluid under different velocity parametric excitations, as shown in Figures 9–13. Figure 9 gives the periodic oscillation of the cantilevered pipe conveying pulsating fluid when the velocity parametric excitation is U0=11.0. Figure 10 indicates that there exists the period-2 oscillation of the cantilevered pipe conveying pulsating fluid when the velocity parametric excitation is U0=19.5. Based on the definition of the coordinate system, we know that the right is positive and the left is negative for the cantilevered pipe conveying pulsating fluid. It is found that, in the cantilevered pipe conveying pulsating fluid, there are the nonlinear oscillations around the left-mode and right-mode positions, respectively. The first-order and third-order oscillation modes of the cantilevered pipe conveying pulsating fluid vibrate twice times around the right-mode position and then vibrate twice times around the left-mode position, as shown in Figures 10(b) and 10(d). Between the left-mode and right-mode nonlinear oscillations, there exists one-time oscillation around the vertical equilibrium position in the cantilevered pipe conveying pulsating fluid. This is a typical nonlinear dynamic phenomenon.
Figure 11 illustrates that the multiperiodic oscillation happens when the velocity parametric excitation is U0=23.2. From Figure 11(b), it is demonstrated that, in the first-order mode of the nonlinear oscillations, there exists one-time oscillation around the right-mode position and then there exists twice-time oscillation around the left-mode position. For the third-order mode, there exists twice-time oscillation around the right-mode position and then there exists one-time oscillation around the left-mode position in the cantilevered pipe conveying pulsating fluid, as shown in Figure 11(d).
Figure 12 indicates that there exist the multitime chaotic oscillations when the velocity parametric excitation is U0=31.5. Around the left-mode and right-mode positions, the chaotic oscillations of the cantilevered pipe conveying pulsating fluid are the irregular motions. Figure 13 indicates that the motion of the cantilevered pipe conveying pulsating fluid again returns to the period-1 oscillation when the velocity parametric excitation is U0=42.5. In this case, it is found from Figure 13(b) that the first-order mode of the cantilevered pipe conveying pulsating fluid is the oscillation around the left-mode position. However, the third-order vibration mode of the cantilevered pipe conveying pulsating fluid is a motion around the right-mode position, as shown in Figure 13(d).
We also find the periodic, multitime periodic, and multitime chaotic oscillations of the cantilevered pipe conveying pulsating fluid under different external excitations, as shown in Figures 14–17. Figure 14 gives the multiperiodic oscillation of the cantilevered pipe conveying pulsating fluid when the external excitation is F=2.9. It is observed from Figure 14(b) that, for the first-order mode, the motions of the cantilevered pipe conveying pulsating fluid vibrate three times around the right-mode position and twice times around the left-mode position. For the third-order mode, the vibrations of the cantilevered pipe conveying pulsating fluid move twice times around the right-mode position and three times around the left-mode position, as shown in Figure 14(d). In addition, it is observed from Figure 14 that the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid also have the motions around the vertical equilibrium position.
Figure 15 indicates that the multitime periodic oscillation happens for the cantilevered pipe conveying pulsating fluid when the external excitation is F=4.5. For the first-order mode, four-time nonlinear oscillations occur around the right-mode position and twice-time nonlinear oscillations occur around the left-mode position in the cantilevered pipe conveying pulsating fluid, as shown in Figure 15(b). It is also found from Figure 15(d) that, for the third-order mode, there are twice-time nonlinear oscillations around the right-mode position and four-time nonlinear oscillations around the left-mode position. Figure 16 represents that the quasi-periodic oscillation occurs when the external excitation is F=6.3. Figure 17 demonstrates that there exist the multitime chaotic oscillations of the cantilevered pipe conveying pulsating fluid when the external excitation is F=6.5. The irregular oscillations around the left-mode and the right-mode positions exist in the cantilevered pipe conveying pulsating fluid under the external harmonic excitation.
Moreover, the oscillations corresponding to different velocity parametric excitations and external excitations are consistent with the bifurcation diagrams of the cantilevered pipe conveying pulsating fluid under the external harmonic force, as shown in Figures 7 and 8. Based on the aforementioned analyses, it can be found that the flow rate and the external excitation can affect the nonlinear oscillations of the cantilevered pipe conveying pulsating fluid under the case of 1 : 3 internal resonance. When the velocity parametric excitation or the external excitation reaches a certain value, the multitime chaotic oscillations will happen in the cantilevered pipe conveying pulsating fluid, and the disorderly motions will lead to the destruction of the cantilevered pipe conveying pulsating fluid under the external harmonic force. Therefore, the velocity parametric excitation and the external excitation causing the multitime chaotic oscillations should be controlled when engineers design and use the infusion tube system.
5. Conclusions
The nonlinear resonant responses and multitime chaotic dynamics of the cantilevered pipe conveying pulsating fluid are investigated under the external harmonic force. Based on the nonlinear partial differential governing equation of motion for the cantilevered pipe conveying pulsating fluid derived by using Hamilton’s principle, the 1 : 3 internal resonance and primary parametric resonance-1/2 subharmonic resonance are considered. A combination method of the method of multiple scales and Galerkin technique is utilized to obtain four-dimensional nonlinear averaged equations. Several amplitude-frequency response curves are obtained corresponding to the certain parameters. From the analysis of the amplitude-frequency response curves, it is found that there exist the hardening-spring type behaviors and the jumping phenomena. The jumping phenomena also occur in the amplitude-force response curves versus the flow velocity and external force.
Moreover, we find that the flow velocity, external force, and coupling degree of two order modes can affect the nonlinear vibrations of the cantilevered pipe conveying pulsating fluid under the external harmonic force. The nonlinear vibrations of the cantilevered pipe conveying pulsating fluid can be excited more easily with the increase of the flow velocity, external force, and coupling degree of two order modes. It is known that the nonlinear dynamic behaviors of the cantilevered pipe conveying pulsating fluid under the external harmonic force will be affected due to the flow rate and external excitation under the case of 1 : 3 internal resonance. It is observed that the multitime chaotic vibrations will occur for the cantilevered pipe conveying pulsating fluid when the velocity parametric excitation or external excitation reaches a certain value. From Figures 9–17, it is found that, for the first-order and third-order vibration modes, several-time nonlinear vibrations simultaneously or respectively happen around the right-mode and left-mode positions in the cantilevered pipe conveying pulsating fluid. Therefore, the velocity parametric excitation and external excitation should be controlled when engineers design and use the infusion tube systems.
Data Availability
All data generated or analyzed during this study are included in this published article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors gratefully acknowledge the support of the National Natural Science Foundation of China (NNSFC) through Grant nos. 11672188, 11832002, and 11427801 and the Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (PHRIHLB).
TangD. M.DowellE. H.Chaotic oscillations of a cantilevered pipe conveying fluid1988226328310.1016/s0889-9746(88)80011-22-s2.0-38249029816PaïdoussisM. P.MoonF. C.Nonlinear and chaotic fluid elastic vibrations of a flexible pipe conveying fluid1988256759110.1016/s0889-9746(88)80023-92-s2.0-0023704905SemlerC.LiG. X.PaïdoussisM. P.The non-linear equations of motion of pipes conveying fluid19945577599PaïdoussisM. P.20142ndLondon, UKAcademic PressHolmesP. J.Bifurcations to divergence and flutter in flow-induced oscillations: a finite dimensional analysis19775347150310.1016/0022-460x(77)90521-12-s2.0-0017522897HolmesP. J.MarsdenJ.Bifurcation to divergence and flutter in flow-induced oscillations: an infinite dimensional analysis19781436738410.1016/0005-1098(78)90036-52-s2.0-0017988271HuoY. L.WangZ. M.Dynamic analysis of a vertically deploying retracting cantilevered pipe conveying fluid201636022423810.1016/j.jsv.2015.09.0142-s2.0-84944339204HosseiniM.BahaadiniR.Size dependent stability analysis of cantilever micro-pipes conveying fluid based on modified strain gradient theory2016101113AskarianA. R.AbtahiH.Firouz-AbadiR. D.HaddadpourH.DowellE. H.Bending-torsional instability of a viscoelastic cantilevered pipe conveying pulsating fluid with an inclined terminal nozzle2018322999300810.1007/s12206-018-0603-02-s2.0-85050120772TubaldiE.AmabiliM.PaïdoussisM. P.Fluid-structure interaction for nonlinear response of shells conveying pulsatile flow201637125227610.1016/j.jsv.2016.01.0242-s2.0-84969348273BaiY. C.XieW. D.GaoX. F.XuW. H.Dynamic analysis of a cantilevered pipe conveying fluid with density variation20188163865510.1016/j.jfluidstructs.2018.06.0052-s2.0-85048835196BajajA. K.SethnaP. R.LundgrenT. S.Hopf bifurcation phenomena in tubes carrying a fluid19803921323010.1137/01390192-s2.0-0019072146Sri NamchchivayaN.Non-linear dynamics of supported pipe conveying pulsating fluid-I, subharmonic resonance19892418519610.1016/0020-7462(89)90037-12-s2.0-0024884789Sri NamchchivayaN.TienW. M.Non-linear dynamics of supported pipe conveying pulsating fluid-II, combination resonance19892419720810.1016/0020-7462(89)90038-32-s2.0-0024928193SarkarA.PaïdoussisM. P.A cantilever conveying fluid: coherent modes versus beam modes20043946748110.1016/s0020-7462(02)00213-52-s2.0-0038148079McDonaldR. J.Sri NamachchivayaN.Pipes conveying pulsating fluid near a 0 : 1 resonance: local bifurcations20052162966410.1016/j.jfluidstructs.2005.07.0192-s2.0-28444446268YoshizawaM.NaoH.HasegawaE.TsujiokaY.Lateral vibration of a flexible pipe conveying fluid with pulsating flow20082922432250HouY.ZengG. H.Research on nonlinear dynamic characteristics of fluid-conveying pipes system2011228-22957457910.4028/www.scientific.net/amr.228-229.5742-s2.0-79956125346SetoodehA. R.AfrahimS.Nonlinear dynamic analysis of FG micro-pipes conveying fluid based on strain gradient theory201411612813510.1016/j.compstruct.2014.05.0132-s2.0-84904539048KheiriM.PaïdoussisM. P.On the use of generalized Hamilton’s principle for the derivation of the equation of motion of a pipe conveying fluid201450182410.1016/j.jfluidstructs.2014.06.0072-s2.0-84907646916GanC. B.GuoS. Q.LeiH.YangS. X.Random uncertainty modeling and vibration analysis of a straight pipe conveying fluid20147750351910.1007/s11071-014-1313-52-s2.0-84905086645ZhangT.OuyangH.ZhangY. O.LvB. L.Nonlinear dynamics of straight fluid-conveying pipes with general boundary conditions and additional springs and masses2016407880790010.1016/j.apm.2016.03.0502-s2.0-84977580725WangZ. M.LiuY. Z.Transverse vibration of pipe conveying fluid made of functionally graded materials using a symplectic method201629814915910.1016/j.nucengdes.2015.12.0072-s2.0-84953727377LiangF.YangX.-D.QianY.-J.ZhangW.Free vibration analysis of pipes conveying fluid based on linear and nonlinear complex modes approach201798175011210.1142/s17588251175011252-s2.0-85039155826LiangF.YangX.-D.ZhangW.QianY.-J.Nonlinear free vibration of spinning viscoelastic pipes conveying fluid2018107185007610.1142/s175882511850076x2-s2.0-85052999417LiangF.YangX. D.QianY. J.ZhangW.Transverse free vibration and stability analysis of spinning pipes conveying fluid201813719520410.1016/j.ijmecsci.2018.01.0152-s2.0-85041397100LiangF.YangX. D.ZhangW.QianY. J.Dynamical modeling and free vibration analysis of spinning pipes conveying fluid with axial deployment2018417657910.1016/j.jsv.2017.12.0052-s2.0-85041541453LiangF.YangX. D.ZhangW.QianY. J.Coupled bi-flexural-torsional vibration of fluid-conveying pipes spinning about an eccentric axis201919195000310.1142/s02194554195000322-s2.0-85052930662LiangF.YangX. D.ZhangW.QianY. J.MelnikR. V. N.Parametric vibration analysis of pipes conveying fluid by nonlinear normal modes and a numerical iterative approach201911385210.4208/aamm.oa-2018-00642-s2.0-85069199064LiG. X.PaïdoussisM. P.Stability, double degeneracy and chaos in cantilevered pipes conveying fluid1994298310710.1016/0020-7462(94)90054-x2-s2.0-0028336585SemlerC.PaïdoussisM. P.Nonlinear analysis of the parametric resonances of a planar fluid-conveying cantilevered pipe19961078782510.1006/jfls.1996.00532-s2.0-13544255708PaïdoussisM. P.SemlerC.Non-linear dynamics of a fluid conveying cantilevered pipe with a small mass attached at the free end199833153210.1016/s0020-7462(97)00002-4JinJ. D.ZouG. S.Bifurcations and chaotic motions in the autonomous system of a restrained pipe conveying fluid200326078380510.1016/s0022-460x(02)00982-32-s2.0-0037422147GhayeshM. H.PaïdoussisM. P.AmabiliM.Nonlinear dynamics of cantilevered extensible pipes conveying fluid20133326405641810.1016/j.jsv.2013.06.0262-s2.0-84883750937WangL.LiuZ. Y.AbdelkefiA.WangY. K.DaiH. L.Nonlinear dynamics of cantilevered pipes conveying fluid: towards a further understanding of the effect of loose constraints20179510.1016/j.ijnonlinmec.2017.05.0122-s2.0-85020539656AskarianA. R.HaddadpourH.FirouzabadiR. D.AbtahiH.Nonlinear dynamics of extensible viscoelastic cantilevered pipes conveying pulsatile flow with an end nozzle201791223510.1016/j.ijnonlinmec.2017.02.0032-s2.0-85013057015LiuT.ZhangW.WangJ. F.Nonlinear dynamics of composite laminated circular cylindrical shell clamped along a generatrix and with membranes at both ends2017901393141710.1007/s11071-017-3734-42-s2.0-85027870110ZhangW.LiuT.XiA.WangY. N.Resonant responses and chaotic dynamics of composite laminated circular cylindrical shell with membranes2018423659910.1016/j.jsv.2018.02.0492-s2.0-85050662973LiuT.ZhangW.MaoJ. J.ZhengY.Nonlinear breathing vibrations of eccentric rotating composite laminated circular cylindrical shell subjected to temperature, rotating speed and external excitations201912746349810.1016/j.ymssp.2019.02.0612-s2.0-85063113370ZhangW.ZhengY.LiuT.GuoX. Y.Multi-pulse jumping double-parameter chaotic dynamics of eccentric rotating ring truss antenna under combined parametric and external excitations20199876180010.1007/s11071-019-05227-82-s2.0-85072174121PandaL. N.KarR. C.Nonlinear dynamics of a pipe conveying pulsating fluid with parametric and internal resonances20074993010.1007/s11071-006-9100-62-s2.0-34248356741NiQ.TangM.LuoY.WangY.WangL.Internal-external resonance of a curved pipe conveying fluid resting on a nonlinear elastic foundation20147686788610.1007/s11071-013-1175-22-s2.0-84899103881ZhangY. L.ChenL. Q.Steady-state response of pipes conveying pulsating fluid near a 2 : 1 internal resonance in the supercritical regime201461450056MaoX. Y.DingH.ChenL. Q.Steady-state response of a fluid-conveying pipe with 3 : 1 internal resonance in supercritical regime20168679580910.1007/s11071-016-2924-92-s2.0-84976870588ZhangY. F.YaoM. H.ZhangW.WenB. C.Dynamical modeling and multi-pulse chaotic dynamics of cantilevered pipe conveying pulsating fluid in parametric resonance20176810.1016/j.ast.2017.05.0272-s2.0-85020935239DingH.JiJ. C.ChenL. Q.Nonlinear vibration isolation for fluid-conveying pipes using quasi-zero stiffness characteristics201912167568810.1016/j.ymssp.2018.11.0572-s2.0-85057580309NayfehA. H.MookD. T.1979Oxford, UKOxford University PressChenJ. E.ZhangW.GuoX. Y.SunM.Theoretical and experimental studies on nonlinear oscillations of symmetric cross-ply composite laminated plates2013731697171410.1007/s11071-013-0896-62-s2.0-84880923091ChenJ. E.ZhangW.LiuJ.SunM.Dynamic properties of truss core sandwich plate with tetrahedral core201513486988210.1016/j.compstruct.2015.08.1092-s2.0-84942515590WeiZ.LiY.SangB.LiuY.ZhangW.Complex dynamical behaviors in a 3D simple chaotic flow with 3D stable or 3D unstable manifolds of a single equilibrium2019297195009510.1142/s02181274195009502-s2.0-85068977176LiuY. J.NazarimehrF.KhalafA. J. M.AhmedA.HayatT.Detecting bifurcation points in a memristive neuron model20192281943195010.1140/epjst/e2019-800205-32-s2.0-85073549911LiuY.KhalafA. J. M.HayatT.AlsaediA.PhamV.-T.JafariS.A complete investigation of the effect of external force on a 3D megastable oscillator2020301205001210.1142/s0218127420500121