Flexural–torsional buckling of misaligned axially moving beams. I. Three-dimensional modeling,equilibria, and bifurcations

The equations of motion for the flexural–flexural–torsional–extensional dynamics of a beam are generalized to the field of axially moving continua by including the effects of translation speed and initial tension. The governing equations are simplified on the basis of physically justifiable assumptions and are shown to reduce to simpler models published in the literature. The resulting nonlinear equations of motion are used to investigate the flexural–torsional buckling of translating continua such as belts and tapes caused by parallel pulley misalignment.The effect of pulley misalignment on the steady motion (equilibrium) solutions and the bifurcation characteristics of the system are investigated numerically. The system undergoes multiple pitchfork bifurcations as misalignment is increased, with out-of-plane equilibria born at each bifurcation. The amount of misalignment to cause buckling and the post-buckled shapes are determined for various translation speeds and ratios of the flexural stiffnesses in the two bending planes. Increasing translation speed decreases the misalignment necessary to cause flexural–torsional buckling. In Part II of the present work, the stability and vibration characteristics of the planar and non-planar equilibria are analyzed.     

Hopf分岔的代数判据及其在车辆动力学中的应用

利用Hurwitz行列式,给出平衡点失稳而发生Hopf分岔的代数判定准则和计算方法,这一方法将Hopf分岔点的求解转化为一个非线性方程的求解问题,从而克服了以前方法在计算Hopf分岔点时,对于参数的每一次变化通过求特征根并判定特征根的实部是否为零的庞大工作量。应用这一方法,我们进行了非线性车辆系统蛇行运动稳定性的研究,得到了轮对系统发生蛇行运动的临界速度的解析表达式。     

高超声速飞行器横侧向失稳非线性分岔分析

针对滑翔式高超声速飞行器大攻角横侧向失稳问题,采用延拓算法和分岔理论,求解并分析了以俯仰舵偏为连续参数的稳态平衡分岔图和以副翼舵偏为连续参数的横侧向机动稳态平衡分岔图,对平衡分支的稳定性和突变点进行了分析,并给出了特征根拓扑结构变化.研究表明,高超声速飞行器存在极限分岔点、Hopf分岔点以及叉型分岔点,且从叉型分岔点延伸出多个平衡分支,引起横侧向的自滚转失稳;从Hopf分岔点延伸出极限环分支,该分支对应较为复杂的极限环运动,其中还包含倍周期分岔、花环分岔、极限环极限点分岔等复杂的分岔现象;在横侧向机动飞行情况下,模型存在横向操作偏离失稳问题,且存在多个不稳定的平衡点.研究结果为实现高超声速飞行器的稳定飞行和控制器的设计提供了极其重要的动力学信息.      

Aeroelastic analysis of helicopter rotor blade in hover using an efficient reduced-order aerodynamic model

This paper presents a coupled flap–lag–torsion aeroelastic stability analysis and response of a hingeless helicopter blade in the hovering flight condition. The boundary element method based on the wake eigenvalues is used for the prediction of unsteady airloads of the rotor blade. The aeroelastic equations of motion of the rotor blade are derived by Galerkin's method. To obtain the aeroelastic stability and response, the governing nonlinear equations of motion are linearized about the nonlinear steady equilibrium positions using small perturbation theory. The equilibrium deflections are calculated through the iterative Newton–Raphson method. Numerical results comprising steady equilibrium state deflections, aeroelastic eigenvalues and time history response about these states for a two-bladed rotor are presented, and some of them are compared with those obtained from a two-dimensional quasi-steady strip aerodynamic theory. Also, the effect of the number of aerodynamic eigenmodes is investigated. The results show that the three-dimensional aerodynamic formulation has considerable impact on the determination of both the equilibrium condition and lead-lag instability.     

Pseudo-nonlinear dynamic analysis of buckled pipes

