Bifurcation and chaos of a simple walking model driven by a rhythmic signal

The walk of animals is achieved by the interaction between the dynamics of their mechanical system and the central pattern generator (CPG). In this paper, we analyze dynamic properties of a simple walking model of a biped robot driven by a rhythmic signal from an oscillator. In particular, we examine the long-term global behavior and the bifurcation of the motion that leads to chaotic motion, depending on the model parameter values. The simple model consists of a hip and two legs connected at the hip through a rotational joint. The joint is driven by a rhythmic signal from an oscillator, which is an open loop. In order to analyze the bifurcation, we first obtained approximate solutions of the walking motion and then constructed discrete dynamics using the Poincaré map. As a result, we found that consecutive period-doubling bifurcations occur as the model parameter values change, and that the walking motion leads to chaotic motion over the critical value of the model parameters. Moreover, we approximately obtained the period-doubling solutions and the critical value by employing a Newton-Raphson method. Our analytical results were verified by the numerical simulations.     

Period-doubling route to chaos for a swirling flow in an open cylindrical container with a rotating disk

Conclusion The period-doubling route to chaos for a swirling vortex flow in an open cylindrical container, at an aspect ratio of 2, driven by a rotating bottom disk was recognized by using laser-Doppler velocimetry. The onset of periodic motion for the flow was found when Re was in the range between 1850 and 1900. The flow was subharmonically bifurcated into a double-period motion when Re was about 2150. When the Reynolds number was in the range from 2300 to 2400, the flow bifurcated again through the period-doubling mode. When the Reynolds number was further increased, the flow eventually showed chaotic motion. The existence of a free surface promotes the onset of periodicity, and the difference of the critical Reynolds number with and without a free surface was estimated to be about 600.This work was supported by the National Science Council of the Republic of China under grant No. NSC-82-0410-E-002-191     

Stability, bifurcation and chaos of a delayed oscillator with negative damping and delayed feedback control

This paper presents a detailed analysis on the dynamics of a delayed oscillator with negative damping and delayed feedback control. Firstly, a linear stability analysis for the trivial equilibrium is given. Then, the direction of Hopf bifurcation and stability of periodic solutions bifurcating from trivial equilibrium are determined by using the normal form theory and center manifold theorem. It shows that with properly chosen delay and gain in the delayed feedback path, this controlled delayed system may have stable equilibrium, or periodic solutions, or quasi-periodic solutions, or coexisting stable solutions. In addition, the controlled system may exhibit period-doubling bifurcation which eventually leads to chaos. Finally, some new interesting phenomena, such as the coexistence of periodic orbits and chaotic attractors, have been observed. The results indicate that delayed feedback control can make systems with state delay produce more complicated dynamics.     

Alternating chaos versus synchronized vibrations of interacting plate with beams

We study networks of coupled oscillators governed by ODEs and yielded by physically validated sets of a few PDEs governing dynamics of structural members (plate and beams), chaos and phase synchronization and contact/no-contact non-linear dynamics of structural members coupled via boundary conditions. We have detected, illustrated and discussed a few novel kinds of hybrid states of the studied plate-beam(s) contact/no-contact interactions as well as novel scenarios of transition into chaos exhibited by the interplay of continuous objects. Classical (time histories, phase portraits, Poincaré maps, FFT, Lyapunov exponents) and non-classical (2D Morlet wavelets) approaches are used while monitoring non-linear dynamics of the interacting spatial structural members. Our results include examples from structural mechanics and the studied objects are modelled by validated mechanical hypotheses and assumptions. Novel non-linear phenomena including switching to different vibration regimes and phase chaotic synchronization are illustrated and discussed.     

Periodicity and chaos in a flexible crank-rocker mechanism

In this paper we introduce a crank-rocker mechanism at which the rocker is flexible. Using Hamilton’s principle we obtain the governing equations of motion for the elastic mode of the rocker. By applying the Bubnov–Galerkin global averaging method, we reduce the governing equations of motion to an ordinary differential equation which is Duffing’s oscillator with time varying coefficients. Through the application of Banach’s fixed-point theorem we predict the periodic solutions. Then we study the geometrical features of the motion near the 1 : 1, 1 : 2 and 2 : 1 commensurabilities. It is also shown that homoclinic and heteroclinic orbits can exist for the system.     

On the response of a preloaded mechanical oscillator with a clearance: Period-doubling and chaos

The dynamic behavior of a harmonically excited, preloaded mechanical oscillator with dead-zone nonlinearity is described quantiatively. The governing strongly nonlinear differential equation is solved numerically. Damping coefficient-force ratio maps for two different values of the excitation frequency have been formed and the boundaries of the regions of different motion types are determined. The results have been compared with the results of the forced Duffing's equation available in the literature in order to identify the differences between cubic and dead-zone nonlinearities. Period-doubling bifurcations, which take place with a change of any of the system parameters, have been found to be the most common route to chaos. Such bifurcations follow the scaling rule of Feigenbaum. b half length of the clearance.     

