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.     

Homogeneous Equilibrium Configurations of a Hyperelastic Compressible Cube under Equitriaxial Dead-Load Tractions

Non-uniqueness, bifurcation and stability of homogeneous solutions to the equilibrium problem of a hyperelastic cube subject to equitriaxial dead-load tractions are investigated. Besides the basic and theoretical questions raised by the analysis, the study is motivated by the somewhat surprising feature of this nonlinear problem for which the symmetric load may give rise to asymmetric stable deformations. In reality, the equilibrium problem, formulated for general homogeneous compressible isotropic materials with polyconvex energy function, may exhibit primary and secondary bifurcations. A primary bifurcation occurs when there exist paths of equilibrium states that bifurcate from the primary path of three equal principal stretches. These bifurcation branches have two coinciding stretches and along them, through secondary bifurcations, other completely asymmetric bifurcation branches, which are characterized by all three stretches different, may risen. In this case, the cube transforms into an oblique parallelepiped. With increasing loads, they are also possible discontinuous paths of equilibria which evince prompt jumps in the deformation process. Of course, the set of asymmetric solutions admitted by the equilibrium problem depends on the specific form of the stored energy function adopted. In this paper, expressions governing the global development of asymmetric equilibrium branches are derived. In particular, conditions to have bifurcation points are individualized. For compressible neo-Hookean and Mooney-Rivlin materials a wide parametric analysis is carried out showing by means of graphs the most interesting branches. Finally, using the energy criterion, a detailed study is performed to assess the stability of the computed solutions.      

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

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

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.     

Global and bifurcation analysis of a structure with cyclic symmetry

This article combines the application of a global analysis approach and the more classical continuation, bifurcation and stability analysis approach of a cyclic symmetric system. A solid disc with four blades, linearly coupled, but with an intrinsic non-linear cubic stiffness is at stake. Dynamic equations are turned into a set of non-linear algebraic equations using the harmonic balance method. Then periodic solutions are sought using a recursive application of a global analysis method for various pulsation values. This exhibits disconnected branches in both the free undamped case (non-linear normal modes, NNMs) and in a forced case which shows the link between NNMs and forced response. For each case, a full bifurcation diagram is provided and commented using tools devoted to continuation, bifurcation and stability analysis.     

Bifurcation of response of a nonlinear viscoelastic spherical membrane

The quasistatic inflation of a nonlinear viscoelastic spherical membrane by monotonically increasing pressure is considered. The deformation is assumed to be spherically symmetric. For the constitutive equation assumed, circumstances are shown to exist when the radius history must either have a jump discontinuity or bifurcate. A necessary condition for bifurcation and its dependence on material properties and radius history is analysed. Examples of bifurcation for various pressure histories are presented. Post-bifurcation branches are constructed and the possibility of secondary bifurcation is discussed.     

Finite element analysis of post-buckling dynamics in plates. Part II: A non-stationary analysis

With the secondary bifurcation and the local post-secondary buckling behavior being analyzed in Part I, Part II of this study consists of developing an adaptive non-stationary load sweeping algorithm to investigate post-buckling dynamics and mode jumping phenomena of generally (mechanically and thermally) loaded thin plates in a global context. The non-stationary sweeping procedure has the merits of adapting large load steps to capture static characteristics of stable equilibrium paths both before and after mode jumping and reduce automatically the step size to ensure a dynamic transition between the two stable branches. Thus, it is computationally effective. Furthermore, by adopting the non-stationary sweeping scheme, this procedure can avoid spurious convergence of the transient response to an unstable equilibrium.Corresponding to different post-secondary bifurcation forms, which are determined using asymptotical finite element analysis developed in Part I, subsequent buckling patterns of various complexity occurring after mode jumping are obtained using the method developed in this article. Qualitative changes in post-buckled patterns are observed after the occurrence of the secondary bifurcation or the mode jumping. Free vibration analysis using the tangent stiffness matrix obtained from the converged static or dynamic solutions shows a vibration modal shifting phenomena occurs during the process of the load sweep. The spurious convergence phenomenon caused by the application of the traditional hybrid static–dynamic method is found and explained.     

Codimension two bifurcation of a vibro-bounce system

