Perturbation-based prediction of vibration phase shift along fluid-conveying pipes due to Coriolis forces,nonuniformity, and nonlinearity

This paper investigates the influence of tread conicity variation on hunting dynamical changing features for railway vehicle. Nonlinear contact force between wheels and rail is estimated by Vermeulen–Johnson creep force law. And a piecewise linear function is employed to appropriate the collision between wheel flange and rail. Hunting asymmetrical motions are analyzed by lateral bifurcation behaviors between maximum and minimum of car body lateral displacement. The result shows that the critical speed decreases with the increase in tread conicity, while the speed gap between linear and nonlinear speeds is narrowing. Compared with wheelsets, lateral amplitudes of the bogies are vulnerable to the tread conicity and decrease gradually. Besides, more asymmetrical motions are put into consideration with regard to tread conicity variation. Similarly, one asymmetrical behavior with small amplitude difference originates from the same chaotic attractor at both sides. And small interactive amplitude jumps in two sides at chaotic or periodic occasions are revealed. Distinguishingly, the other new asymmetrical type is found at a certain tread conicity that amplitudes of the hunting motion in positive and reverse directions no longer coincide and go away in opposite directions when the tread conicity increases to a certain value. And this particular asymmetrical motion disappears with further growth of tread conicity.

     

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.

     

Influences of aerodynamic loads on hunting stability of high-speed railway vehicles and parameter studies

The influences of steady aerodynamic loads on hunting stability of high-speed railway vehicles were investigated in this study.A mechanism is suggested to explain the change of hunting behavior due to actions of aerodynamic loads：the aerodynamic loads can change the position of vehicle system（consequently the contact relations）,the wheel/rail normal contact forces,the gravitational restoring forces/moments and the creep forces/moments.A mathematical model for hunting stability incorporating such influences was developed.A computer program capable of incorporating the effects of aerodynamic loads based on the model was written,and the critical speeds were calculated using this program.The dependences of linear and nonlinear critical speeds on suspension parameters considering aerodynamic loads were analyzed by using the orthogonal test method,the results were also compared with the situations without aerodynamic loads.It is shown that the most dominant factors a ff ecting linear and nonlinear critical speeds are different whether the aerodynamic loads considered or not.The damping of yaw damper is the most dominant influencing factor for linear critical speeds,while the damping of lateral damper is most dominant for nonlinear ones.When the influences of aerodynamic loads are considered,the linear critical speeds decrease with the rise of cross wind velocity,whereas it is not the case for the nonlinear critical speeds.The variation trends of critical speeds with suspension parameters can be significantly changed by aerodynamic loads.Combined actions of aerodynamic loads and suspension parameters also a ff ect the critical speeds.The effects of such joint action are more obvious for nonlinear critical speeds.     

Effect of Lateral Linear Stiffness on Nonlinear Characteristics of Hunting Motion of a Railway Wheelset

Railway wheelset experiences the problem of hunting above a critical speed, which is a kind of self-excited oscillation. At the critical speed, it is known that the system undergoes a subcritical Hopf bifurcation. Therefore, for clarifying the nonlinear characteristics of hunting it is very important to detect, for example, the nonlinear forces in the wheelset due to the creep forces acting between the wheels and rails, and the nonlinear component of the resorting forces by the suspensions. However, it is impossible to determine each force quantitatively. In the present paper, it is first shown, by using the center manifold theory and the method of normal form, that the nonlinear characteristics of the bifurcation in a wheelset model with two degrees of freedom are governed by a single parameter, hence each nonlinear force need not be detected when examining the nonlinear characteristics. Also, a method of determining the governing parameter from experimentally observed radiuses of the unstable limit cycle is proposed. Next, we experimentally investigate the variation of the parameter due to the presence of linear spring suspensions in the lateral direction and discuss the variation of the nonlinear characteristics of the hunting motion, which depends on the lateral stiffness. As a result, the improvement of the stability of the wheelset against the disturbance by the linear spring suspensions is clarified.     

轮对非线性动力学模型稳定性和分岔特性研究

