Supercritical as well as subcritical Hopf bifurcation in nonlinear flutter systems

The Hopf bifurcations of an airfoil flutter system with a cubic nonlinearity are investigated, with the flow speed as the bifurcation parameter. The center manifold theory and complex normal form method are Used to obtain the bifurcation equation. Interestingly, for a certain linear pitching stiffness the Hopf bifurcation is both supercritical and subcritical. It is found, mathematically, this is caused by the fact that one coefficient in the bifurcation equation does not contain the first power of the bifurcation parameter. The solutions of the bifurcation equation are validated by the equivalent linearization method and incremental harmonic balance method.     

Nonlinear flutter behavior of a plate with motion constraints in subsonic flow

This paper is aimed at presenting the nonlinear flutter peculiarities of a cantilevered plate with motion-limiting constraints in subsonic flow. A non-smooth free-play structural nonlinearity is considered to model the motion constraints. The governing nonlinear partial differential equation is discretized in space and time domains by using the Galerkin method. The equilibrium points and their stabilities are presented based on qualitative analysis and numerical studies. The system loses its stability by flutter and undergoes the limit cycle oscillations (LCOs) due to the nonlinearity. A heuristic analysis scheme based on the equivalent linearization method is applied to theoretical analysis of the LCOs. The Hopf and two-multiple semi-stable limit cycle bifurcation bifurcations are supercritical or subcritical, which is dependent on the location of the motion constraints. For some special cases the bifurcations are, interestingly, both supercritical and subcritical. The influence of varying parameters on the dynamics is discussed in detail. The results predicted by the analysis scheme are in good agreement with the numerical ones.     

BIFURCATIONS OF AIRFOIL IN INCOMPRESSIBLE FLOW

Bifurcations of an airfoil with nonlinear pitching stiffness in incompressible flow are investigated. The pitching spring is regarded as a spring with cubic stiffness. The motion equations of the airfoil are written as the four dimensional one order differential equations. Taking air speed and the linear part of pitching stiffness as the parameters, the analytic solutions of the critical boundaries of pitchfork bifurcations and Hopf bifurcations are obtained in 2 dimensional parameter plane. The stabilities of the equilibrium points and the limit cycles in different regions of 2 dimensional parameter plane are analyzed. By means of harmonic balance method, the approximate critical boundaries of 2-multiple semi-stable limit cycle bifurcations are obtained, and the bifurcation points of supercritical or subcritical Hopf bifurcation are found. Some numerical simulation results are given.     

Limit cycle oscillation suppression of 2-DOF airfoil using nonlinear energy sink

This paper presents a novel mechanical attachment, i.e., nonlinear energy sink （NES）, for suppressing the limit cycle oscillation （LCO） of an airfoil. The dynamic responses of a two-degree-of-freedom （2-DOF） airfoil coupled with an NES are studied with the harmonic balance method. Different structure parameters of the NES, i.e., mass ratio between the NES and airfoil, NES offset, NES damping, and nonlinear stiffness in the NES, are chosen for studying the effect of the LCO suppression on an aeroelastic system with a supercritical Hopf bifurcation or subcritical Hopf bifurcation, respectively. The results show that the structural parameters of the NES have different influence on the supercritical Hopf bifurcation system and the subcritical Hopf bifurcation system.     

Modified post-bifurcation dynamics and routes to chaos from double-Hopf bifurcations in a hyperchaotic system