Chaotic motions of a rigid rotor in short journal bearings

In the present paper the conditions that give rise to chaotic motions in a rigid rotor on short journal bearings are investigated and determined. A suitable symmetry was given to the rotor, to the supporting system, to the acting system of forces and to the system of initial conditions, in order to restrict the motions of the rotor to translatory whirl. For an assigned distance between the supports, the ratio between the transverse and the polar mass moments of the rotor was selected conveniently small, with the aim of avoiding conical instability. Since the theoretical analysis of a system's chaotic motions can only be carried out by means of numerical investigation, the procedure here adopted by the authors consists of numerical integration of the rotor's equations of motion, with trial and error regarding the three parameters that characterise the theoretical model of the system: m, the half non-dimensional mass of the rotor, , the modified Sommerfeld number relating to the lubricated bearings, and , the dimensionless value of rotor unbalance. In the rotor's equations of motion, the forces due to the lubricating film are written under the assumption of isothermal and laminar flow in short bearings. The number of numerical trials needed to find the system's chaotic responses has been greatly reduced by recognition of the fact that chaotic motions become possible when the value of the dimensionless static eccentricity % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbnL2yY9% 2CVzgDGmvyUnhitvMCPzgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqe% fqvATv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0d% Xdh9vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9% pgea0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaaca% qabeaadaabauaaaOqaaiabew7aLnaaBaaaleaacaWGZbaabeaaaaa!4046!\[\varepsilon _s \] is greater than 0.4. In these conditions, non-periodic motions can be obtained even when rotor unbalance values are not particularly high (=0.05), whereas higher values (>0.4) make the rotor motion periodic and synchronous with the driving rotation. The present investigation has also identified the route that leads an assigned rotor to chaos when its angular speed is varied with prefixed values of the dimensionless unbalance . The theoretical results obtained have then been compared with experimental data. Both the theoretical and the experimental data have pointed out that in the circumstances investigated chaotic motions deserve more attention, from a technical point of view, than is normally ascribed to behaviours of this sort. This is mainly because such behaviours are usually considered of scarce practical significance owing to the typically bounded nature of chaotic evolution. The present analysis has shown that when the rotor exhibits chaotic motions, the centres of the journals describe orbits that alternate between small and large in an unpredictable and disordered manner. In these conditions the thickness of the lubricating film can assume values that are extremely low and such as to compromise the efficiency of the bearings, whereas the rotor is affected by inertia forces that are so high as to determine severe vibrations of the supports.Nomenclature C radial clearance of bearing (m) - D diameter of bearing (m) - e dimensional eccentricity of journal (m) - e _s value of e corresponding to the static position of the journal - E dimensional static unbalance of rotor (m) - f _x, f _y =F _x/(P), F _y/(P): non-dimensional components of fluid film force - F _x, F _y dimensional components of fluid film force (N) - g acceleration of gravity (m/s²) - L axial length of bearing (m) - m % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbnL2yY9% 2CVzgDGmvyUnhitvMCPzgarmWu51MyVXgaruWqVvNCPvMCG4uz3bqe% fqvATv2CG4uz3bIuV1wyUbqee0evGueE0jxyaibaieYlf9irVeeu0d% Xdh9vqqj-hEeeu0xXdbba9frFf0-OqFfea0dXdd9vqaq-JfrVkFHe9% pgea0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaaca% qabeaadaabauaaaOqaaiabg2da9maalaaabaGaeqyYdC3aaWbaaSqa% beaacaaIYaaaaaGcbaGaeqyYdC3aa0baaSqaaGabciaa-bdaaeaaca% WFYaaaaaaakiabg2da9maalaaabaGaeqyYdC3aaWbaaSqabeaacaaI% YaaaaOGaam4qaaqaaiabeo8aZjaadEgaaaaaaa!4C14!\[ = \frac{{\omega ^2 }}{{\omega _0^2 }} = \frac{{\omega ^2 C}}{{\sigma g}}\]: half non-dimensional mass of rotor - M half mass of rotor (kg) - n angular speed of rotor (in r.p.m.=60/2) - t time     

Bifurcation and chaos of a flag in an inviscid flow