为了探究轮对系统的横向失稳问题,考虑了陀螺效应和一系悬挂阻尼的影响作用,建立非线性轮轨接触关系的轮对动力学模型,研究轮对系统的蛇行稳定性、Hopf分岔特性及迁移转化机理.通过稳定性判据获得了轮对系统失稳临界速度.采用中心流形定理和规范型方法对轮对动力学模型进行化简,得到与轮对系统分岔特性相同的一维复变量方程,理论推导求得轮对系统的第一Lyapunov系数的表达式,根据其符号即可判断轮对系统的Hopf分岔类型.讨论了不同参数对轮对系统Hopf分岔临界速度的影响,探究了轮对系统的超临界、亚临界Hopf分岔域在二维参数空间的分布规律.利用数值模拟得到轮对系统的3种典型Hopf分岔图,验证了轮对系统超临界、亚临界Hopf分岔域分布规律的正确性.结果表明,轮对系统的临界速度随着等效锥度的增大而减小,随着一系悬挂的纵向刚度和纵向阻尼的增大而增大,随着纵向蠕滑系数的增大呈先增大后减小.系统参数变化会引起轮对系统Hopf分岔类型发生改变,即亚临界与超临界Hopf分岔相互迁移转化.轮对系统Hopf分岔域在二维参数空间的分布规律对于轮对系统参数匹配和优化设计具有一定的指导意义.     

Modeling and analysis of piezoaeroelastic energy harvesters

This work investigates the influence of structural and aerodynamic nonlinearities on the dynamic behavior of a piezoaeroelastic system. The system is composed of a rigid airfoil supported by nonlinear torsional and flexural springs in the pitch and plunge motions, respectively, with a piezoelectric coupling attached to the plunge degree of freedom. The analysis shows that the effect of the electrical load resistance on the flutter speed is negligible in comparison to the effects of the linear spring coefficients. The effects of aerodynamic nonlinearities and nonlinear plunge and pitch spring coefficients on the system’s stability near the bifurcation are determined from the nonlinear normal form. This is useful to characterize the effects of different parameters on the system’s output and ensure that subcritical or “catastrophic” bifurcation does not take place. Numerical solutions of the coupled equations for two different configurations are then performed to determine the effects of varying the load resistance and the nonlinear spring coefficients on the limit-cycle oscillations (LCO) in the pitch and plunge motions, the voltage output and the harvested power.     

Two-parameter bifurcation analysis of limit cycles of a simplified railway wheelset model

The effect of the nonlinear terms on bifurcation behaviors of limit cycles of a simplified railway wheelset model is investigated. At first, the stable equilibrium state loses its stability via a Hopf bifurcation. The bifurcation curve is divided into a supercritical branch and a subcritical one by a generalized Hopf point, which plays a key role in determining the occurrence of flange contact and derailment of high-speed railway vehicles, and the occurrence of this critical situation is an important decision-making criteria for design parameters. Secondly, bifurcations of limit cycles are discussed by comparing the bifurcation behavior of cycles for two different nonlinear parameters. Unlike local Hopf bifurcation analysis based on a single bifurcation parameter in most papers, global bifurcation analysis of limit cycles based on two bifurcation parameters is investigated, simultaneously. It is shown that changing nonlinear parameter terms can affect bifurcation types of cycles and division of parameter domains. In particular, near the branch points of cycles, two symmetrical limit cycles are created by a pitchfork bifurcation and then two symmetrical cycles both undergo a period-doubling bifurcation to form two stable period-two cycles. Around the resonant points, period orbits can make several turns, whose number of turns corresponds to the ratio of resonance. Thirdly, near the Neimark–Sacker bifurcation of cycles, a stable torus is created by a supercritical Neimark–Sacker bifurcation, which shows that the orbit of the model exhibits modulated oscillations with two frequencies near the limit cycle. These results demonstrate that nonlinear parameter terms can produce very complex global bifurcation phenomena and make obvious effects on possible hunting motions even though a simple railway wheelset model is concerned.

     

轮对非线性动力学系统蛇行运动的解析解

