Hopf Bifurcation and Hunting Behavior in a Rail Wheelset with Flange Contact

An analytical investigation of Hopf bifurcation and hunting behavior of a rail wheelset with nonlinear primary yaw dampers and wheel-rail contact forces is presented. This study is intended to complement earlier studies by True et al., where they investigated the nonlinearities stemming from creep-creep force saturation and nonlinear contacts between a realistic wheel and rail profile. The results indicate that the nonlinearities in the primary suspension and flange contact contribute significantly to the hunting behavior. Both the critical speed and the nature of bifurcation are affected by the nonlinear elements. Further, the results show that in some cases, the critical hunting speed from the nonlinear analysis is less than the critical speed from a linear analysis. This indicates that a linear analysis could predict operational speeds that in actuality include hunting.     

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.     

Passive control of vibration of thin-walled gears: advanced modelling of ring dampers

Vibrations on gears are mainly induced by the gear mesh contact. Resonance conditions of the gear may occur during service if the mesh frequency is close to the natural frequencies of the system at the designed speed of the shaft. Since detuning is not always possible in gears, the response level must be reduced by increasing the damping of the system. In this paper, a passive approach based on the application of a ring damper to reduce the vibration level is presented. The ring damper is placed in a groove underneath the outer rim of the gear. The contact is guaranteed by the preload due to the elasticity of the ring damper itself and above all by the centrifugal force that presses the damper against the groove during rotation. The relative motion of the two components at the contact interface dissipates energy by friction, and hence damping is generated. The vibration amplitude is reduced by optimizing the material and geometrical properties of the ring damper. One of the most important parameters in the determination of the amount of damping due to friction phenomena is the static normal load at the contact, which depends on the mass, the shape, and the material of the ring damper. A numerical method is presented, which couples the static and dynamic equilibrium equations of the assembly. The core of the proposed method is the contact element that takes into account local stick–slip–lift off of the contact and determines the contact forces in terms of static and dynamic loads, which are then used to solve the coupled static and dynamic equilibrium. Since the ring damper has a cut that breaks its continuous circular shape in order to be fitted on the groove, the hypothesis of cyclic symmetry for the gear/ring–damper assembly fails. As a consequence, an appropriate reduced-order modeling is presented to allow the forced response calculations. The algorithm is applied to a dummy bevel gear and to a ring damper having a flat punch contact area. The forced response calculations are performed to highlight the nonlinear interaction between the gear and damper by varying the parameters that mainly affect the amount and distribution of the contact forces and therefore the response level.     

Dynamics of a spacecraft with large flexible appendage constrained by multi-strut passive damper

This paper is concerned with the dynamics of a spacecraft with multi-strut passive damper for large flexible appendage.The damper platform is connected to the spacecraft by a spheric hinge,multiple damping struts and a rigid strut.The damping struts provide damping forces while the rigid strut produces a motion constraint of the multibody system.The exact nonlinear dynamical equations in reducedorder form are firstly derived by using Kane’s equation in matrix form.Based on the assumptions of small velocity and small displacement,the nonlinear equations are reduced to a set of linear second-order differential equations in terms of independent generalized displacements with constant stiffness matrix and damping matrix related to the damping strut parameters.Numerical simulation results demonstrate the damping effectiveness of the damper for both the motion of the spacecraft and the vibration of the flexible appendage,and verify the accuracy of the linear equations against the exact nonlinear ones.     

粘弹性夹芯层合旋转圆板的气动弹性动力稳定性分析

基于哈密尔顿原理,考虑作用于圆板上的空气动力载荷以及粘弹性复模量本构模型,本文建立了含有粘弹性夹芯层的旋转圆板的空气动力学模型,并给出了相应的稳定性的控制方程及边界条件。运用Galerkin方法,文中数值求解并给出了夹芯层合旋转圆板的动力学基本特征,以及结构几何参数、材料参数对于其前、后行波振动频率、阻尼和动力稳定性的影响。结果表明：通过优化合理选取粘弹性夹芯层合旋转圆板的几何参数和材料参数,可提高旋转圆板的临界转速和颤振速度,从而增强其稳定性。     

Nonlinear dynamic responses of high-speed railway vehicles under combined self-excitation and forced excitation considering the influence of unsteady aerodynamic loads