In this study, the post-divergence behavior of fluid-conveying pipes supported at both ends is investigated using the nonlinear equations of motion. The governing equation exhibits a cubic nonlinearity arising from mid-plane stretching. Exact solutions for post-buckling configurations of pipes with fixed–fixed, fixed–hinged, and hinged–hinged boundary conditions are investigated. The pipe is stable at its original static equilibrium position until the flow velocity becomes high enough to cause a supercritical pitchfork bifurcation, and the pipe loses stability by static divergence. In the supercritical fluid velocity regime, the equilibrium configuration becomes unstable and bifurcates into multiple equilibrium positions. To investigate the vibrations that occur in the vicinity of a buckled equilibrium position, the pseudo-nonlinear vibration problem around the first buckled configuration is solved precisely using a new solution procedure. By solving the resulting eigenvalue problem, the natural frequencies and the associated mode shapes of the pipe are calculated. The dynamic stability of the post-buckling configurations obtained in this manner is investigated. The first buckled shape is a stable equilibrium position for all boundary conditions. The buckled configurations beyond the first buckling mode are unstable equilibrium positions. The natural frequencies of the lowest vibration modes around each of the first two buckled configurations are presented. Effects of the system parameters on pipe behavior as well as the possibility of a subcritical pitchfork bifurcation are also investigated. The results show that many internal resonances might be activated among the vibration modes around the same or different buckled configurations.     

Using Newton-GMRES for viscoelastic flow time-steppers

Kinetic theory models exhibit dynamics that depend on a few low-order moments of the underlying conformational distribution function. This dependence is exhibited in a compact spectrum of eigenvalues for the Jacobian matrix associated with the dynamical system. We take advantage of this spectrum of eigenvalues through Newton-GMRES iterations to enable dynamic viscoelastic simulators (time-steppers) to obtain stationary states and perform stability/bifurcation analysis. Results are presented for three example problems: (1) the equilibrium behavior of the Doi model with the Onsager excluded volume potential, (2) pressure-driven flow of non-interacting rigid dumbbells in a planar channel, and (3) pressure-driven flow of non-interacting rigid dumbbells through a planar channel with a linear array of equally spaced cylinders.     

Stability and bifurcation for limit cycle oscillations of an airfoil with external store

In this paper, stability and local bifurcation behaviors for the nonlinear aeroelastic model of an airfoil with external store are investigated using both analytical and numerical methods. Three kinds of degenerated equilibrium points of bifurcation response equations are considered. They are characterized as (1) one pair of purely imaginary eigenvalues and two pairs of conjugate complex roots with negative real parts; (2) two pairs of purely imaginary eigenvalues in nonresonant case and one pair of conjugate complex roots with negative real parts; (3) three pairs of purely imaginary eigenvalues in nonresonant case. With the aid of Maple software and normal form theory, the stability regions of the initial equilibrium point and the explicit expressions of the critical bifurcation curves are obtained, which can lead to static bifurcation and Hopf bifurcation. Under certain conditions, 2-D tori motion may occur. The complex dynamical motions are considered in this paper. Finally, the numerical solutions achieved by the fourth-order Runge–Kutta method agree with the analytic results.     

Stability and bifurcation for a flexible beam under a large linear motion with a combination parametric resonance

Stability and bifurcation behaviors for a model of a flexible beam undergoing a large linear motion with a combination parametric resonance are studied by means of a combination of analytical and numerical methods. Three types of critical points for the bifurcation equations near the combination resonance in the presence of internal resonance are considered, which are characterized by a double zero and two negative eigenvalues, a double zero and a pair of purely imaginary eigenvalues, and two pairs of purely imaginary eigenvalues in nonresonant case, respectively. The stability regions of the initial equilibrium solution and the critical bifurcation curves are obtained in terms of the system parameters. Especially, for the third case, the explicit expressions of the critical bifurcation curves leading to incipient and secondary bifurcations are obtained with the aid of normal form theory. Bifurcations leading to Hopf bifurcations and 2-D tori and their stability conditions are also investigated. Some new dynamical behaviors are presented for this system. A time integration scheme is used to find the numerical solutions for these bifurcation cases, and numerical results agree with the analytic ones.     

On the effect of nonsmooth Coulomb friction on Hopf bifurcations in a 1-DoF oscillator with self-excitation due to negative damping

This article deals with the question, to what extent damping due to nonsmooth Coulomb friction may affect the stability and bifurcation behavior of vibrational systems with self-excitation due to negative effective damping which??for the smooth case??is related to a Hopf bifurcation of the steady state. Without damping due to Coulomb friction, the stability of the trivial solution is controlled by the effective viscous damping of the system: as the damping becomes negative, the steady state loses stability at a Hopf point. Adding Coulomb friction changes the trivial solution into a set of equilibria, which??for oscillatory systems??is asymptotically stable for all values of effective viscous damping. The Hopf point vanishes and an unstable limit cycle appears which borders the basin of attraction of the equilibrium set. Moreover, the influence of nonlinear damping terms is discussed. The effect of Coulomb frictional damping may be seen as adding an imperfection to the classical smooth Hopf scenario: as the imperfection vanishes, the behavior of the smooth problem is recovered.     