A three-degree-of-freedom vibro-bounce system is considered. The disturbed map of period one single-impact motion is derived analytically. A center manifold theorem technique is applied to reduce the Poincaré map to a three-dimensional one, and the normal form map associated with Hopf-flip bifurcation is obtained. Dynamical behavior of the system, near the point of codimension two bifurcation, is investigated by using qualitative analysis and numerical simulation. It is found that near the point of Hopf-flip bifurcation there exists not only Hopf bifurcation of period one single-impact motion, but also Hopf bifurcation of period two double-impact motion. The results from simulation show that there exists an interesting torus doubling bifurcation near the codimension two bifurcation. The torus doubling bifurcation makes the quasi-periodic attractor associated with period one single-impact motion transform to the other quasi-periodic attractor represented by two attracting closed circles. The torus bifurcation is qualitatively different from the typical torus doubling bifurcation occurring in the vibro-impact systems. Different routes from period one single-impact motion to chaos are observed by numerical simulation.The project supported by the National Natural Science Foundation of China (10172042, 50475109) and the Natural Science Foundation of Gansu Province Government of China (ZS-031-A25-007-Z (key item))     

Global analysis of secondary bifurcation of an elastic bar

In a three dimensional framework of finite deformation configurations, this paper investigates the secondary bifurcation of a uniform, isotropic and linearly elastic bar under compression in a large range of parameters. The governing differential equations and finite dimensional equations of this problem are discussed. It is found that, for a bar with two ends hinged, usually many secondary bifurcation points appear on the primary branches which correspond to the maximum bending stiffness. Results are shown on parameter charts. Secondary modes and branches are also calculated with numerical methods. The project supported in part by the National Natural Science Foundation of China     

Bifurcation and chaos of airfoil with multiple strong nonlinearities

The bifurcation and chaos phenomena of two-dimensional airfoils with multiple strong nonlinearities are investigated.First,the strongly nonlinear square and cubic plunging and pitching stiffness terms are considered in the airfoil motion equations,and the fourth-order Runge-Kutta simulation method is used to obtain the numerical solutions to the equations.Then,a post-processing program is developed to calculate the physical parameters such as the amplitude and the frequency based on the discrete numerical solutions.With these parameters,the transition of the airfoil motion from balance,period,and period-doubling bifurcations to chaos is emphatically analyzed.Finally,the critical points of the period-doubling bifurcations and chaos are predicted using the Feigenbaum constant and the first two bifurcation critical values.It is shown that the numerical simulation method with post-processing and the prediction procedure are capable of simulating and predicting the bifurcation and chaos of airfoils with multiple strong nonlinearities.     

Adapting simulation codes to compute unstable solutions and bifurcation behaviour

Simulation codes for solving large systems of ordinary differential equations suffer from the disadvantage that bifurcation‐theoretic results about the underlying dynamical system cannot be obtained from them easily, if at all. Bifurcation behaviour typically can be inferred only after significant computational effort, and even then the exact location and nature of the bifurcation cannot always be determined definitively. By incorporating relatively minor changes to an existing simulation code for the Taylor–Couette problem, specifically, by implementing the Newton–Picard method, we have developed a computational structure that enables us to overcome some of the inherent limitations of the simulation code and begin to perform bifurcation‐theoretic tasks. While a complete bifurcation picture was not developed, three distinct solution branches of the Taylor–Couette problem were analysed. These branches exhibit a wide variety of behaviours, including Hopf bifurcation points, symmetry‐breaking bifurcation points, turning points and bifurcation to motion on a torus. Unstable equilibrium and time‐periodic solutions were also computed. Copyright © 2007 John Wiley & Sons, Ltd.     

Stability and bifurcation analysis of a nonlinear transversally loaded rotating shaft

The dynamic stability and self-excited posteritical whirling of rotating transversally loaded shaft made of a standard material with elastic and viscous nonlinearities are analyzed in this paper using the theory of bifurcations as a mathematical tool. Partial differential equations of motion are derived under assumption that von Karman's nonlinearity is absent but geometric curvature nonlinearity is included. Galerkin's first-mode discretization procedure is then applied and the equations of motion are transformed to two third-order nonlinear equations that are analyzed using the theory of bifurcation. Condition for nontrivial equilibrium stability is determined and a bifurcating periodic solution of the second-order approximation is derived. The effects of dimensionless stress relaxation time and cubic elastic and viscous nonlinearities as well as the role of the transverse load are studied in the exemplary numerical calculations. A strongly stabilizing influence of the relaxation time is found that may eliminate self-excited vibration at all. Transition from super- to subcritical bifurcation is observed as a result of interaction between system nonlinearities and the transverse load.     

Subharmonic resonance of a clamped-clamped buckled beam with 1:1 internal resonance under base harmonic excitations