In order to understand the onset of hyperchaotic behavior recently observed in many systems, we study bifurcations in the modified Chen system leading from simple dynamics into chaotic regimes. In particular, we demonstrate that the existence of only one fixed point of the system in all regions of parameter space implies that this simple point attractor may only be destabilized via a Hopf or double Hopf bifurcation as system parameters are varied. Saddle-node, transcritical and pitchfork bifurcations are precluded. The normal form immediately following double Hopf bifurcations is constructed analytically by the method of multiple scales. Analysis of this generalized double Hopf normal form along standard lines reveals possible regimes of periodic solutions, two-period tori, and three-period tori in parameter space. However, considering these more carefully, we find that only certain combinations or sequences of these dynamical regimes are possible, while others derived and considered in earlier work are in fact mathematically impossible. We also discuss the post-bifurcation dynamics in the context of two intermittent routes to chaos (routes following either (i) subcritical or (ii) supercritical Hopf or double Hopf bifurcations). In particular, the route following supercritical bifurcations is somewhat subtle. Such behavior following repeated Hopf bifurcations is well-known and widely observed, including in the classical Ruelle?CTakens and quasiperiodic routes to chaos. However, to the best of our knowledge, it has not been considered in the context of the double-Hopf normal form, although it has been numerically observed and tracked in the post-double Hopf regime. Numerical simulations are employed to corroborate these various predictions from the normal form. They reveal the existence of stable periodic and toroidal attractors in the post-supercritical-Hopf cases, and either attractors at infinity or bounded chaotic dynamics following subcritical Hopf bifurcations. Future work will map out the remainder of the routes into the chaotic regimes, including further bifurcations of the post-supercritical-Hopf two- and three-tori via either torus doubling or breakdown.     

Subcritical, nontypical and period-doubling bifurcations of a delta wing in a low speed wind tunnel

Limit Cycle Oscillations (LCOs) involving Delta wings are an important area of research in modern aeroelasticity. Such phenomena can be the result of geometric or aerodynamic nonlinearity. In this paper, a flexible half-span Delta wing is tested in a low speed wind tunnel in order to investigate its dynamic response. The wing is designed to be more flexible than the models used in previous research on the subject in order to expand the airspeed range in which LCOs occur. The experiments reveal that this wing features a very rich bifurcation behavior. Three types of bifurcation are observed for the first time for such an aeroelastic system: subcritical bifurcations, period-doubling/period-halving and nontypical bifurcations. They give rise to a great variety of LCOs, even at very low angles of attack. The LCOs resulting from the nontypical bifurcation display Hopf-type behavior, i.e. having fundamental frequencies equal to one of the linear modal frequencies. All of the other LCOs have fundamental frequencies that are unrelated to the underlying linear system modes.     

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

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

An incremental method for limit cycle oscillations of an airfoil with an external store

This paper proposes an incremental method, which is based on the harmonic balance method, to analyze the nonlinear aeroelastic problem of an airfoil with an external store. The governing equations of limit cycle oscillations (LCOs) of the airfoil are deduced by the harmonic balancing procedure. Different from usual procedures, the harmonic balance equations are not solved directly but instead transformed into an equivalent minimization problem. The minimization problem is solved using the Levenberg–Marquardt method. Numerical examples show that the LCOs obtained by the presented method are in excellent agreement with numerical solutions. The bifurcation of the LCOs is further analyzed using the Floquet theory. It is found that the LCOs exhibit saddle-node, symmetry breaking and period-doubling bifurcations with the wind speed as control parameter. Compared with the harmonic balance method, the presented method has a wider convergence region and hence makes it easier to choose a proper initial guess for iterations.     

轮对系统的Hopf分岔研究

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

Criticality of Hopf bifurcation in state-dependent delay model of turning processes

In this paper the non-linear dynamics of a state-dependent delay model of the turning process is analyzed. The size of the regenerative delay is determined not only by the rotation of the workpiece, but also by the vibrations of the tool. A numerical continuation technique is developed that can be used to follow the periodic orbits of a system with implicitly defined state-dependent delays. The numerical analysis of the model reveals that the criticality of the Hopf bifurcation depends on the feed rate. This is in contrast to simpler constant delay models where the criticality does not change. For small feed rates, subcritical Hopf bifurcations are found, similar to the constant delay models. In this case, periodic orbits coexist with the stable stationary cutting state and so there is the potential for large amplitude chatter and bistability. For large feed rates, the Hopf bifurcation becomes supercritical for a range of spindle speeds. In this case, stable periodic orbits instead coexist with the unstable stationary cutting state, removing the possibility of large amplitude chatter. Thus, the state-dependent delay in the model has a kind of stabilizing effect, since the supercritical case is more favorable from a practical viewpoint than the subcritical one.     