挠性联结双体航天器的稳定性与分岔

研究圆轨道内受万有引力矩作用的挠性联结双体航天器在轨道平面内的姿态运动,讨论其相对轨道坐标系统平衡状态的稳定性与分岔。提出判平衡方程非平凡解存在性的几何方法,并应用Ｌｉａｐｕｎｏｖ直接法、Ｌｉａｐｕｎｏｖ－Ｓｃｈｍｉｄｔ约化方法和奇异性理论导出解析形式的稳定性与分岔的充要条件,从而对系统的全局运动性态作出定性的描述。     

Supercritical equilibrium solutions of axially moving beams with hybrid boundary conditions

In this paper supercritical equilibria and critical speeds of axially moving beams constrained by sleeves with torsion springs are deduced. Transverse vibration of the beams is governed by a nonlinear integro-partial-differential equation. In the supercritical regime, the corresponding static equilibrium equation for the hybrid boundary conditions is analytically solved for the equilibria and the critical speeds. In the view of the non-trivial equilibrium, comparisons are made among the integro-partial-differential equation, a nonlinear partial-differential equation for transverse vibration, and coupled equations for planar motion under the hybrid boundary conditions.     

A delayed ratio-dependent predator–prey model of interacting populations with Holling type III functional response

In this article, we study a ratio-dependent predator–prey model described by a Holling type III functional response with time delay incorporated into the resource limitation of the prey logistic equation. This investigation includes the influence of intra-species competition among the predator species. All the equilibria are characterized. Qualitative behavior of the complicated singular point (0,0) in the interior of the first quadrant is investigated by means of a blow-up transformation. Uniform persistence, stability, and Hopf bifurcation at the positive equilibrium point of the system are examined. Global asymptotic stability analyses of the positive equilibrium point by the Bendixon–Dulac criterion for non-delayed model and by constructing a suitable Lyapunov functional for the delayed model are carried out separately. We perform a numerical simulation to validate the applicability of the proposed mathematical model and our analytical findings.     

A timestepper approach for the systematic bifurcation and stability analysis of polymer extrusion dynamics

We discuss how matrix-free/timestepper algorithms can efficiently be used with dynamic non-Newtonian fluid mechanics simulators in performing systematic stability/bifurcation analysis. The timestepper approach to bifurcation analysis of large-scale systems is applied to the plane Poiseuille flow of an Oldroyd-B fluid with non-monotonic slip at the wall, in order to further investigate a mechanism of extrusion instability based on the combination of viscoelasticity and non-monotonic slip. Due to the non-monotonicity of the slip equation the resulting steady-state flow curve is non-monotonic and unstable steady states appear in the negative-slope regime. It has been known that self-sustained oscillations of the pressure gradient are obtained when an unstable steady state is perturbed [M.M. Fyrillas, G.C. Georgiou, D. Vlassopoulos, S.G. Hatzikiriakos, A mechanism for extrusion instabilities in polymer melts, Polymer Eng. Sci. 39 (1999) 2498–2504].Treating the simulator of a distributed parameter model describing the dynamics of the above flow as an input–output “black-box” timestepper of the state variables, stable and unstable branches of both equilibrium and periodic oscillating solutions are computed and their stability is examined. It is shown for the first time how equilibrium solutions lose stability to oscillating ones through a subcritical Hopf bifurcation point which generates a branch of unstable limit cycles and how the stable periodic solutions lose their stability through a critical point which marks the onset of the unstable limit cycles. This implicates the coexistence of stable equilibria with stable and unstable periodic solutions in a narrow range of volumetric flow rates.     

Bifurcations and chaotic dynamics in suspended cables under simultaneous parametric and external excitations