A two-dimensional model is developed to study the flutter instability of a flag immersed in an inviscid flow. Two dimensionless parameters governing the system are the structure-to-fluid mass ratio M^⁎ and the dimensionless incoming flow velocity U^⁎. A transition from a static steady state to a chaotic state is investigated at a fixed M^⁎=1 with increasing U^⁎. Five single-frequency periodic flapping states are identified along the route, including four symmetrical oscillation states and one asymmetrical oscillation state. For the symmetrical states, the oscillation frequency increases with the increase of U^⁎, and the drag force on the flag changes linearly with the Strouhal number. Chaotic states are observed when U^⁎ is relatively large. Three chaotic windows are observed along the route. In addition, the system transitions from one periodic state to another through either period-doubling bifurcations or quasi-periodic bifurcations, and it transitions from a periodic state to a chaotic state through quasi-periodic bifurcations.     

The global bifurcation and chaos for the coupling between longitudinal and transverse vibration of a thin elastic plate in large overall motion

This paper presents the global bifurcation and chaotic behavior for the coupling of longitudinal and transverse vibration of a thin elastic plate in large overall motion. First the parametric equations of the homoclinic orbits of such a system is obtained. Then, by using the Melnikov method and digital computer simulation. the behavior of bifurcation and chaos of this vibration system is investigated in the cases of different resonances. The obvious difference between the transverse vibration and the coupling of transverse and longitudinal vibration is also shown.The project supported by the National Natural Science Foundation of China.     

Dynamical models and the onset of chaos in space debris

The increasing threat raised by space debris led to the development of different mathematical models and approaches to investigate the dynamics of small particles orbiting around the Earth. The choice of such models and methods strongly depend on the altitude of the objects above Earth's surface, since the strength of the different forces acting on an Earth orbiting object (geopotential, atmospheric drag, lunar and solar attractions, solar radiation pressure, etc.) varies with the altitude of the debris.In this review, our focus is on presenting different analytical and numerical approaches employed in modern studies of the space debris problem. We start by considering a model including the geopotential, solar and lunar gravitational forces and the solar radiation pressure. We summarize the equations of motion using different formalisms: Cartesian coordinates, Hamiltonian formulation using Delaunay and epicyclic variables, Milankovitch elements. Some of these methods lead in a straightforward way to the analysis of resonant motions. In particular, we review results found recently about the dynamics near tesseral, secular and semi-secular resonances.As an application of the above methods, we proceed to analyze a timely subject, namely the possible causes for the onset of chaos in space debris dynamics. Precisely, we discuss the phenomenon of overlapping of resonances, the effect of a large area-to-mass ratio, the influence of lunisolar secular resonances.We conclude with a short discussion about the effect of the dissipation due to the atmospheric drag and we provide a list of minor effects, which could influence the dynamics of space debris.     

On a New Route to Chaos in Railway Dynamics

Cooperrider's mathematical model of a railway bogie running on a straight track has been thoroughly investigated due to its interesting nonlinear dynamics (see True [1] for a survey). In this article a detailed numerical investigation is made of the dynamics in a speed range, where many solutions exist, but only a couple of which are stable. One of them is a chaotic attractor.Cooperrider's bogie model is described in Section 2, and in Section 3 we explain the method of numerical investigation. In Section 4 the results are shown. The main result is that the chaotic attractor is created through a period-doubling cascade of the secondary period in an asymptotically stable quasiperiodic oscillation at decreasing speed. Several quasiperiodic windows were found in the chaotic motion.This route to chaos was first described by Franceschini [9], who discovered it in a seven-mode truncation of the plane incompressible Navier–Stokes equations. The problem investigated by Franceschini is a smooth dynamical system in contrast to the dynamics of the Cooperrider truck model. The forcing in the Cooperrider model includes a component, which has the form of a very stiff linear spring with a dead band simulating an elastic impact. The dynamics of the Cooperrider truck is therefore non-smooth.The quasiperiodic oscillation is created in a supercritical Neimark bifurcation at higher speeds from an asymmetric unstable periodic oscillation, which gains stability in the bifurcation. The bifurcating quasiperiodic solution is initially unstable, but it gains stability in a saddle-node bifurcation when the branch turns back toward lower speeds.The chaotic attractor disappears abruptly in what is conjectured to be a blue sky catastrophe, when the speed decreases further.     

Delamination of a high-temperature sandwich plate

This paper describes the development of a method for determining the fracture toughness of the core/faceplate bond in high-temperature sandwich plates. The tensile deformation behavior of a sandwich element was also determined. The results from the latter experiment were used in a beam on elastic foundation analysis of the fracture specimen. The faceplate/core toughness was determined at 23 and 180°C. The room temperature toughness was slightly higher and, in both cases, the toughness decreased with crack length. The higher toughness was associated with a greater degree of interlaminar failure in the faceplates, as opposed to core-pullout.     

Low speed flutter and limit cycle oscillations of a two-degree-of-freedom flat plate in a wind tunnel

