Singularity analysis of response bifurcation for a coupled pitch–roll ship model with quadratic and cubic nonlinearity

Nonlinear Dynamics - Based on the coupling of roll and pitch motion of ships, a mathematical model with quadratic and cubic nonlinear terms is presented. Primary resonance is discussed by the...     

A hyperbolic Lindstedt–Poincaré method for homoclinic motion of a kind of strongly nonlinear autonomous oscillators

A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.     

Wave–wave interactions in a periodic chain with quadratic nonlinearity

Two waves are studied using perturbation analysis for their interactions in an one-dimensional periodic structure with quadratic nonlinearity. A first-order multiple-scales analysis along with numerical simulations on the full chain are used to understand the interaction of two waves when one is the sub- or super-harmonic of the other. The strength of quadratic nonlinearity affects the rate at which the energy is exchanged between the two waves. Depending on parameters and energy states, the interactions between the waves are periodic or whirling and result in quasi-periodic combined propagating waves with either phase drifts or weakly phase-locking properties. The analysis suggests the possibility of the existence of emergent wave harmonics. Due to quadratic nonlinearity, a very small amplitude subharmonic or superharmonic wave mode can drift in its phase, and then burst out with a larger amplitude as it circumnavigates a separatrix. Depending on the parameters and wave numbers, the amplitude of this emergent wave burst can have varying significance.     

Homoclinic Connections in Strongly Self-Excited Nonlinear Oscillators: The Melnikov Function and the Elliptic Lindstedt–Poincaré Method

A criterion to predict bifurcation of homoclinic orbits instrongly nonlinear self-excited one-degree-of-freedom oscillator

Bifurcation and chaos of a cable–beam coupled system under simultaneous internal and external resonances

The bifurcation and chaos of a cable–beam coupled system under simultaneous internal and external resonances are investigated. The combined effects of the nonlinear term due to the cable’s geometric and coupled behavior between the modes of the beam and the cable are considered. The nonlinear partial-differential equations are derived by the Hamiltonian principle. The Galerkin method is applied to truncate the governing equation into a set of ordinary differential equations. The bifurcation diagrams in three separate loading cases, namely, excitation acting on the cable, on the beam and simultaneously on the beam and cable, are analyzed with changing forcing amplitude. Based on careful numerical simulations, bifurcations and possible chaotic motions are represented to reveal the combined effects of nonlinearities on the dynamics of the beam and the cable when they act as an overall structure.     

Poincaré oscillations and geostrophic adjustment in a rotating paraboloid

Free liquid oscillations (Poincaré oscillations) in a rotating paraboloid are investigated theoretically and experimentally. Within the framework of shallow-water theory, with account for the centrifugal force, expressions for the free oscillation frequencies are obtained and corrections to the frequencies related with the finiteness of the liquid depth are found. It is shown that in the rotating liquid, apart from the wave modes of free oscillations, a stationary vortex mode is also generated, that is, a process of geostrophic adjustment takes place. Solutions of the shallow-water equations which describe the wave dynamics of the adjustment process are presented. In the experiments performed the wave and vortex modes were excited by removing a previously immersed hemisphere from the central part of the paraboloid. Good agreement between theory and experiment was obtained.     

Poincaré and the shortcomings of the Hamilton-Jacobi method for classical or quantized mechanics

A great theorem was proven by H. Poincaré in celestial mechanics. It states that, in the most general problems of mechanics, the total energy of the system is the only well behaved first integral of the system, while other so-called integrals cannot be represented by uniform and convergent series. This very important result can be explained and visualized by comparison with standard methods of discussion, as, for example, the Hamilton-Jacobi procedure. The discussion shows that there are serious limitations to the use of this procedure, which collapses in the most general problems (Poincaré theorem) and can be used only for “almost separated” variables. The Poincaré theorem appears to provide the distinction between determinism in mechanics and statistical mechanics according to Boltzmann. The research presented here done under Contract Nonr 266(56) and was first described in a Quarterly Report dated July 31, 1959.     

Stability of Bragg grating solitons in a semilinear dual-core system with cubic–quintic nonlinearity

The existence and stability of quiescent Bragg grating solitons in a dual-core fiber, where one core contains a Bragg grating with cubic–quintic nonlinearity, and the other is a linear are studied. The model admits two disjoint bandgaps when the relative group velocity in the linear core, c, is zero: one in the upper half and the other in the lower half of the system’s linear spectrum. In the general case (i.e., \(c\ne 0\)), a central gap (which is a genuine gap) is formed, while the lower and upper gaps overlap with one branch of continuous spectrum, and therefore, they are not genuine bandgaps. For quiescent solitons, exact analytical solutions are found in implicit form for \(c=0\). For nonzero c, soliton solutions are obtained numerically. The system supports two disjoint families (referred to as Type 1 and Type 2) of zero-velocity soliton solutions, separated by a border. Both Type 1 and Type 2 soliton solutions exist throughout the upper and lower gaps but not in the central gap. The stability of both soliton families is investigated by means of systematic numerical simulations. It is found that Type 2 solitons are always unstable and are destroyed upon propagation. On the other hand, unstable Type 1 solitons may either decay into radiation or radiate some energy and evolve into a moving Type 1 soliton. Also, in the case of Type 1 solitons, we have identified stable regions in the plane of quintic nonlinearity and frequency. The influence of coupling coefficient and the relative group velocity in the linear core on the stability of solitons are analyzed.     