The bifurcations and chaotic dynamics of parametrically and externally excited suspended cables are investigated in this paper. The equations of motion governing such systems contain quadratic and cubic nonlinearities, which may result in two-to-one and one-to-one internal resonances. The Galerkin procedure is introduced to simplify the governing equations of motion to ordinary differential equations with two-degree-of-freedom. The case of one-to-one internal resonance between the modes of suspended cables, primary resonant excitation, and principal parametric excitation of suspended cables is considered. Using the method of multiple scales, a parametrically and externally excited system is transformed to the averaged equations. A pseudo arclength scheme is used to trace the branches of the equilibrium solutions and an investigation of the eigenvalues of the Jacobian matrix is used to assess their stability. The equilibrium solutions experience pitchfork, saddle-node, and Hopf bifurcations. A detailed bifurcation analysis of the dynamic (periodic and chaotic) solutions of the averaged equations is presented. Five branches of dynamic solutions are found. Three of these branches that emerge from two Hopf bifurcations and the other two are isolated. The two Hopf bifurcation points, one is supercritical Hopf bifurcation point and another is primary Hopf bifurcation point. The limit cycles undergo symmetry-breaking, cyclic-fold, and period-doubling bifurcations, whereas the chaotic attractors undergo attractor-merging, boundary crises. Simultaneous occurrence of the limit cycle and chaotic attractors, homoclinic orbits, homoclinic explosions and hyperchaos are also observed.     

Stability and Hamiltonian Hopf bifurcation for a nonlinear symmetric bladed rotor

The modal interaction which leads to Hamiltonian Hopf bifurcation is studied for a nonlinear rotating bladed-disk system. The model, which is discussed in the paper, is a Jeffcott rotor carrying a number of planar blades which bend in the plane of the motion. The rigid rotating disk is supported on nonlinear bearings. It is supposed that this dynamical system is a Hamiltonian system which is perturbed by small dissipative and nonlinear forces. Krein’s theorem is employed for obtaining a stability criterion. The nonlinear eigenvalue equations on the stability boundary are turned into ordinary differential equations (ODEs) by differentiating them over the rotating speed. By solving these ODEs, the eigenmodes and the eigenvalues on the stability boundary are obtained. The bifurcation analysis is performed by applying multiple scales method around the boundary. The rotor nonlinear behavior and damping effects are studied for different conditions on the rotating speed and nonlinearity type by the bifurcation equation. It is shown that the damping distribution between the blades and bearings may shift the unstable mode. Depending on the nonlinearity type, subcritical and supercritical Hopf bifurcation are possible.     

Bifurcation Equations Through Multiple-Scales Analysis for a Continuous Model of a Planar Beam

The Multiple-Scale Method is applied directly to a one-dimensional continuous model to derive the equations governing the asymptotic dynamic of the system around a bifurcation point. The theory is illustrated with reference to a specific example, namely an internally constrained planar beam, equipped with a lumped viscoelastic device and loaded by a follower force. Nonlinear, integro-differential equations of motion are derived and expanded up to cubic terms in the transversal displacements and velocities of the beam. They are put in an operator form incorporating the mechanical boundary conditions, which account for the lumped viscoelastic device; the problem is thus governed by mixed algebraic-integro-differential operators. The linear stability of the trivial equilibrium is first studied. It reveals the existence of divergence, Hopf and double-zero bifurcations. The spectral properties of the linear operator and its adjoint are studied at the bifurcation points by obtaining closed-form expressions. Notably, the system is defective at the double-zero point, thus entailing the need to find a generalized eigenvector. A multiple-scale analysis is then performed for the three bifurcations and the relevant bifurcation equations are derived directly in their normal forms. Preliminary numerical results are illustrated for the double-zero bifurcation.     

Equilibria of axially moving beams in the supercritical regime

Equilibria of axially moving beams are computationally investigated in the supercritical transport speed ranges. In the supercritical regime, the pattern of equilibria consists of the straight configuration and of non-trivial solutions that bifurcate with transport speed. The governing equations of coupled planar is reduced to a partial-differential equation and an integro-partial-differential equation of transverse vibration. The numerical schemes are respectively presented for the governing equations and the corresponding static equilibrium equation of coupled planar and the two governing equations of transverse motion for non-trivial equilibrium solutions via the finite difference method and differential quadrature method under the simple support boundary. A steel beam is treated as example to demonstrate the non-trivial equilibrium solutions of three nonlinear equations. Numerical results indicate that the three models predict qualitatively the same tendencies of the equilibrium with the changing parameters and the integro-partial-differential equation yields results quantitatively closer to those of the coupled equations.     