This paper explores the dynamical response of a two-degree-of-freedom flat plate undergoing classical coupled-mode flutter in a wind tunnel. Tests are performed at low Reynolds number (Re~2.5×10⁴), using an aeroelastic set-up that enables high amplitude pitch–plunge motion. Starting from rest and increasing the flow velocity, an unstable behaviour is first observed at the merging of frequencies: after a transient growth period the system enters a low amplitude limit-cycle oscillation regime with slowly varying amplitude. For higher velocity the system transitions to higher-amplitude and stable limit cycle oscillations (LCO) with amplitude increasing with the flow velocity. Decreasing the velocity from this upper LCO branch the system remains in stable self-sustained oscillations down to 85% of the critical velocity. Starting from rest, the system can also move toward a stable LCO regime if a significant perturbation is imposed. Those results show that both the flutter boundary and post-critical behaviour are affected by nonlinear mechanisms. They also suggest that nonlinear aerodynamic effects play a significant role.     

Chaotic attitude motion of a magnetic rigid spacecraft and its control

This paper deals with chaotic attitude motion of a magnetic rigid spacecraft with internal damping in a circular orbit near the equatorial plane of the earth. The dynamical model of the problem is established. The Melnikov analysis is carried out to prove the existence of a complicated non-wandering Cantor set. The dynamical behaviors are numerically investigated by means of time history. Poincare map, power spectrum and Lyapunov exponents. Numerical simulations indicate that the onset of chaos is characterized by the intermittency as the increase of the torque of the magnetic forces and decrease of the damping. The input-output feedback linearization method is applied to control chaotic attitude motions to the given fixed point and periodic motion.     

Solitary waves and chaos in nonlinear visco-elastic rod

The dynamics behavior of a nonlinear visco-elastic rod subjected to axially periodic load is investigated theoretically and numerically. The weak longitudinal periodic load is distributed uniformly along the rod. Firstly, equation of motion of the rod is derived. Utilizing perturbation technique, we acquire Kdv type equation describing strain wave in the rod. By use traveling wave method, the elliptic cosine wave solution and the solitary wave solution in the rod are provided. Then, Melnikov method is applied to analyze the dynamic behaviour of the rod qualitatively. The explicit conditions for the onset of chaotic dynamics are yielded. With the help of the Poincare map method, phase trajectory and time-displacement history diagrams, the theoretical results obtained are checked.     

Hysteresis modelling and chaos prediction in one- and two-DOF hysteretic models

In the present work hysteresis is simulated by means of internal variables. The analytical models of different types of hysteresis loops allow the reproduction of major and minor loops and provide a high degree of correspondence with experimental data. In models of this type adding an external periodic excitation or increasing the number of dimensions can lead to the occurrence of chaotic behaviour. Using an effective algorithm based on numerical analysis of the wandering trajectories [1–7], the evolution of the chaotic behaviour regions of oscillators with hysteresis is presented in various parametric planes. The substantial influence of a hysteretic dissipation value on the form and location of these regions, as well as the restraining and generating effects of hysteretic dissipation on the occurrence of chaos, are ascertained. Conditions for pinched hysteresis are defined. Furthermore, autonomous coupled hysteretic oscillators under sliding friction are investigated. Conditions for the occurrence of chaotic behaviour in a two-degree-of-freedom (two-DOF) hysteretic system are found in the plane of maximal static friction forces of both oscillators versus belt velocity.     

Regular oscillations or chaos in a fractional order system with any effective dimension

This paper introduces a fractional order system which can generate regular oscillations or create chaos. It shows that this system is capable to create regular or nonregular oscillations under suitable conditions. These necessary conditions are achieved by violation of the no-chaos criteria. The effective dimension of the proposed system can be chosen any order less than three. Therefore, this system is a good example for limit cycle or chaos generation via fractional-order systems with low orders. Numerical simulations illustrate behavior of the proposed system in different situations.     

磁场中旋转运动圆板的分叉与混沌研究

应用数值模拟方法研究磁场中旋转运动圆板的分叉与混沌问题。首先,基于薄板理论和麦克斯韦电磁场方程组,给出了动能、应变势能、外力虚功以及电磁力的表达式,再利用哈密顿原理,得到磁场中旋转运动圆板横向振动的非轴对称非线性磁弹性振动微分方程组。其次,采用贝塞尔函数作为圆板的振型函数进行伽辽金积分,得到了轴对称情况下横向振动的常微分方程组表达式。最后,针对主共振,取周边夹支边界条件的圆板作为算例,得到了当振型函数取一阶时,将磁感应强度、外激励振幅和激励频率作为控制参数的分叉图及庞加莱映射图等计算结果,并讨论了分叉参数对系统的分叉与混沌的影响。数值计算结果表明,这些控制参数的变化影响系统稳定性,在分叉参数逐渐变化的过程中,系统经历从混沌到多倍周期运动再到混沌的往复过程。