An analytical study of time-delayed control of friction-induced vibrations in a system with a dynamic friction model

We investigate the control of friction-induced vibrations in a system with a dynamic friction model which accounts for hysteresis in the friction characteristics. Linear time-delayed position feedback applied in a direction normal to the contacting surfaces has been employed for the purpose. Analysis shows that the uncontrolled system loses stability via. a subcritical Hopf bifurcation making it prone to large amplitude vibrations near the stability boundary. Our results show that the controller achieves the dual objective of quenching the vibrations as well as changing the nature of the bifurcation from subcritical to supercritical. Consequently, the controlled system is globally stable in the linearly stable region and yields small amplitude vibrations if the stability boundary is crossed due to changes in operating conditions or system parameters. Criticality curve separating regions on the stability surface corresponding to subcritical and supercritical bifurcations is obtained analytically using the method of multiple scales (MMS). We have also identified a set of control parameters for which the system is stable for lower and higher relative velocities but vibrates for the intermediate ones. However, the bifurcation is always supercritical for these parameters resulting in low amplitude vibrations only.     

Criticality of bifurcation in the tuned axial–torsional rotary drilling model

We study the bifurcation characteristics of a lumped-parameter model of rotary drilling with 1:1 internal resonance between the axial and the torsional modes which leads to the largest stability thresholds. For this special case, the two-degree-of-freedom model for the drill-string reduces to an effectively single-degree-of-freedom system facilitating further analysis. The regenerative effect of the cutting action due to the axial vibrations is incorporated through a delayed term in the cutting force with the delay depending on the torsional oscillations. This state dependency of the delay introduces nonlinearity in the current model. Steady drilling loses stability via a Hopf bifurcation, and the nature of the bifurcation is determined by an analytical study using the method of multiple scales. We find that both subcritical and supercritical Hopf bifurcations are present in this system depending on the choice of operating parameters. Hence, the nonlinearity due to the state-dependent delay term could both be stabilizing or destabilizing in nature, and the self-interruption nonlinearity is essential to capture the global behavior. Numerical bifurcation analysis of a global axial–torsional model of rotary drilling further confirms the analytical results from the method of multiple scales. Further exploration of the rotary drilling dynamics unravels more complex phenomena including grazing bifurcations and possibly chaotic solutions.     

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.

     

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.     

Design of piezoaeroelastic energy harvesters

We design a piezoaeroelastic energy harvester consisting of a rigid airfoil that is constrained to pitch and plunge and supported by linear and nonlinear torsional and flexural springs with a piezoelectric coupling attached to the plunge degree of freedom. We choose the linear springs to produce the minimum flutter speed and then implement a linear velocity feedback to reduce the flutter speed to any desired value and hence produce limit-cycle oscillations at low wind speeds. Then, we use the center-manifold theorem to derive the normal form of the Hopf bifurcation near the flutter onset, which, in turn, is used to choose the nonlinear spring coefficients that produce supercritical Hopf bifurcations and increase the amplitudes of the ensuing limit cycles and hence the harvested power. For given gains and hence reduced flutter speeds, the harvested power is observed to increase, achieve a maximum, and then decrease as the wind speed increases. Furthermore, the response undergoes a secondary supercritical Hopf bifurcation, resulting in either a quasiperiodic motion or a periodic motion with a large period. As the wind speed is increased further, the response becomes eventually chaotic. These complex responses may result in a reduction in the generated power. To overcome this adverse effect, we propose to adjust the gains to increase the flutter speed and hence push the secondary Hopf bifurcation to higher wind speeds.     

Hopf bifurcation analysis of asymmetrical rotating shafts