在对铁路车辆系统的极限环幅值和非线性临界速度进行分析时通常采用数值方法, 不便于研究其随系统参数的变化规律. 轮对系统保留了影响车辆系统动力学性能的几个关键要素: 如轮轨几何非线性约束、轮轨接触蠕滑关系和悬挂系统等, 可以反映铁路车辆系统蛇行运动的本质特性. 轮对系统自由度少、参数少, 可以采用解析方法进行分析. 本文选取合适的特征量把轮对非线性动力学方程无量纲化, 得到了带有小参数的两自由度微分方程; 采用多尺度方法对该方程进行了解析求解; 给出了轮对系统极限环幅值的解析表达式并对其稳定性进行了判定; 给出了轮对系统的分岔速度解析表达式, 并进而获得系统的非线性临界速度的解析表达式. 在对得到的解析解用数值结果进行验证后, 用得到的解析解进行了系统参数影响分析. 传统的分岔图计算方法(如降速法、路径跟踪法等)需对微分方程进行大量数值积分计算方可求解系统的非线性临界速度值, 而通过本文获得的解析表达式可直接给出系统的非线性临界速度值和极限环幅值, 便于研究轮对系统动力学特性随参数的变化规律,进行快速方案比对和筛选, 为转向架结构优化设计提供参考.      

Subcritical Hopf and saddle-node bifurcations in hunting motion caused by cubic and quintic nonlinearities: experimental identification of nonlinearities in a roller rig

Nonlinear Dynamics - Railway vehicles suffer from hunting motion, even when traveling below the critical speed obtained by linear analysis, due to the nonlinear characteristics of the wheel system....     

Control of cross-flow-induced vibrations of square cylinders using linear and nonlinear delayed feedbacks

We investigate the effectiveness of linear and nonlinear time-delay feedback controls to suppress high amplitude oscillations of an elastically mounted square cylinder undergoing galloping oscillations. A representative model that couples the transverse displacement and the aerodynamic force is used. The quasi-steady approximation is used to model the galloping force. A linear analysis is performed to investigate the effect of linear time-delay controls on the onset speed of galloping and natural frequencies. It is demonstrated that a linear time-delay control can be used to delay the onset speed of galloping. The normal form of the Hopf bifurcation is then derived to characterize the type of the instability (supercritical or subcritical) and to determine the effects of the linear and nonlinear time-delay parameters on their outputs near the bifurcation. The results show that the nonlinear time-delay control can be efficiently implemented to significantly reduce the galloping amplitude and suppress any dangerous behavior by converting any subcritical Hopf bifurcation into a supercritical one.     

A discrete approach for a generalized Beck’s column in parametric resonance

A Galerkin projection based on non-standard bases is conceived to derive an equivalent discrete model of a continuous system under non-conservative forces. The problem of deriving a discrete model capable of describing critical and post-critical scenarios for non-selfadjoint systems is discussed and an heuristic rule for a proper choice of trial functions is given. The procedure is utilized to analyze the effect of non-conservative autonomous and non-autonomous (pulsating) forces acting on a linearly damped Beck’s column involving geometrical nonlinearities. The linear and nonlinear behaviours arising from the analysis of the proposed discrete model are in good agreement with those observed through the unavoidably more involved direct continuous approach. Critical scenarios for the autonomous and non-autonomous cases are investigated and the multiple scales method is applied in order to obtain the bifurcation equations in the neighborhood of a Hopf bifurcation point in primary parametric resonance. A comparison between critical and post-critical continuous and discrete model is performed adopting two control parameters, namely the amplitudes of the static and dynamic components of the forces, playing the role of detuning and bifurcation parameters, respectively.     

Integration of uniform design and quantum-behaved particle swarm optimization to the robust design for a railway vehicle suspension system under different wheel conicities and wheel rolling radii