The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated. The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton's principle. A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method. A high-dimensional model of the buckled beam is derived, concerning nonlinear coupling. The incremental harmonic balance (IHB) method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve, and the Floquet theory is used to analyze the stability of the periodic solutions. Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited. Bifurcations including the saddle-node, Hopf, perioddoubling, and symmetry-breaking bifurcations are observed. Furthermore, quasi-periodic motion is observed by using the fourth-order Runge-Kutta method, which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.     

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.     

Three routes to chaos for a three-degree-of-freedom articulated cylinder system subjected to annular flow and impacting on the outer pipe

In this paper, the dynamics of a cantilevered articulated system of rigid cylinders interconnected by rotational springs, within a pipe containing fluid flow is studied. Although the formulation is generalized to any number of degrees-of-freedom (articulations), the present work is restricted to three-degree-of-freedom systems. The motions are considered to be planar, and the equations of motion, apart from impacting terms, are linearized. Impacting of the articulated cylinder system on the outer pipe is modelled by either a cubic spring (for analytical convenience) or, more realistically, by a trilinear spring model. The critical flow velocities, for which the system loses stability, by flutter (Hopf bifurcation) or divergence (pitchfork bifurcation) are determined by an eigenvalue analysis. Beyond these first bifurcations, it is shown that, for different values of the system parameters, chaos is obtained through three different routes as the flow is incremented: a period-doubling cascade, the quasiperiodic route, and type III intermittency. The dynamical behaviour of the system and differing routes to chaos are illustrated by phase-plane portraits, bifurcation diagrams, power spectra, Poincaré sections, and Lyapunov exponent calculations.     

Asymptotic Resonance, Interaction of Modes and Subharmonic Bifurcation

We study the existence of small amplitude oscillations near elliptic equilibria of autonomous systems, which mix different normal modes. The reference problem is the Fermi-Pasta-Ulam β-model: a chain of nonlinear oscillators with nearest-neighborhood interaction. We develop a new bifurcation approach that locates secondary bifurcations from the unimodal primary branches. Two sufficient conditions for bifurcation are given: one involves only the arithmetic properties of the eigenvalues of the linearized system (asymptotic resonance), while the other takes into account the nonlinear character of the interaction between normal modes (nonlinear coupling). Both conditions are checked for the Fermi-Pasta-Ulam problem.     

多孔介质中的双稳热对流

对矩形横截面多孔介质中热对流的复杂分岔行为──二次分岔进行研究．使用Ｌｉａｐｕｎｏｖ－Ｓｃｈｍｉｄｔ约化并充分利用问题本身的对称性，研究了于最低的两个不同临界Ｒａｙｌｅｉｇｈ数处从平凡的静态传热解产生的热对流主分岔解之间的相互作用；揭示了主分岔解的二次分岔并给出了主分岔解及二次分岔解的渐近展开．稳定性分析表明从第二临界Ｒａｙｌｅｉｇｈ数产生的主分岔解经二次分岔后由不稳定变得稳定，从而与由最小临界Ｒａｙｌｅｉｇｈ数产生的主分岔解组成双稳定热对流．文中理论分析可较恰当地解释已有的数值模拟结果．     

Complex dynamics of perfect discrete systems under partial follower forces

Equilibrium points, primary and secondary static bifurcation branches, and periodic orbits with their bifurcations of discrete systems under partial follower forces and no initial imperfections are examined. Equilibrium points are computed by solving sets of simultaneous, non-linear algebraic equations, whilst periodic orbits are determined numerically by solving 2- or 4-dimensional non-linear boundary value problems. A specific application is given with Ziegler's 2-DOF cantilever model. Numerous, complicated static bifurcation paths are computed as well as complicated series of periodic orbit bifurcations of relatively large periods. Numerical simulations indicate that chaotic-like transient motions of the system may appear when a forcing parameter increases above the divergence state. At these forcing parameter values, there co-exist numerous branches of bifurcating periodic orbits of the system; it is conjectured that sensitive dependence on initial conditions due to the large number of co-existing periodic orbits causes the chaotic-like transients observed in the numerical simulations.     

Nonlinear normal modes and their superposition in a two degrees of freedom asymmetric system with cubic nonlinearities

This paper investigates nonlinear normal modes and their superposition in a two degrees of freedom asymmetric system with cubic nonlinearities for all nonsingular conditions, based on the invariant subspace in nonlinear normal modes for the nonlinear equations of motion. The focus of attention is to consider relation between the validity of superposition and the static bifurcation of modal dynamics. The numerical results show that the validity has something to do not only with its local restriction, but also with the static bifurcation of modal dynamics. Project Supported by the National Natural Science Foundation and PSF of China