The Hopf and double Hopf bifurcations analysis of asymmetrical rotating shafts with stretching nonlinearity are investigated. The shaft is simply supported and is composed of viscoelastic material. The rotary inertia and gyroscopic effect are considered, but, shear deformation is neglected. To consider the viscoelastic behavior of the shaft, the Kelvin–Voigt model is used. Hopf bifurcations occur due to instability caused by internal damping. To analyze the dynamics of the system in the vicinity of Hopf bifurcations, the center manifold theory is utilized. The standard normal forms of Hopf bifurcations for symmetrical and asymmetrical shafts are obtained. It is shown that the symmetrical shafts have double zero eigenvalues in the absence of external damping, but asymmetrical shafts do not have. The asymmetrical shaft in the absence of external damping has a saddle point, therefore the system is unstable. Also, for symmetrical and asymmetrical shafts, in the presence of external damping at the critical speeds, supercritical Hopf bifurcations occur. The amplitude of periodic solution due to supercritical Hopf bifurcations for symmetrical and asymmetrical shafts for the higher modes would be different, due to shaft asymmetry. Consequently, the effect of shaft asymmetry in the higher modes is considerable. Also, the amplitude of periodic solutions for symmetrical shafts with rotary inertia effect is higher than those of without one. In addition, the dynamic behavior of the system in the vicinity of double Hopf bifurcation is investigated. It is seen that in this case depending on the damping and rotational speed, the sink, source, or saddle equilibrium points occur in the system.     

Transonic limit cycle oscillation behavior of an aeroelastic airfoil with free-play

The limit cycle oscillation (LCO) behaviors of an aeroelastic airfoil with free-play for different Mach numbers are studied. Euler equations are adopted to obtain the unsteady aerodynamic forces. Aerodynamic and structural describing functions are employed to deal with aerodynamic and structural nonlinearities, respectively. Then the flutter speed and flutter frequency are obtained by V-g method. The LCO solutions for the aeroelastic airfoil obtained by using dynamically linear aerodynamics agree well with those obtained directly by using nonlinear aerodynamics. Subsequently, the dynamically linear aerodynamics is assumed, and results show that the LCOs behave variously in different Mach number ranges. A subcritical bifurcation, consisting of both stable and unstable branches, is firstly observed in subsonic and high subsonic regime. Then in a narrow Mach number range, the unstable LCOs with small amplitudes turn to be stable ones dominated by the single degree of freedom flutter. Meanwhile, these LCOs can persist down to very low flutter speeds. When the Mach number is increased further, the stable branch turns back to be unstable. To address the reason of the stability variation for different Mach numbers at small amplitude LCOs, we find that the Mach number freeze phenomenon provides a physics-based explanation and the phase reversal of the aerodynamic forces will trigger the single degree of freedom flutter in the narrow Mach number range between the low and high Mach numbers of the chimney region. The high Mach number can be predicted by the freeze Mach number, and the low one can be estimated by the Mach number at which the aerodynamic center of the airfoil lies near its elastic axis. Influence of angle of attack and viscous effects on the LCO behavior is also discussed.     

随机干扰与随机参数激励联合作用下的Hopf分叉

本文研究van der Pol-Duffing型的非线性振子在随机干扰和随机参数联合作用下的Hopf分叉现象。本文所得结果证实了当系统处在于Hopf分叉点附近时,对系统的参数的变化具有敏感性。在研究过程中,我们利用Markov扩散过程逼近系统的随机响应,得到了沿稳定矩的概率1稳定和矩稳定的条件。对于非线性振子,我们得到了振幅过程的稳态概论密度函数。研究发现,确定性系统的Hopf分叉点在随机参数作用下具有漂移现象,这种漂移是由系统的性质所决定的,当分叉点为超临界的,分叉点向前漂移;而当分叉点为亚临界时,这种漂移是向后的。当系统处在外部随机干扰作用下时,系统出现非零响应。另外我们发现,稳态矩的分叉与其阶数无关。     

Periodic Delay Effects on Cutting Dynamics

We consider a delay equation whose delay is perturbed by a small periodic fluctuation. In particular, it is assumed that the delay equation exhibits a Hopf bifurcation when its delay is unperturbed. The periodically perturbed system exhibits more delicate bifurcations than a Hopf bifurcation. We show that these bifurcations are well explained by the Bogdanov-Takens bifurcation when the ratio between the frequencies of the periodic solution of the unperturbed system (Hopf bifurcation) and the external periodic perturbation is 1:2. Our analysis is based on center manifold reduction theory.     

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.