Spectral Analysis of the Neumann–Poincaré Operator and Characterization of the Stress Concentration in Anti-Plane Elasticity

When holes or hard elastic inclusions are closely located, stress which is the gradient of the solution to the anti-plane elasticity equation can be arbitrarily large as the distance between two inclusions tends to zero. It is important to precisely characterize the blow-up of the gradient of such an equation. In this paper we show that the blow-up of the gradient can be characterized by a singular function defined by the single layer potential of an eigenfunction corresponding to the eigenvalue 1/2 of a Neumann–Poincaré type operator defined on the boundaries of the inclusions. By comparing the singular function with the one corresponding to two disks osculating to the inclusions, we quantitatively characterize the blow-up of the gradient in terms of explicit functions. In electrostatics, our results apply to the electric field, which is the gradient of the solution to the conductivity equation, in the case where perfectly conducting or insulating inclusions are closely located.     

Dynamics of a cylindrical shell with elastic bottom filled with a gas–liquid mixture and subject to two-frequency excitation

The paper reports experimental data on the nonlinear dynamic deformation of the elastic bottom of a cylindrical shell and the formation of bubbles and their clusters under two-frequency excitation     

Bifurcation of periodic solution in a Predator–Prey type of systems with non-monotonic response function and periodic perturbation

A Predator–Prey type of dynamical systems with non-monotonic response function and time-periodic perturbation is considered in this paper. We present a proof for the number of equilibria in the unperturbed system at some parts of the parameter space. The perturbed system is a dynamical system defined by a periodic vector field. We present an alternative proof for a classical result on the period of the periodic solution. By using a numerical continuation method AUTO (Doedel et al., 1986 [9]), we present a bifurcation analysis for periodic solution of the perturbed system where we found fold, cusp and Swallowtail bifurcations.     

Aeroelastic analysis of wind turbines using a tightly coupled CFD–CSD method

This paper presents aeroelastic analyses of wind turbines, using the compressible flow Helicopter Multi-Block (HMB2) solver of Liverpool University, coupled with a Computational Structural Dynamics method. For this study, the MEXICO and NREL Phase VI wind turbines were employed. A static aeroelastic method was first employed for the analysis of the MEXICO blade and the effect of the torsional stiffness was studied at 10, 15 and 24 m/s axial wind speeds. The torsional deformations showed strong dependency on this parameter and the blade region from mid-span to the tip was the most susceptible to aeroelastic effects. The work progressed by studying both the static and dynamic response on the NREL wind turbine, where the nacelle and the tower were considered. Mean deflections between the static and dynamic methods showed consistency and, due to the structural properties of this blade, flapping modes were dominant. The dynamic aeroelastic method enabled an assessment of the effect of flapping on the blade loads, in conjunction with the effect of tower. Aeroelastic effects were found to be secondary for the MEXICO blade, but had a stronger effect on the larger NREL Phase VI blade.     

On occurrence of sudden increase of spiking amplitude via fold limit cycle bifurcation in a modified Van der Pol–Duffing system with slow-varying periodic excitation

Nonlinear Dynamics - The main purpose of the paper is to reveal the mechanism of certain special phenomena in bursting oscillations such as the sudden increase of the spiking amplitude. When...     

Simulation of wave–body interaction: A desingularized method coupled with acceleration potential

A two-dimensional, potential-theory based, fully nonlinear numerical wave tank (NWT) is developed for the simulation of wave–body interaction. In this NWT, the concept of acceleration potential is adopted in addition to the velocity potential. Both potentials are solved using the desingularized boundary integral equation method (DBIEM). By tapping the strength of the DBIEM, a new acceleration-potential solving method is proposed, which turns the originally implicit kinematic boundary condition on the surface of a passively moving body into an explicit one. Unlike the other existing methods such as the mode decomposition method, the indirect method and the iterative method, the present method requires solving of only one boundary value problem to determine the acceleration potential, and hence significantly enhances the computational efficiency. Using this NWT, the nonlinear interaction between a freely floating barge and various incident waves is investigated. It is confirmed that the new acceleration potential solving method outperforms other existing methods, saving at least 45% of the computational time.     

Dynamic analysis of an aeroelastic airfoil with freeplay nonlinearity by precise integration method based on Padé approximation

Nonlinear Dynamics - An improved precise integration method (PIM) incorporated with Padé approximation (PadéPIM) is proposed, and the aeroelastic behavior of an aeroelastic airfoil with...     

