This paper investigates the dynamic behavior of a cantilevered microtube conveying fluid, undergoing large motions and subjected to motion-limiting constraints. Based on the modified couple stress theory and the von Kármán relationship, the strain energy of the microtube can be deduced and then the governing equation of motion is derived by using the Hamilton principle. The Galerkin method is applied to produce a set of ordinary differential equations. The effect of the internal material length scale parameter on the critical flow velocity is investigated. By using the projection method, the Hopf bifurcation is demonstrated. The results show that size effect on the vibration properties is significant.
Guizhou Institute of TechnologyXJGC201909091. Introduction
The dynamics of microtube conveying fluid is of considerable interest in engineering and has been studied widely over the past few decades [1–9]. Holmes [10] proved that pipes supported at both ends cannot flutter. The Hopf bifurcation of a cantilevered tube carrying fluid was investigated theoretically by Bajaj et al. [11], revealing that the type of the Hopf bifurcation being supercritical or subcritical depends on the mass ratio and friction factor in a complicated manner. The fluid-conveying tube was modified to include motion-limiting restraints and a full nonlinear analysis was given by Paidoussis and Semler [12]. In [13] by Jin, a linear spring support was attached to the system of Paidoussis et al. [8], and the critical velocity and the chaotic motions of this system were investigated. A variant of the basic system, that is, a standing pipe conveying fluid was studied by Wang and Ni [14], and the dynamics of this variant was also compared with that of a hinging system developed by Jin [13]. By applying the multidimensional Lindstedt–Poincaré method, a clamped-clamped fluid-conveying pipe under external periodic excitation was studied by Liang and Wen [15]. Based on the Galerkin truncation of order 8 and the numerical analysis, the dynamic behaviors of fluid-conveying pipes for three types of supports with lateral motion were discussed by Li et al. [16]. Zhang et al. [17] investigated the nonlinear forced vibration of fluid-conveying pipes in the supercritical regime, that is, vibration about a curved equilibrium.
Due to recent technological development in science and technology, the characteristic size of pipes becomes smaller and smaller. In [18] by Rinaldi et al., the inside diameter of the circular microtube ranges from 1 to 100 μm. However, regarding its application in microelectronic-mechanical-systems and engineering microfluidic devices, as an example, microtube conveying fluid can be applied to a class of microresonators [19, 20] and be considered to deliver drugs in targeted cancer therapy [21]. In recent years, there has been a great deal of interest in the static and dynamic behavior of pipe (or beam) at micro- and nanoscales [22–27]. Since the size-dependent behavior of microscale structure has been observed experimentally [28–30], we cannot directly extend the analysis of macroscale structures to that of microscale structure. Thus, nonclassical elasticity theories have to been introduced in order to incorporate the size dependence. Recently, a modified couple stress theory capable of capturing the size effect was developed by Yang et al. [31], in which the couple stress tensor is symmetrical and only one internal material length scale parameter is involved. Park and Gao [25] applied the theory to study the static mechanical properties of a Euler–Bernoulli beam. Asghari et al. [27] presented a nonlinear size-dependent Timoshenko beam model based on the modified couple stress theory. In their work, the geometric nonlinearity represented by the von Kármán relationship was adopted to derive the governing equation of motion. For fluid-conveying pipe with microscale, a new theoretical model was developed by Wang [22] for the linear vibration analysis, in which the Euler–Bernoulli model assumption and the modified couple stress theory were employed. In another paper by Xia and Wang [26], the size-dependent vibration analysis of microtube was extended to the Timoshenko model. Yang et al. [23] investigated the microfluid-induced nonlinear free vibration of microtube with both immovable ends by using the modified couple stress theory. In their paper, the geometry nonlinearity arising from the mid-plane stretching was taken into account, and the static postbuckling problem was discussed. Based on the nonlocal strain gradient theory and the Rayleigh beam theory, the wave dispersion in viscoelastic lipid nanotubes conveying protein solution, such as the relationship between wave frequency and wave number and the effect of the wave number and damping coefficient on the stability of tubes, was investigated by Cao and Wang [32]. In contrast to the past literature using the beam theory to model the pipe structures, Wang et al. [33] as well as Li and Wang [34] investigated the vibration behavior of fluid-conveying nanotubes using the shell model.
Perhaps the first study of the dynamics of cantilevered micropipes conveying fluid was contributed by Hosseini and Bahaadini [35], who derived the linear governing equation of motion based on the modified strain gradient theory and then carried out an analysis of eigenvalues with a parametric study in order to examine the effect of the length scale parameter. In another paper by Bahaadini and Hosseini [36], the effect of the fluid slip condition on the free vibration and the flutter instability of viscoelastic cantilevered carbon nanotubes (CNTs) conveying fluid was investigated. The material property of the CNT was simulated by the Kelvin–Voigt viscoelastic constitutive relationship. The equations derived by Hosseini and Bahaadini are linear. Hu et al. [37] developed a nonlinear two-dimensional model for cantilevered fluid-conveying micropipes and explored the possible size-dependent nonlinear responses based on the modified couple stress theory. Guo et al. [38] applied the center manifold theory, normal form method and symmetry to reduce rigorously the equations of motion to a two-degree-of-freedom dynamic system and calculate the corresponding coefficients. Furthermore, the averaging method was employed to investigate the reduced equations for analyzing two types of periodic motions, that is, the planar periodic motions and the spatial periodic motions, along with their stabilities.
However, to the best of our knowledge, no nonlinear vibration model of cantilevered microtube, undergoing large motions and subjected to motion-limiting constraints, is available to date. In this paper, we present a nonlinear model to account for the size effect on the dynamics of cantilevered microtube. The system is modeled as Euler–Bernoulli microtube. The geometry nonlinearity is taken into account by introducing the von Kármán relationship. By the use of modified couple stress theory, the strain energy of the microtube can be calculated. The kinetic expression of the system is available. The Hamilton principle is utilized to derive the governing equation of motion. The application of the Galerkin expansion and modal truncation techniques leads to a set of ordinary differential equations. With the projection method, the Hopf bifurcation is analyzed technically. Finally, the size effect on the dynamics of the system is discussed.
2. Derivation of the Governing Equation of Motion
The system subjected to motion constraints consists of a microtube of length L, cross-sectional area Ap, flexural rigidity EI, and mass per unit length m, conveying an incompressible fluid of mass M per unit length, flowing axially with velocity V not varying with time. The cross section of the microtube is assumed to be symmetric, either circular or rectangular. Using the rectangular Cartesian coordinate system x,y,z, where the x-axis is coincident with the centroidal axis of the undeformed tube, the y-axis is the neutral axis and the z-axis is the symmetric axis as shown in Figure 1.
Schematic of a fluid-conveying microtube with one end free and subjected to motion constraints.
For a slender microtube with a large aspect ratio, the Euler–Bernoulli beam theory can be employed to obtain the displacement components of a point s,y,z on a pipe cross section, with s being the arc length along the centroidal axis, usually written as [25](1)ux=u−zψ,uy=0,uz=w.
In the above equations, u and w represent the displacements of the centerline of the pipe with the coordinate (s, 0, 0) in the x-direction and the z-direction, respectively. ψ is the rotation angle of the cross section with respect to the z-axis. For the case of a cantilevered tube without axial force, it is reasonable to regard the centerline to be inextensible, which leads to [39](2)1+∂u∂s2+∂w∂s2=1.
Based on the Euler–Bernoulli beam theory, the cross sections of the pipe remain perpendicular to the deformed centroidal line. Hence, we find that [38](3)sinψ=∂w∂s.
Before calculating the strain energy, the utilization of the modified couple stress theory will be reviewed first. For more details on this topic, the interested reader may refer to [31, 40]. According to the modified couple theory, the strain energy density is a function of both the strain tensor and the curvature tensor. Furthermore, the strain energy U in a deformed isotropic linear elastic material occupying region Ω can be written as [23, 27](4)U=12∫Ωσijεij+mijχijdv.where the stress tensor σij, the strain tensor εij, the deviatoric part of the couple stress tensor mij, and the symmetric curvature tensor χij are given by(5)σ=λtrεδ+2Gε,(6)ε=12∇u+∇uΤ,(7)m=2l2Gχ,(8)χ=12∇θ+∇θΤ,respectively. In the above equations, λ and G are Lamé’s constants (G is also known as the shear modulus), δ is Kronecker’s delta function, l is a material length scale parameter capturing the size effect on the dynamic behavior, u is the displacement vector with the components given by (1), and θ is the rotation vector given by(9)θ=12curlu.
In the present paper, the geometry nonlinearity will be taken into account, so the axial strain can be expressed by the von Kármán relationship as [27](10)εxx=∂ux∂s+12∂w∂s2=∂u∂s−z∂ψ∂s+12∂w∂s2.
The other nonzero components of the strain tensor can be obtained as(11)εxz=12∂w∂s−ψ.
Then, we have (6) in the form of(12)ε=∂u∂s−z∂ψ∂s+12∂w∂s2012∂w∂s−ψ00012∂w∂s−ψ00.
By combining (12) and (5), one can get(13)σ=∂u∂s−z∂ψ∂s+12∂w∂s2λ+2G0G∂w∂s−ψ0λ∂u∂s−z∂ψ∂s+12∂w∂s20G∂w∂s−ψ0λ∂u∂s−z∂ψ∂s+12∂w∂s2.
The combination of (1) and (9) yields(14)θy=−12ψ+∂w∂s,θx=θz=0.
The substitution of (14) into (8) leads to(15)χ=0−14∂ψ∂s+∂2w∂s20−14∂ψ∂s+∂2w∂s200000.
The utilization of (7) yields(16)m=2l2G0−14∂ψ∂s+∂2w∂s20−14∂ψ∂s+∂2w∂s200000.
By substituting (12), (13), (15), and (16) into equation (4), the expression of the strain energy of the pipe can be obtained as(17)U=12∫Ω∂u∂s−z∂ψ∂s+12∂w∂s22λ+2G+G∂w∂s−ψ2+14l2G∂ψ∂s+∂2w∂s22dv.
Large motions imply that terms of higher order than the linear ones have to be kept in the equation. In the current paper, although the deflection of the pipe can be considered to be large, only cubic-nonlinear terms will be retained in the final equations. Because the variational technique always requires one order higher than the one to be sought, all expressions under the integrand in (17) have to be exact to quartic-nonlinear terms. It is so for the expression of kinetic energy below.
It follows from (2) that(18)u′+12w′2=Ow′4,where the primes denote ∂/∂s. Equation (3) can be written as(19)w′=sinψ=ψ−13!ψ3+Oψ5,which leads to(20)w′−ψ=Oψ3,(21)ψ=w′+13!w′3+Ow′5.
After differentiating (21) with respect to s, we obtain(22)ψ′=w″+12w′2w″+Ow″w′4.
The substitution of (18), (20), and (22) into (17) yields(23)U=12∫Ωz2w″21+w′2λ+2G+14l2Gw″24+2w′2dv=12∫0LIw″21+w′2λ+2G+14Apl2Gw″24+2w′2ds,where I is the second moment of the cross section of the pipe.
By the use of (1), the kinetic energy of the pipe can be expressed as(24)Tp=12∫Ωρ∂ux∂t−z∂ψ∂t2+∂w∂t2dv,where ρ represents the mass per unit volume of the tube.
For a slender pipe, the second moment I without multiplying by E(i.e., the elastic modulus) or G can be neglected. Additionally, the symmetry of the cross section leads to(25)∫Ωρz∂ψ∂t∂ux∂tdv=0.
The kinetic energy of the fluid was given by Semler et al. [39]:(27)Tf=12M∫0L∂u∂t+V∂u∂s+12+∂w∂t+V∂w∂s2ds.
Hence, according to the formulation of Benjamin [1], Hamilton’s principle for a fluid-conveying tube can be written as(28)δ∫t1t2Tp+Tf−U−MV2dt−∫t1t2MVw˙L+VwL′δwLdt=0,where the subscript L represents the values of the corresponding quantities at s=L, and the dots and primes denote ∂/∂t and ∂/∂s, respectively.
After several times of transformations and manipulations, the governing equation of motion is found to be(29)m+Mw¨+2MVw˙′1+w′2+w″MV21+w′2−w″∫sL∫0sm+Mw˙′2+w′w¨′dsds+∫sL2MVw′w˙′+MV2w′w″ds+w′∫0sm+Mw˙′2+w′w¨′ds+EIw4+4w′w″w‴+w″3+w4w′2+Apl2Gw4+2w′w″w‴+12w″3+12w4w′2=0,where the Poisson effect is neglected, and λ+2GI is then replaced by EI.
By assuming that the internal dissipation of the pipe material is viscoelastic and of the Kelvin–Voigt type with coefficient a, according to the analysis given by Semler et al. [39], we only need to add this term aEIw˙4 to (29) to take the dissipation into account.
Furthermore, with the introduction of the impact force modeled as a cubic spring with stiffness K, that is, fw=Kw3 [12], the resulting governing equation of motion can be written as(30)m+Mw¨+2MVw˙′1+w′2+w″MV21+w′2−w″∫sL∫0sm+Mw˙′2+w′w¨′dsds+∫sL2MVw′w˙′+MV2w′w″ds+w′∫0sm+Mw˙′2+w′w¨′ds+EIw4+4w′w″w‴+w″3+w4w′2+Apl2Gw4+2w′w″w‴+12w″3+12w4w′2+aEIw˙4+Kw3δs−sb=0.
By introducing the following nondimensional quantities,(31)η=wL,ξ=sL,τ=EIm+ML41/2t,υ=MEI1/2VL,β=MM+m,κ=KL5EI,ξb=sbL,α=EIm+ML41/2a,one can write (30) in its dimensionless form as follows:(32)η¨+2υβη˙′1+η′2+υ2η″1+η′2−η″∫ξ1∫0ξη˙′2+η′η¨′dξdξ+∫ξ12υβη′η˙′+υ2η′η″dξ+η′∫0ξη˙′2+η′η¨′dξ+η4+4η′η″η‴+η″3+η4η′2+Apl2GEIη4+2η′η″η‴+12η″3+12η4η′2+αη˙4+κη3δξ−ξb=0.
It is observed that there exists two nonlinear inertial terms in (32), which are the same as in [39, 41]. The appearance of nonlinear inertial terms render the equation nonstandard and call the straightforward application of conventional dynamic system theory into question. Because η represents the lateral displacement divided by the length L, η could be assumed to be small. Following the approach used by Li and Paidoussis [41] and Semler et al. [39], the nonlinear inertial terms can be eliminated through a perturbation technique, by which (32) can be reduced to the form without nonlinear inertial terms:(33)η¨+2υβη˙′+υ2η″+1+2l0η4+αη˙4+κη3δξ−ξb+2υβη˙′η′2+υ2η″η′2+3+2l0η′η″η‴+1+l0η″3−l0η4η′2+η′∫0ξη˙′2−2υβη˙″η′−υ2η′η‴+1+2l0η″η4dξ−η″∫ξ12υβη′η˙′+υ2η′η″+1+2l0η″η‴dξ−η″∫ξ1∫0ξη˙′2−2υβη˙″η′−υ2η′η‴+1+2l0η″η4dξdξ=0.
In the above equation, we let l0=1/2Apl2G/EI for simplicity.
3. Discretization of the Governing Equation
By letting(34)ηξ,τ=∑r=1Nφrξqrτ,where(35)φrξ=coshλrξ−cosλrξ−σrsinhλrξ−sinλrξ,σr=sinhλr−sinλrcoshλr+cosλr,r=1,2,...,N,φrξ is the eigenfunction of the cantilevered beam, qrτ is the corresponding generalized coordinate, and λr is the beam eigenvalue. According to Galerkin’s technique, after substituting (34) into (33), multiplication by the factor φrξ followed by integration from 0 to 1 leads to [12](36)q¨i+cijq˙j+kijq˙j+αijklqjqkql+βijklqjqkq˙l+γijklqjq˙kq˙l=0,where cij=δijαλi4+2υβ∫01ϕiϕj′dξ, kij=δij1+2l0∫01ϕi4ϕjdξ+υ2∫01ϕiϕj″dξ, αijkl=∫01ϕi⋅l02ϕj′ϕk″ϕk″′+ϕj″ϕk″ϕl″−ϕj′ϕk′ϕl″′+2ϕj′∫0ξϕk″ϕl4dξ−2ϕj″∫ξ1ϕk″ϕl″′dξ−2ϕj″∫ξ1∫0ξϕk″ϕl4dξdξ+υ2ϕj′ϕk′ϕl″−ϕj′∫0ξϕk′ϕl″′dξ−ϕj″∫ξ1ϕk′ϕl″dξ+ϕj″∫ξ1∫0ξϕk′ϕl″′dξdξ+3ϕj′ϕk″ϕl″′+ϕj″ϕk″ϕl″+ϕj′∫0ξϕk″ϕl4dξ−ϕj″∫ξ1ϕk″ϕl″′dξ−ϕj″∫ξ1∫0ξϕk″ϕl4dξdξdξ, βijkl=∫01ϕi⋅2υβϕj′ϕk′ϕl′−ϕj′∫0ξϕk′ϕl″dξ−ϕj″∫ξ1ϕk′ϕl′dξ+ϕj″∫ξ1∫0ξϕk′ϕl″dξdξdξ, γijkl=∫01ϕi⋅ϕj′∫0ξϕk′ϕl′dξ−ϕj″∫ξ1∫0ξϕk′ϕl′dξdξdξ,where the dots and primes represent ∂/∂τ and ∂/∂ξ, respectively, and the indices i,j,k, and l range from 1 to N.
As shown by Paidoussis and Semler [12], the dynamics in the two-mode version of the analytical model is in good qualitative agreement with the experimental observations. Since the main purpose of this paper is to investigate part of the qualitative behaviors of the present system, the two-mode expansion, that is, N=2 in (34), is adopted to discretize (33).
In order to apply the available tools from dynamic system theory, we render (36) into the first-order form by introducing the following transformations:(37)q1=x1,q2=x2,q˙1=x3,q˙2=x4.
Equation (36) can be rewritten in matrix form as(38)X˙=LX+NX,where(39)L=00100001−k11−k12−c11−c12−k21−k22−c21−c22,NX represents the cubic-nonlinear terms of (36).
4. Theoretical Analysis and Numerical Computations
It should be noted that X=0,0,0,0 is always a solution to (38). The equilibrium is stable for small internal flow velocity and may lose stability when the flow velocity crossed a critical value. As will be shown below, the pipe becomes unstable by flutter and the size effect on the critical velocity is significant. The present work mainly investigates the effect of the internal material length scale parameter on the dynamic behaviors of micropipes and the bifurcation diagrams with increasing flow velocity. Hence, among the five dimensionless quantities υ,β,κ,ξb,α in (32) and the dimensionless quantity l0 in (33), only υ and l0 are considered to be variable parameters. In what follows, we choose the values of the fixed parameters to be(40)α=0.005,β=0.2,ξb=0.82,κ=105,as was also utilized in [12]. υ and l0 are kept as controllable parameters for the determination of critical flow velocity. We study first the loss of stability of the equilibrium, that is, the occurrence of the Hopf bifurcation. Based on the numerical simulation, the graph of υ versus l0, representing the size effect on the critical velocity, is shown in Figure 2. It can be seen that by increasing the value of l0, the critical flow velocity that was predicted from the modified couple stress theory increases.
For the system defined by (38) with α=0.005,β=0.2,ξb=0.82,κ=105,the critical flow velocity versus the nondimensional quantity l0 defined by l0=Apl2G/2EI.
For each point (l˜0,υ˜) of the curve, the Jacobian matrix of (38) at the equilibrium possesses a pair of pure imaginary eigenvalues and two eigenvalues with negative real part, known as a degenerate point. Up to this point, it is potential for the Hopf bifurcation to occur. Technically, we need to verify another two inequalities, that is, με′0≠0 and l10≠0, called the “transversality” and “the first Lyapunov’s coefficient not being zero,” respectively, to complete the analysis, where με′0 and l10 will be defined below.
Physically, for a given microtube, l0 is undoubtedly a constant, but υ remains a controllable parameter. Hence, it is convenient to write (38) near a degenerate point as(41)X˙=Al˜0,υ˜+εX+FX,l˜0,υ˜+ε,where ε represents the perturbation of the critical velocity υ˜.
Letting λε=με±iω0ε be the eigenvalues of matrix Al˜0,υ˜+ε satisfying λ0=±iω0 with ω0>0. qε and pε are chosen as being defined by(42a)Al˜0,υ˜+εqε=λεqε,(42b)AΤl˜0,υ˜+εpε=λ¯εpε,(42c)<pε,qε>=1,where <⋅,⋅> represents the scalar product defined by <p,q>=p¯Τq.
By differentiating (42a) with respect to ε, we find that(43)Aε′l˜0,υ˜q0+Al˜0,υ˜qε′0=λε′0q0+λ0qε′0.
By multiplying by p¯0Τ, one obtains(44)λε′0=<p0,Aε′l˜0,υ˜q0>.
Thus(45)με′0=Re<p0,Aε′l˜0,υ˜q0>.
The calculation of l10 is related to two aspects: one associated with the center manifold theory and the other associated with the normal form method. In what follows, we only give the formulation of l10. For more details on this topic, the interested reader is referred to [42, 43]. For this purpose, we write (41) with ε=0 as(46)X˙=Al˜0,υ˜X+FX,l˜0,υ˜=Al˜0,υ˜X+12BX,X,l˜0,υ˜+13!CX,X,X,l˜0,υ˜,where BX,X,l˜0,υ˜ and CX,X,X,l˜0,υ˜ are the symmetric multilinear vector functions given by(47)BiX,X,l˜0,υ˜=∑j,k=14∂2FiY,l˜0,υ˜∂yj∂ykY=0xjxk,i=1,2,3,4,(48)CiX,X,X,l˜0,υ˜=∑j,k,l=14∂3FiY,l˜0,υ˜∂yj∂yk∂ylY=0xjxkxl,i=1,2,3,4,respectively.
For the symmetry of the current system, we have BX,X,l˜0,υ˜=0.
The utilization of the projection method [42] leads to(49)l10=12ω0Re<p0,Cq0,q0,q¯0,l˜0,υ˜>.
With the aid of a computer, the plots of μ′0 and l10 that are associated with each value of l˜0,υ˜ on the curve in Figure 2 can be obtained, as shown in Figures 3 and 4.
For the system of Figure 2, the curve of με′0 associated with each point (l˜0,υ˜).
For the system of Figure 2, the values of the first Lyapunov’s coefficient l10 at each point (l˜0,υ˜).
As can be seen from the two figures above, μ′0 and l0 are not equal to zero and then the Hopf bifurcation can take place. The stability of this bifurcated solution depends on the sign of l10. The periodic solution is stable in case l10 is negative while it is unstable when l10 is positive, as said to be “supercritical” and “subcritical” bifurcation, respectively. A variant of the macrotube system that was modified to contain a “friction factor” was studied by Bajaj et al. [11], and the subcritical Hopf bifurcation might occur for certain regions of the parameter space. For the present system, it should be noted that l10 is always negative and then the Hopf bifurcation is supercritical. In what follows, the dynamics will be explored numerically.
For various values of l0, that is, at 0, 0.2, 0.5, and 1, the calculations have produced the bifurcation diagrams shown in Figures 5–8. In these figures the displacement plotted in the ordinate is the amplitude x1 and the variable parameter is the dimensionless flow velocity υ. Having discarded the transient solutions, whenever the velocity x4 was zero, we recorded the displacement x1, constructing the bifurcation diagrams with positive and negative values, as shown in Figures 5–8. It should be mentioned that the model will reduce to the classical tube model when the size effect is suppressed by letting l0=0.
The bifurcation diagrams of the nonlinear system defined by (38) for α=0.005,β=0.2,ξb=0.82,κ=105, and l0=0, as υ is varied. (a) For 5.4≤υ≤8 and (b), (c), (d), and (e) for smaller ranges of υ to show finer details.
The bifurcation diagrams of the nonlinear system defined by (38) for α=0.005,β=0.2,ξb=0.82,κ=105, and l0=0.2, as υ is varied. (a) For 6.4≤υ≤9.5 and (b), (c), (d), and (e) for smaller ranges of υ, showing more details.
The bifurcation diagram for the nonlinear system of Figure 6 with l0=0.5, as υ is varied.
The bifurcation diagram for the nonlinear system of Figure 7 with l0=1, as υ is varied.
In these figures, it is observed that the qualitative dynamics of the system is similar to each other. The Hopf bifurcation occurs first and then the pipe undergoes limit cycle motion with symmetric trajectory, which is followed by a pitchfork bifurcation thereafter the limit cycle lost its symmetry. As υ increases further, a sequence of period-doubling bifurcation arises, leading the motions to chaos. Beyond the chaotic band, the system reverts to the periodic motion of period 3. From the point of view of singularity theory, in spite of the different physical assumptions made for the system, when neglecting the terms due to gravity and introducing the material length scale parameter l, the pipe behaviors are qualitatively like those of Paidoussis et al. [8] as well as Paidoussis and Semler [12]. In fact, the similar bifurcation sequence can be also found in other vibroimpact system with symmetry [44]. However, for the reliable and optical design of microtube, it is vital for us to distinguish critical flow velocity υ, the threshold of chaotic motion, and so forth, for different values of l0. We then compare these bifurcation diagrams in order to explore the effect of the parameter l0 on the system under consideration.
For brevity, we only choose two figures, corresponding to l0=0 and l0=0.2, respectively, to complete this analysis. For this purpose, it is convenient to amplify these two bifurcation diagrams over a smaller range of υ as shown in Figures 5(b)–5(e) and Figures 6(b)–6(e). Thus, the finer details can be observed. The Hopf bifurcation is noted at υ=5.6784138 for l0=0, while it was noted at υ=6.6766354, a higher flow velocity, for l0=0.2. A series of υ, that is, 6.905075, 7.245702, 7.303842, 7.316114, and so forth, at which the pitchfork bifurcation and a sequence of period-doubling bifurcations occur for the case of l0=0, is lower than that for l0=0.2, that is, 8.141000, 8.545417, 8.614248, 8.628766, and so forth. It is so for the occurrence of chaos and the resumption of periodic motions. Based on the results, it seems that the internal material length scale parameter makes the microtube more stable.
5. Conclusions
In this paper, the full nonlinear equation of a fluid-conveying cantilevered microtube undergoing large motions and subjected to motion-limiting constraints was deduced based on the Hamilton principle and the modified couple stress theory. The governing equation of motion was discretized into a set of ordinary differential equations via Galerkin’s method. The effect of the internal material length scale parameter on the critical flow velocity was investigated. By using the projection method, the supercritical Hopf bifurcation was demonstrated. The results show that the size effect on the vibration properties is significant. The following outcomes were found from the results:
The size effect on the governing equation, especially on the nonlinear terms in the governing equation was obtained.
The results showed that, by increasing the value of the nondimensional internal material length scale parameter l0, the critical flow velocity υ predicted from the modified couple stress theory increases.
At each degenerate point l˜0,υ˜ of the curve of the critical flow velocity versus the nondimensional quantity l0, the real part of the derivative of the eigenvalues with respect to the flow velocity με′0 is positive, and the first Lyapunov’s coefficient l10 is negative. Therefore, the existence of flutter (i.e., the supercritical Hopf bifurcation) was demonstrated rigorously.
By comparing the bifurcation diagrams with different l0, it is found that the flow velocities at which various bifurcations and chaos occur, increased with increasing the value of the nondimensional internal material length scale parameter l0; that is, the internal material length scale parameter makes the microtube more stable.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The author declares no conflicts of interest.
Acknowledgments
This work was supported by Guizhou Institute of Technology under Grant XJGC20190909.
BenjaminT. B.Dynamics of a system of articulated pipes conveying fluid-I. Theory1961261130745748610.1098/rspa.1961.0090GregoryR. W.PaidoussisM. P.Unstable oscillation of tubular cantilevers conveying fluid I. Theory1966293143551252710.1098/rspa.1966.0187PaïdoussisM. P.IssidN. T.Dynamic stability of pipes conveying fluid197433326729410.1016/s0022-460x(74)80002-72-s2.0-0000163696HolmesP. J.Bifurcations to divergence and flutter in flow-induced oscillations: a finite dimensional analysis197753447150310.1016/0022-460x(77)90521-12-s2.0-0017522897RousseletJ.HerrmannG.Dynamic behavior of continuous cantilevered pipes conveying fluid near critical velocities198148494394710.1115/1.31577602-s2.0-33745530376NamchchivayaN. S.Non-linear dynamics of supported pipe conveying pulsating fluid-I. Subharmonic resonance198924185196NamchchivayaN. S.TienW. M.Non-linear dynamics of supported pipe conveying pulsating fluid-II. Combination resonance198924319720810.1016/0020-7462(89)90038-32-s2.0-0024928193PaïdoussisM. P.LiG. X.MoonF. C.Chaotic oscillations of the autonomous system of a constrained pipe conveying fluid1989135111910.1016/0022-460x(89)90750-52-s2.0-0024767278WangL.A further study on the non-linear dynamics of simply supported pipes conveying pulsating fluid200944111512110.1016/j.ijnonlinmec.2008.08.0102-s2.0-56949108308HolmesP. J.Pipes supported at both ends cannot flutter197845361962210.1115/1.34243712-s2.0-0018018406BajajA. K.SethnaP. R.LundgrenT. S.Hopf bifurcation phenomena in tubes carrying a fluid198039221323010.1137/01390192-s2.0-0019072146PaidoussisM. P.SemlerC.Nonlinear and chaotic oscillations of a constrained cantilevered pipe conveying fluid: a full nonlinear analysis19944665567010.1007/bf001622362-s2.0-0027850767JinJ.-D.Stability and chaotic motions of a restrained pipe conveying fluid1997208342743910.1006/jsvi.1997.11952-s2.0-0031553192WangL.NiQ.A note on the stability and chaotic motions of a restrained pipe conveying fluid20062964-51079108310.1016/j.jsv.2006.03.0162-s2.0-33644786858LiangF.WenB.Forced vibrations with internal resonance of a pipe conveying fluid under external periodic excitation201124647748310.1016/s0894-9166(11)60047-52-s2.0-84855332988LiB.WangZ.JingL.Dynamic response of pipe conveying fluid with lateral moving supports2018201817329578710.1155/2018/32957872-s2.0-85057022362ZhangY.-L.FengH.-R.ChenL.-Q.Supercritical nonlinear vibration of a fluid-conveying pipe subjected to a strong external excitation2016201621390749810.1155/2016/39074982-s2.0-84955158770RinaldiS.PrabhakarS.VengallatoreS.PaïdoussisM. P.Dynamics of microscale pipes containing internal fluid flow: damping, frequency shift, and stability201032981081108810.1016/j.jsv.2009.10.0252-s2.0-72049096868NajmzadehM.HaaslS.EnokssonP.A silicon straight tube fluid density sensor20071781657166310.1088/0960-1317/17/8/0322-s2.0-34547683510SparksD.SmithR.CruzV.Dynamic and kinematic viscosity measurements with a resonating microtube20091491384110.1016/j.sna.2008.09.0132-s2.0-58149191366BhirdeA. A.PatelV.GavardJ.Targeted killing of cancer cells in vivo and in vitro with EGF-directed carbon nanotube-based drug delivery20093230731610.1021/nn800551s2-s2.0-63449106550WangL.Size-dependent vibration characteristics of fluid-conveying microtubes201026467568410.1016/j.jfluidstructs.2010.02.0052-s2.0-77952668476YangT.-Z.JiS.YangX.-D.FangB.Microfluid-induced nonlinear free vibration of microtubes201476475510.1016/j.ijengsci.2013.11.0142-s2.0-84891720570TangM.NiQ.WangL.LuoY.WangY.Nonlinear modeling and size-dependent vibration analysis of curved microtubes conveying fluid based on modified couple stress theory20148411010.1016/j.ijengsci.2014.06.0072-s2.0-84904165150ParkS. K.GaoX.-L.Bernoulli-Euler beam model based on a modified couple stress theory200616112355235910.1088/0960-1317/16/11/0152-s2.0-33750597812XiaW.WangL.Microfluid-induced vibration and stability of structures modeled as microscale pipes conveying fluid based on non-classical Timoshenko beam theory201094-595596210.1007/s10404-010-0618-z2-s2.0-77957557954AsghariM.KahrobaiyanM. H.AhmadianM. T.A nonlinear Timoshenko beam formulation based on the modified couple stress theory201048121749176110.1016/j.ijengsci.2010.09.0252-s2.0-78649635287FleckN. A.MullerG. M.AshbyM. F.HutchinsonJ. W.Strain gradient plasticity: theory and experiment199442247548710.1016/0956-7151(94)90502-92-s2.0-0028370983LamD. C. C.YangF.ChongA. C. M.WangJ.TongP.Experiments and theory in strain gradient elasticity20035181477150810.1016/s0022-5096(03)00053-x2-s2.0-0038721068McFarlandA. W.ColtonJ. S.Role of material microstructure in plate stiffness with relevance to microcantilever sensors20051551060106710.1088/0960-1317/15/5/0242-s2.0-24144462483YangF.ChongA. C. M.LamD. C. C.TongP.Couple stress based strain gradient theory for elasticity200239102731274310.1016/s0020-7683(02)00152-x2-s2.0-0037150391CaoD. Y.WangY. Q.Wave dispersion in viscoelastic lipid nanotubes conveying viscous protein solution202013512410.1140/epjp/s13360-019-00074-3WangY. Q.LiH. H.ZhangY. F.ZuJ. W.A nonlinear surface-stress-dependent model for vibration analysis of cylindrical nanoscale shells conveying fluid201864557010.1016/j.apm.2018.07.0162-s2.0-85050677186LiZ. H.WangY. Q.Vibration and stability analysis of lipid nanotubes conveying fluid201923111210.1007/s10404-019-2290-22-s2.0-85073573511HosseiniM.BahaadiniR.Size dependent stability analysis of cantilever micro-pipes conveying fluid based on modified strain gradient theory201610111310.1016/j.ijengsci.2015.12.0122-s2.0-84955491274BahaadiniR.HosseiniM.Effects of nonlocal elasticity and slip condition on vibration and stability analysis of viscoelastic cantilever carbon nanotubes conveying fluid201611415115910.1016/j.commatsci.2015.12.0272-s2.0-84952888252HuK.WangY. K.DaiH. L.WangL.QianQ.Nonlinear and chaotic vibrations of cantilevered micropipes conveying fluid based on modified couple stress theory20161059310710.1016/j.ijengsci.2016.04.0142-s2.0-84970997360GuoY.XieJ.WangL.Three-dimensional vibration of cantilevered fluid-conveying micropipes-types of periodic motions and small-scale effect201810211213510.1016/j.ijnonlinmec.2018.04.0012-s2.0-85045394648SemlerC.LiG. X.Paı̈doussisM. P.The non-linear equations of motion of pipes conveying fluid1994169557759910.1006/jsvi.1994.10352-s2.0-0028762577KongS.ZhouS.NieZ.WangK.The size-dependent natural frequency of Bernoulli-Euler micro-beams200846542743710.1016/j.ijengsci.2007.10.0022-s2.0-41749123679LiG. X.PaidoussisM. P.Stability, double degeneracy and chaos in cantilevered pipes conveying fluid19942918310710.1016/0020-7462(94)90054-x2-s2.0-0028336585KuznetsovY. A.1998SecondBerlin, GermanySpringer-VerlagIoossG.JosephD. D.1990SecondBerlin, GermanySpringer-VerlagYueY.XieJ.-H.Symmetry, cusp bifurcation and chaos of an impact oscillator between two rigid sides20072881109111710.1007/s10483-007-0813-z2-s2.0-38149126391