The dynamic characteristics of a railway vehicle system under unsteady aerodynamic loads are examined in this study. A dynamic analysis model of the railway vehicle considering the influences of aerodynamic loads was established. The model not only considers the forced excitation effect of unsteady aerodynamic loads but also accounts for the effect of unsteady aerodynamic loads on the change of the wheel–rail contact normal forces as well as changes of the wheelset creep coefficients and creep forces/moments. Therefore, this model also considers the influences of unsteady aerodynamic loads on the self-excited vibration characteristics of the vehicle system. The time-history curves, phase trajectory diagrams, Poincaré sections, and Lyapunov exponents of the vehicle system running on a smooth straight track under unsteady aerodynamic loads were determined. The results show that when the critical speed is exceeded, the vehicle system usually performs quasi-periodic motion under unsteady aerodynamic loads, which is significantly different from the periodic motion under steady aerodynamic loads. In different cases, the amplitude and phase of motion are significantly different. The amplitude of the motions can be increased by more than 159%, and the difference of phase can be up to 173°. (The phase is almost reversed.) The dynamic responses of the vehicle system under unsteady aerodynamic loads contain abundant frequency components, including the frequency of the self-excited vibration, the frequency of the forced excitation, and combinations of their integer multiples. The vibration forms corresponding to the main harmonic components under unsteady and steady aerodynamic loads were compared, and the self-excited vibration component of the vehicle system under unsteady aerodynamic loads was identified. The variations in the critical speed with various parameter combinations were computed. The variation range of the critical velocity can reach 73%.

     

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.     

Rheological–dynamical analogy: Prediction of damping parameters of hysteresis damper

This research aims to predict the damping parameters of hysteresis damper based on an analytical rheological–dynamical (RDA) visco-elasto-plastic solution of one-dimensional longitudinal continuous vibrations of a bar. A visco-elasto-plastic bar or damper is an energy dissipation device. An attempt is made to estimate quantitatively the influence of material physical parameters of materials on the damping ratio in both the linear visco-elastic analysis and the nonlinear visco-elasto-plastic analysis of damper subjected to external vibration forces. Two types of damping are considered: viscous damping in the case of linear analysis, defined as stiffness and/or mass proportional and, in the case of nonlinear analysis, hysteresis damping caused by inelastic deformations of damper. Owing to the visco-elastic nature of the materials of the damper and the frequency dependence of the viscous damping ratio ξ, it is useful to consider separately the situations arising when ξ is positive (the system is stable) and when it is negative. A negative damping ratio means that the complementary solution of the response would not die away (the system is unstable because of factor e^{ξ · ω · t}). In the case of nonlinear analysis, the force–displacement relation is nonlinear, so it is very difficult to predict the actual damping and stiffness coefficients, even if the force–displacement characteristic is simply perfect elasto-plastic. Using the RDA method, which takes into account the rate of release of visco-elasto-plastic energy of the dissipation devices; nonlinear behaviors are linearized, enabling to obtain the equivalent damping and stiffness coefficients and the effective period for the damper.     

Nonlinear Stiffness of a Magneto-Rheological Damper

The last decade has witnessed an important role of magneto-rheological dampers in the semi-active vibration control on the basis of empirical models. Those models established by fitting experimental data, however, do not offer any explicit expressions for the stiffness and the damping of magneto-rheological dampers. Hence, it is not easy for engineers to get any intuitive information about the effects of stiffness and damping of a magneto-rheological damper on the dynamic performance of a controlled system. To manifest the nonlinear properties of a magneto-rheological damper, this paper presents the hysteretic phenomena and the additional nonlinear stiffness of a typical magneto-rheological damper in terms of equivalent linear stiffness and equivalent linear damping. Then, it gives a brief discussion about the effect of nonlinear stiffness on the vibration control through the numerical simulations and an experiment for the semi-active suspension of a quarter car model with a magneto-rheological damper installed. Both numerical simulations and experimental results show that the additional nonlinear stiffness in the magneto-rheological damper is remarkable, and should be taken into consideration in the design of vibration control.     

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.     

磁流变阻尼器汽车悬架的非线性振动特性分析

通过对磁流变阻尼器sigmoid模型进行数值计算分析,揭示了各参数对该模型阻尼力的影响,证明了通过合理有效地选择参数可设计出符合工程实际应用的磁流变阻尼器.对确定性激励下磁流变阻尼器的Sigmoid模型进行了理论分析.指出了采用非线性的阻尼器模型能更精确地描述汽车悬架的振动特性.从而证明了磁流变阻尼器有良好的使用价值和Sigmoid模型磁流变阻尼器在汽车控制工程应用中的可行性.     