This paper proposes a systematic method, inte-grating the uniform design(UD)of experiments and quantum-behaved particle swarm optimization(QPSO),to solve the problem of a robust design for a railway vehicle suspension system. Based on the new nonlinear creep model derived from combining Hertz contact theory, Kalker's linear the-ory and a heuristic nonlinear creep model,the modeling and dynamic analysis of a 24 degree-of-freedom railway vehi-cle system were investigated.The Lyapunov indirect method was used to examine the effects of suspension parameters, wheel conicities and wheel rolling radii on critical hunting speeds.Generally,the critical hunting speeds of a vehicle sys-tem resulting from worn wheels with different wheel rolling radii are lower than those of a vehicle system having origi-nal wheels without different wheel rolling radii.Because of worn wheels, the critical hunting speed of a running rail-way vehicle substantially declines over the long term. For safety reasons,it is necessary to design the suspension sys-tem parameters to increase the robustness of the system and decrease the sensitive of wheel noises.By applying UD and QPSO,the nominal-the-best signal-to-noise ratio of the sys-tem was increased from?48.17 to?34.05 dB.The rate of improvement was 29.31%.This study has demonstrated that the integration of UD and QPSO can successfully reveal the optimal solution of suspension parameters for solving the robust design problem of a railway vehicle suspension sys-tem.     

Hopf bifurcation analysis of railway bogie

The railway bogie, the most important running component, has direct association with the dynamic performance of the whole vehicle system. The bifurcation type of the bogie that is affected by vehicle parameters will decide the behavior of the vehicle hunting stability. This paper mainly analyzes the effect of the yaw damper and wheel tread shape on the stability and bifurcation type of the railway bogie. The center manifold theorem is adopted to reduce the dimension of the bogie dynamical model, and the symbolic expression for determining the bifurcation type at the critical speed is obtained by the method of normal form. As a result, the influence of yaw damper on the bifurcation type of the bogie is given qualitatively in contrast to typical wheel profiles with high and low wheel tread effective conicities. Besides, the discriminant of bifurcation type for the wheel tread parameter variation is given which depicts the variation tendency of dynamics characteristics. Finally, numerical analysis is given to exhibit corresponding bifurcation diagrams.     

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

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

A universal bifurcation mechanism arising from progressive hydroelastic waves

A unidirectional, weakly dispersive nonlinear model is proposed to describe the supercritical bifurcation arising from hydroelastic waves in deep water. This model equation, including quadratic, cubic, and quartic nonlinearities, is an extension of the famous Whitham equation. The coefficients of the nonlinear terms are chosen to match with the key properties of the full Euler equations, precisely, the associated cubic nonlinear Schr?dinger equation and the amplitude of the solitary wave at the bifurcation point. It is shown that the supercritical bifurcation, rich with Stokes, solitary, generalized solitary, and dark solitary waves in the vicinity of the phase speed minimum, is a universal bifurcation mechanism. The newly developed model can capture the essential features near the bifurcation point and easily be generalized to other nonlinear wave problems in hydrodynamics.     

Post-critical behavior of Beck's column with a tip mass

This study examines how a tip mass with rotary inertia affects the stability of a follower-loaded cantilevered column. Using nonlinear modeling and perturbation analysis, expressions are set up for determining the stability of the straight column and the amplitude of post-critical flutter oscillations. Bifurcation diagrams are given, showing how the vibration amplitude changes with follower load and other parameters. These results agree closely with numerical simulation. It is found that sufficiently large values of tip mass rotary inertia can change the primary bifurcation from supercritical into subcritical. This can imply very large motions for follower loads just beyond critical, contrasting the finite amplitude motions accompanying supercritical bifurcations. Also, the straight column may be destabilized by a sufficiently strong disturbance at loads far below the value of critical load predicted by linear theory. A similar change in bifurcation is found to occur with increased external (as compared to internal) damping, and with a shortening in column length. These effects are not revealed by linear modeling and analysis, which may consequently fail to predict even qualitatively the real critical load for a column with tip mass.     

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一. 目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生. 但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少. 本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌. 考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性 振动的积分--偏微分控制方程. 首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为. 基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混 沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为. 结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为. 随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大. 较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动. 最后,给出时程图、频谱图、相图以及Poincaré 映射图来确定梁的混沌运动.      

轮对系统的Hopf分岔研究