Nonlinear dynamics of a non-autonomous network with coupled discrete–continuum oscillators

A network model of a multi-modular floating platform incorporated with a runway structure, viewed as a non-autonomous network with discrete–continuum oscillators, is developed for a general purpose of dynamic analysis. Numerical analysis shows the coupling effect between the two different types of oscillators on various complex dynamics, including sudden leaps, torus motions, beating vibrations, the synergetic effect of phase lock and anti-phase synchronizations. The amplitude death phenomenon, a suppressed weak oscillation state, is studied by using the fundamental solution derived by the averaging method. The parametric domain of the onset of amplitude death is illustrated to show the great significance to the stability design of the floating platform. The effect of the flexural rigidity of the runway on the distribution of amplitude death state is also discussed.     

Stochastic analysis of a predator–prey model with modified Leslie–Gower and Holling type II schemes

Nonlinear Dynamics - This article studies a predator–prey model with modified Leslie–Gower and Holling type II schemes under white-noise disturbances. The sensitivity analysis of the...     

Multiple scales and normal forms in a ring of delay coupled oscillators with application to chaotic Hindmarsh–Rose neurons

We study the appearance and stability of spatiotemporal periodic patterns like phase-locked oscillations, mirror-reflecting waves, standing waves, in-phase or antiphase oscillations, and coexistence of multiple patterns, in a ring of bidirectionally delay coupled oscillators. Hopf bifurcation, Hopf–Hopf bifurcation, and the equivariant Hopf bifurcation are studied in the viewpoint of normal forms obtained by using the method of multiple scales which is a kind of perturbation technique, thus a clear bifurcation scenario is depicted. We find time delay significantly affects the dynamics and induces rich spatiotemporal patterns. With the help of the unfolding system near Hopf–Hopf bifurcation, it is confirmed in some regions two kinds of stable oscillations may coexist. These phenomena are shown for the delay coupled limit cycle oscillators as well as for the delay coupled chaotic Hindmarsh–Rose neurons.     

Bifurcation of a modified railway wheelset model with nonlinear equivalent conicity and wheel–rail force

This paper presents the bifurcation behaviors of a modified railway wheelset model to explore its instability mechanisms of hunting motion. Equivalent conicity data measured from China high-speed railway vehicle are used to modify the wheelset model. Firstly, the relationships between longitudinal stiffness, lateral stiffness, equivalent conicity and critical speed are taken into account by calculating the real parts of the eigenvalues of the Jacobian matrix and Hurwitz criterion for the corresponding linear model. Secondly, measured equivalent conicity data are fitted by a nonlinear function of the lateral displacement rather than are considered as a constant as usual. Nonlinear wheel–rail force function is used to describe the wheel–rail contact force. Based on these modifications, a modified railway wheelset model with nonlinear equivalent conicity and wheel–rail force is set up, and then, some instability mechanisms of China high-speed train vehicle are investigated based on Hopf bifurcation, fold (limit point) bifurcation of cycles, cusp bifurcation of cycles, Neimark–Sacker bifurcation of cycles and 1:1 resonance. In particular, fold bifurcation of cycles can produce a vast effect on the hunting motion of the modified wheelset model. One of the main reasons leading to hunting motion is due to the fold bifurcation structure of cycles, in which stable limit cycles and unstable limit cycles may coincide, and multiple nested limit cycles appear on a side of fold bifurcation curve of cycles. Unstable hunting motion mainly depends on the coexistence of equilibria and limit cycles and their positions; if the most outward limit cycle is stable, then the motion of high-speed vehicle should be safe in a reasonable range. Otherwise, if the initial values are chosen near the most outward unstable limit cycle or the system is perturbed by noises, the high-speed vehicle will take place unstable hunting motion and even lead to serious train derailment events. Therefore, in order to control hunting motions, it may be the easiest way in theory to guarantee the coexistence of the inner stable equilibrium and the most outward stable limit cycle in a wheelset system.

     

Determination of the full-field stress and displacement using photoelasticity and sampling moiré method in a 3D-printed model

The quantitative characterization of the full-field stress and displacement is significant for analyzing the failure and instability of engineering materials. Various optical measurement techniques such as photoelasticity, moiré and digital image correlation methods have been developed to achieve this goal. However, these methods are difficult to incorporate to determine the stress and displacement fields simultaneously because the tested models must contain particles and grating for displacement measurement; however, these elements will disturb the light passing through the tested models using photoelasticity. In this study, by combining photoelasticity and the sampling moiré method, we developed a method to determine the stress and displacement fields simultaneously in a three-dimensional (3D)-printed photoelastic model with orthogonal grating. Then, the full-field stress was determined by analyzing 10 photoelastic patterns, and the displacement fields were calculated using the sampling moiré method. The results indicate that the developed method can simultaneously determine the stress and displacement fields.