Effects of nonlinear damping suspension on nonperiodic motions of a flexible rotor in journal bearings

Nonlinear damping suspension is a promising method to be used in a rotor-bearing system for vibration isolation between the bearing and environment. However, the nonlinearity of the suspension may influence the stability of the rotor-bearing system. In this paper, the motions of a flexible rotor in short journal bearings with nonlinear damping suspension are studied. A computational method is used to solve the equations of motion, and the bifurcation diagrams, orbits, Poincaré maps, and amplitude spectra are used to display the motions. The results show that the effect of the nonlinear damping suspension on the motions of the rotor-bearing system depends on the speed of rotor: (a) For low speeds, the rotor- bearing system presents the same motion pattern under the nonlinear damping ( \(p=0.5, 2, 3\) ) suspension as for the linear damping ( \(p=1\) ) suspension; (b) For high speeds, the effect of nonlinear damping depends on a combination of the damping exponent and damping coefficient. The square root damping model ( \(p=0.5\) ) shows a wider stable speed range than the linear damping for large damping coefficients. The quadratic damping ( \(p=2\) ) shows similar results to linear damping with some special damping coefficients. The cubic damping ( \(p=3\) ) shows more stable response than the linear damping in general.     

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

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

黏弹性环形板的临界载荷及动力稳定性

利用线性黏弹性力学的Boltzmann叠加原理,在考察位移单值性条件的基础上,给出黏弹性环形板非线性动力学分析的初边值问题。通过Galerkin方法和引进新的状态变量,将其化归为四维非线性非自治常微分方程组,从而得到黏弹性环形板的四种临界载荷,同时考察了几何缺陷对黏弹性薄板临界载荷的影响。根据Floquet理论,得出黏弹性形板在周期激励下的线性动力稳定性判据。综合使用非线性动力学中的数值分析方法,研究了参数对黏弹性环形板非线性动力稳定性的影响。     

Aeroelastic dynamics stability of rotating sandwich annular plate with viscoelastic core layer

A dynamic model for a rotating sandwich annular plate with a viscoelastic core layer is developed. All fundamental equations and boundary conditions are established based on Hamilton’s principle, and the rotation effect and viscoelastic properties of the sandwich structure are taken into account. The aerodynamics force acting on the plate is described by a rotating damping model, and the constitutive behavior of the viscoelastic core layer is formulated by the frequency-dependent complex modulus. The effects of geometrical and material parameters on frequencies and damping of forward and backward traveling waves and the dynamic stability for the rotating sandwich plate are numerically analyzed by means of Galerkin’s method. The results show that the critical and flutter speeds of the rotating plate can be increased at some certain parameters of the viscoelastic core layer.     

Transverse vibrations of rotating shafts: Probability density and first-passage time of whirl radius

Transverse vibrations are considered for a single mass/two-degrees-of-freedom rotating shaft with linear internal or “rotating” damping and nonlinear external damping. The shaft is excited by external random forces. Analysis of resulting random vibrations is based on stochastic averaging method which yields separated (in the linear approximation) equations for complex amplitudes of forward and backward whirling motions. The former of these motions is shown to be dominant at rotation speeds in the vicinity of the instability threshold. Using this approximation an analytical solution is obtained for probability density of squared radius of the shaft's whirl. This solution can be used to detect on-line shaft's instability from its observed response. Solution is also obtained for expected time for reaching given level by the squared whirl radius of the shaft.     

Nonlinear airship aeroelasticity