针对轮对系统中的非线性动力学问题, 本文基于Hopf分岔代数判据得到考虑陀螺效应的轮对系统Hopf分岔点解析表达式, 即轮对系统蛇形失稳的线性临界速度解析表达式. 基于分岔理论得到轮对系统的第一、第二Lyapunov系数表达式, 并结合打靶法分别得到不同纵向刚度下, 考虑陀螺效应与不考虑陀螺效应的轮对系统分岔图. 通过对比有无陀螺效应的轮对系统分岔图发现, 在同一纵向刚度下, 考虑陀螺效应的轮对系统线性临界速度和非线性临界速度均大于不考虑陀螺效应的轮对系统, 即陀螺效应可以提高轮对系统的运动稳定性. 基于Bautin分岔理论, 以纵向刚度和纵向速度作为参数, 分别得到考虑陀螺效应和不考虑陀螺效应的轮对系统, 从亚临界Hopf分岔到超临界Hopf分岔, 再从超临界Hopf分岔到亚临界Hopf分岔的迁移机理拓扑图. 通过对比有、无陀螺效应的轮对系统Bautin分岔拓扑图发现, 陀螺效应将改变轮对系统的退化Hopf分岔点, 但对于轮对系统Bautin分岔拓扑图的影响不大.      

Nonlinear parametric modeling of suspension bridges under aeroelastic forces: torsional divergence and flutter

A fully nonlinear model of suspension bridges parameterized by one single space coordinate is proposed to describe overall three-dimensional motions. The nonlinear equations of motion are obtained via a direct total Lagrangian formulation and the kinematics, for the deck-girder and the suspension cables, feature the finite displacements of the associated base lines and the flexural and torsional rotations of the deck cross-sections assumed rigid in their own planes. The strain-displacement relationships for the generalized strain parameters, the elongations in the cables, the deck elongation, and the three curvatures, retain the full geometric nonlinearities. The proposed nonlinear model with its full extensional-flexural-torsional coupling is employed to study the torsional divergence caused by the static part of the wind-induced forces. Two suspension bridges are considered as case studies: the Runyang bridge (main span 1,490?m) and the Hu Men bridge (main span 888?m) in China. The evaluation of the onset of the static instability and the post-critical behavior takes into account the prestressed condition of the bridge subject to dead loads. The dynamic bifurcation that occurs at the onset of flutter is also studied accounting for the prestressed equilibrium state about which the equations of motion are obtained via an updated Lagrangian formulation. Such a bifurcation is investigated in the context of the parametric nonlinear model considering the model parameters of the Runyang Suspension Bridge together with its aeroelastic derivatives. The calculated critical wind speeds for the onset of the static and dynamic bifurcations are compared with the results obtained via linear analysis and the main differences are highlighted. Parametric sensitivity studies are carried out to assess the influence of the design parameters on the instabilities associated with the bridge aeroelastic response.     

Identification of the nonlinear vibration characteristics in hydropower house using transfer entropy

The nonlinear behavior is the primary attribute of the dynamic system stability. In this study, the time-delayed transfer entropy method is proposed to identify the nonlinear dynamic behavior of hydropower house and civil construction, including the transport directionality of information, cracked damage location, and degree of dynamic nonlinearity. Differing from the objects investigated in currently available literature, large-scale civil engineering structures such as hydropower house are more complex and larger size, which stimulate the demand for special identification techniques. Due to the fact that nonlinearity of hydropower house can be induced by many types of interactions between structure and mechanism, a simple similarity model of a multistory building, including damaged contact nonlinearity is studied first following by a discussion of the method for identifying information transmission directionality of the linear or nonlinear structure. The method for identifying the source and degrees of structural nonlinearity vibration is described. Furthermore, the procedure for identification of nonlinearity dynamic behavior in the hydropower house structure based on transfer entropy is studied based on a prototype field experiment under various load cases. Rather than the traditional linear signal processing tools and identification methods, the advantage of this proposed method is to identify the nonlinearity dynamic characteristics of hydropower house structure. This study provides a valuable reference for identifying the damage-induced nonlinearities in civil engineering structures as well as studying the nonlinear dynamic characteristics of hydropower house.