On the equilibrium configurations of an elastically constrained rotating disk: An analytical approach

According to the linear theory of vibration for spinning disks, the backward traveling waves of some of the modes may have zero natural frequency at what are called the critical speeds. At these speeds, the linear equations of motion cannot properly predict the amplitude response of the spinning disk, and nonlinear equations of motion must be used. In this paper, geometrical nonlinear equations of motion based on Von Karman plate theory are employed to study the dynamics of an elastically constrained disk near its critical speeds. A one-mode approximation is used to examine the effect of elastic constraint on the amplitude response. Presenting the equations in a space-fixed coordinate system, this study aims to find closed-form solutions for some of the equilibrium configurations of an elastically constrained spinning disk. Also, the stability of these configurations is studied using analytical techniques. It is shown that below the critical speed, one neutrally stable equilibrium solution exists, while above it, a bifurcation occurs. In this situation, two more branches of equilibrium configurations emerge, one of which is neutrally stable and the other unstable. Closed-form expressions for the bifurcation points are obtained. Due to the effect of an elastic constraint, a bifurcation occurs and the previously neutrally stable equilibrium configuration turns unstable. Also at this bifurcation point, two more branches of equilibrium solutions emerge.     

CO-DIMENSION 2 BIFURCATIONS AND CHAOS IN CANTILEVERED PIPE CONVEYING TIME VARYING FLUID WITH THREE-TO-ONE INTERNAL RESONANCES

The nonlinear behavior of a cantilevered fluid conveying pipe subjected to principal parametric and internal resonances is investigated in this paper. The flow velocity is divided into constant and sinusoidai parts. The velocity value of the constant part is so adjusted such that the system exhibits 3:1 internal resonances for the first two modes. The method of multiple scales is employed to obtain the response of the system and a set of four first-order nonlinear ordinary-differential equations for governing the amplitude of the response. The eigenvalues of the Jacobian matrix are used to assess the stability of the equilibrium solutions with varying parameters. The codimension 2 derived from the double-zero eigenvaiues is analyzed in detail. The results show that the response amplitude may undergo saddle-node, pitchfork, Hopf, homoclinic loop and period-doubling bifurcations depending on the frequency and amplitude of the sinusoidal flow. When the frequency of the sinusoidal flow equals exactly half of the first-mode frequency, the system has a route to chaos by period-doubling bifurcation and then returns to a periodic motion as the amplitude of the sinusoidal flow increases.     

Bifurcation of a modified railway wheelset model with nonlinear equivalent conicity and wheel–rail force

This paper presents the bifurcation behaviors of a modified railway wheelset model to explore its instability mechanisms of hunting motion. Equivalent conicity data measured from China high-speed railway vehicle are used to modify the wheelset model. Firstly, the relationships between longitudinal stiffness, lateral stiffness, equivalent conicity and critical speed are taken into account by calculating the real parts of the eigenvalues of the Jacobian matrix and Hurwitz criterion for the corresponding linear model. Secondly, measured equivalent conicity data are fitted by a nonlinear function of the lateral displacement rather than are considered as a constant as usual. Nonlinear wheel–rail force function is used to describe the wheel–rail contact force. Based on these modifications, a modified railway wheelset model with nonlinear equivalent conicity and wheel–rail force is set up, and then, some instability mechanisms of China high-speed train vehicle are investigated based on Hopf bifurcation, fold (limit point) bifurcation of cycles, cusp bifurcation of cycles, Neimark–Sacker bifurcation of cycles and 1:1 resonance. In particular, fold bifurcation of cycles can produce a vast effect on the hunting motion of the modified wheelset model. One of the main reasons leading to hunting motion is due to the fold bifurcation structure of cycles, in which stable limit cycles and unstable limit cycles may coincide, and multiple nested limit cycles appear on a side of fold bifurcation curve of cycles. Unstable hunting motion mainly depends on the coexistence of equilibria and limit cycles and their positions; if the most outward limit cycle is stable, then the motion of high-speed vehicle should be safe in a reasonable range. Otherwise, if the initial values are chosen near the most outward unstable limit cycle or the system is perturbed by noises, the high-speed vehicle will take place unstable hunting motion and even lead to serious train derailment events. Therefore, in order to control hunting motions, it may be the easiest way in theory to guarantee the coexistence of the inner stable equilibrium and the most outward stable limit cycle in a wheelset system.