The aeroelastic derivatives for today's aircraft are calculated in the concept phase using a standard procedure. This scheme has to be extended for large airships, due to various nonlinearities in structural and aerodynamic behaviour. In general, the structural model of an airship is physically as well as geometrically nonlinear. The main sources of nonlinearity are large deformations and the nonlinear material behaviour of membranes. The aerodynamic solution is also included in the nonlinear problem, because the deformed airship influences the surrounding flow. Due to these nonlinearities, the aeroelastic problem for airships can only be solved by an iterative procedure. As one possibility, the coupled aerodynamic and structural dynamic problem was handled using linked standard solvers. On the structural side, the Finite-Element program package ABAQUS was extended with an interface to the aerodynamic solver VSAERO. VSAERO is based on the aerodynamic panel method using potential flow theory. The equilibrium of the internal structural and the external aerodynamic forces leads to the structural response and a trimmed flight state for the specified flight conditions (e.g. speed, altitude). The application of small perturbations around a trimmed state produces reaction forces and moments. These constraint forces are then transferred into translational and rotational acceleration fields by performing an inertia relief analysis of the disturbed structural model. The change between the trimmed flight state and the disturbed one yields the respective aeroelastic derivatives. By including the calculated derivatives in the linearised equation of motion system, it is possible to judge the stability and controllability of the investigated airship.     

高速铁路连续箱型梁桥的动力响应与减振分析

根据列车具体的轴距和轴重,建立了和谐号动车组CRH380AL型列车简化模型;对高速铁路两跨连续梁桥采用多自由度欧拉伯努利梁单元进行主梁的模拟,并将液体黏滞阻尼器模拟为有限元阻尼单元;采用Newmark直接积分方法求解了高速列车作用下的连续梁桥运动方程,数值分析了列车车速以及液体黏滞阻尼器的阻尼系数对于高速铁路连续梁桥振动响应的影响。结果表明:黏滞阻尼器对于桥梁具有明显的减振效果,阻尼力不仅与阻尼系数有关还与列车时速有关;同一黏滞阻尼器条件下,桥梁的最大加速度并不随列车速度的增加而单调增加,而是在某些特定列车车速下桥梁的最大加速度出现了峰值,且随着黏滞阻尼器的阻尼系数增大,桥梁振动响应峰值处的最大加速度减幅不同;同一列车时速的条件下,桥梁的减振效果并不是随着阻尼系数的递增呈正比递增,而是随着阻尼系数的增大,阻尼器的减振效果增幅在减小。     

A numerical and experimental hybrid approach for the investigation of aerodynamic forces on stay cables suffering from rain-wind induced vibration

The aerodynamic forces on a stay cable under a rain-wind induced vibration (RWIV) are difficult to measure directly in a wind tunnel test. This paper presents a hybrid approach that combines an experiment with computational fluid dynamics (CFD) for the investigation on aerodynamic forces of a stay cable under a RWIV. The stay cable and flow field were considered as two substructures of the system. The oscillation of the stay cable was first measured by using a wind tunnel test of a RWIV under an artificial rainfall condition. The oscillation of the cable was treated as a previously known moving boundary condition and applied to the flow field. Only the flow field with the known moving cable boundary was then numerically simulated by using a CFD method (such as Fluent 6.3). The transient aerodynamic forces of the stay cable with a predetermined cable oscillation were obtained from numerical calculations. The characteristics of the aerodynamic forces in the time domain and frequency domain were then analysed for various cases. To verify the feasibility and accuracy of the proposed hybrid approach, the transient aerodynamic forces were applied to a single-degree-of-freedom model (SDOF) of the stay cable to calculate the RWIV of the cable. A comparison was performed between the oscillation responses of the stay cable obtained from the calculated (SDOF model) and experimental results, and the results indicate that the hybrid approach accurately simulates the transient aerodynamic forces of the stay cable. The equivalent damping ratios induced by the aerodynamic forces were obtained for various wind speeds. Furthermore, a nonlinear model of the aerodynamic force is proposed based on the calculation results, and the coefficients in the model were identified by a nonlinear least-squares technique.     

Bifurcation and pull-in voltages of primary resonance of electrostatically actuated SWCNT cantilevers to include van der Waals effect

In this work the voltage response of primary resonance of electrostatically actuated single wall carbon nano tubes (SWCNT) cantilevers over a parallel ground plate is investigated. Three forces act on the SWCNT cantilever, namely electrostatic, van der Waals and damping. While the damping is linear, the electrostatic and van der Waals forces are nonlinear. Moreover, the electrostatic force is also parametric since it is given by AC voltage. Under these forces the dynamics of the SWCNT is nonlinear parametric. The van der Waals force is significant for values less than 50 nm of the gap between the SWCNT and the ground substrate. Reduced order model method (ROM) is used to investigate the system under soft excitation and weak nonlinearities. The voltage-amplitude response and influences of parameters are reported for primary resonance (AC near half natural frequency).