Three-Period Quasi-Periodic Solutions in the Self-Excited Quasi-Periodic Mathieu Oscillator

In this work, we investigate three-period quasi-periodic (QP) oscillations in the vicinity of 2:2:1 resonance in a self-excited QP Mathieu equation using perturbation method. Two successive averaging are performed to reduce the original QP equation to an autonomous amplitude and phase system describing the modulation of the slow flow dynamic. Approximation of three-period QP solution is obtained via the study of limit cycle of the reduced autonomous system. The efficiency of the method is illustrated by comparison between analytical approximations and numerical integration. The double reduction procedure, applied in previous works to construct two-period QP solution, can be implemented to approximate excplicit analytical three-period QP solutions.     

Asymptotic solutions for a damped non-linear quasi-periodic Mathieu equation

Quasi-periodic (QP) solutions of a weakly damped non-linear QP Mathieu equation are investigated near a double primary parametric resonance. A double multiple scales method is applied to reduce the original QP oscillator to an autonomous system performing two successive reduction. The problem for approximating QP solutions of the original system is then transformed to the study of stationary regimes of the induced autonomous system. Explicit analytical approximations to QP oscillations are obtained and comparisons to numerical integration of the original QP oscillator are provided.     

Quasi-Periodic Oscillations,Chaos and Suppression of Chaos in a Nonlinear Oscillator Driven by Parametric and External Excitations

An analysis is given of the dynamic of a one-degree-of-freedom oscillator with quadratic and cubic nonlinearities subjected to parametric and external excitations having incommensurate frequencies. A new method is given for constructing an asymptotic expansion of the quasi-periodic solutions. The generalized averaging method is first applied to reduce the original quasi-periodically driven system to a periodically driven one. This method can be viewed as an adaptation to quasi-periodic systems of the technique developed by Bogolioubov and Mitropolsky for periodically driven ones. To approximate the periodic solutions of the reduced periodically driven system, corresponding to the quasi-periodic solution of the original one, multiple-scale perturbation is applied in a second step. These periodic solutions are obtained by determining the steady-state response of the resulting autonomous amplitude-phase differential system. To study the onset of the chaotic dynamic of the original system, the Melnikov method is applied to the reduced periodically driven one. We also investigate the possibility of achieving a suitable system for the control of chaos by introducing a third harmonic parametric component into the cubic term of the system.     

Generation and evolution of oblique solitary waves in supercritical flows

It is considered that a thin strut sits in a supercritical shallow water flow sheet over a homogeneous or very mildly varying topography. This stationary 3-D problem can be reduced from a Boussinesq-type equation into a KdV equation with a forcing term due to uneven topography, in which the transverse coordinate Y plays a same role as the time in original KdV equation. As the first example a multi-soliton wave pattern is shown by means of N-soliton solution. The second example deals with the generation of solitary wave-train by a wedge-shaped strut on an even bottom. Whitham's average method is applied to show that the shock wave jump at the wedge vertex develops to a cnoidal wave train and eventually to a solitary wavetrain. The third example is the evolution of a single oblique soliton over a periodically varying topography. The adiabatic perturbation result due to Karpman & Maslov (1978) is applied. Two coupled ordinary differential equations with periodic disturbance are obtained for the soliton amplitude and phase. Numerical solutions of these equations show chaotic patterns of this perturbed soliton.     

Averaging and periodic solutions in the plane and parametrically excited pendulum

We first approximate the solutions of the nonautonomous oscillating suspension point pendulum equation by the solutions of a second order autonomous differential equation. Using the strict monotonicity of the periodic solutions of the approximating equation, we prove the existence of a large number of subharmonic periodic solutions of the plane pendulum when its point of suspension is excited parametrically.     

Delay equations with fluctuating delay related to the regenerative chatter

The mathematical models representing machine tool chatter dynamics have been cast as differential equations with delay. In this paper, non-linear delay differential equations with periodic delays which model the machine tool chatter with continuously modulated spindle speed are studied. The explicit time-dependent delay terms, due to spindle speed modulation, are replaced by state-dependent delay terms by augmenting the original equations. The augmented system of equations is autonomous and has two pairs of pure imaginary eigenvalues without resonance. The reduced bifurcation equation is obtained by making use of Lyapunov-Schmidt Reduction method. By using the reduced bifurcation equations, the periodic solutions are determined to analyze the tool motion. Analytical results show both modest increase of stability and existence of periodic solutions near the new stability boundary.     

Stationary response of Duffing oscillator with hardening stiffness and fractional derivative

The stationary response of Duffing oscillator with hardening stiffness and fractional derivative under Gaussian white noise excitation is studied. First, the term associated with fractional derivative is separated into the equivalent quasi-linear dissipative force and quasi-linear restoring force by using the generalized harmonic balance technique, and the original system is replaced by an equivalent nonlinear stochastic system without fractional derivative. Then, the stochastic averaging method of energy envelope is applied to the equivalent nonlinear stochastic system to yield the averaged Itô equation of energy envelope, from which the corresponding Fokker–Planck–Kolmogorov (FPK) equation is established and solved to obtain the stationary probability densities of the energy envelope and the amplitude envelope. The accuracy of the analytical results is validated by those from the Monte Carlo simulation of original system.     

Bifurcation analysis of reduced rotor model based on nonlinear transient POD method

The singularity theory is applied to study the bifurcation behaviors of a reduced rotor model obtained by nonlinear transient POD method in this paper. A six degrees of freedom (DOFs) rotor model with cubically nonlinear stiffness supporting at both ends is established by the Newton's second law. The nonlinear transient POD method is used to reduce a six-DOFs model to a one-DOF one. The reduced model reserves the dynamical characteristics and occupies most POM energy of the original one. The singularity of the reduced system is analyzed, which replaces the original system. The bifurcation equation of the reduced model indicates that it is a high co-dimension bifurcation problem with co-dimension 6, and the universal unfolding (UN) is provided. The transient sets of six unfolding parameters, the bifurcation diagrams between the bifurcation parameter and the state variable are plotted. The results obtained in this paper present a new kind of method to study the UN theory of multi-DOFs rotor system.     

Response and stability analysis of piecewise-linear oscillators under multi-forcing frequencies

Steady-state solutions of a piecewise-linear oscillator under multi-forcing frequencies are obtained using the fixed point algorithm (FPA). Stability analysis is also performed using the same technique. For the periodic solutions of a piecewise-linear oscillator with single forcing frequency, the harmonic balance method (HBM) is also used along with the FPA. Although both FPA and HBM generate accurate solutions, it is observed that the HBM failed to converge to solutions in the superharmonic range of the forcing frequency.The fixed point algorithm was also applied to the oscillator under multifrequency excitation. The algorithm proved to be very effective in obtaining torus solutions and in locating corresponding bifurcation thresholds. A piecewise-linear oscillator model of an offshore articulated loading platform (ALP) subjected to two incommensurate wave frequencies is found to exhibit chaotic behavior. A second order Poincaré mapping technique reveals the hidden fractal-like nature of the resulting chaotic response. A parametric study is performed for the response of the ALP.     

Stationary points iteration method for periodic solution to nonlinear system

The value method which is used to obtain the periodic solution to nonlinear system is mentioned in this article. Different point reflection is defined in the nonlinear autonomous and nonautonomous system firstly and then that linear reflection obtained from the inserting value of nonlinear reflection is asymptotic to original nonlinear reflection. The stationary points obtained by linear reflection are regarded as the asymptotic solution of the stationary points of original system. If this asymptotic solution of the stationary points is not satisfactorily accurate it can be used as the initial point of the next reflection. In addition, a corresponding method of researching the stability of periodic solution is put forward in this article.     

Application of first-order canonical perturbation method with dissipative Hori-like kernel

Lie-Hori canonical perturbation theory provides asymptotic solutions for conservative Hamiltonian systems. This restriction prevents the canonical method from being applied directly to dissipative mechanical systems. There are, however, two main alternatives to overcome this difficulty, enabling the application of canonical perturbation methods. The first one consists in constructing a time-dependent Hamiltonian, through a generating function, related to the energy dissipation pattern of the system. The second embeds the original phase space into a double dimensional one where the dynamics of the system can be formulated in a Hamiltonian way. In this paper, a modified Lie-Hori method that avoid the disadvantages of the former approaches is proposed. Namely, it is not necessary to find out a time-dependent generating function, nor doubling the number of the canonical variables of the original problem. The new algorithm provides first order analytical solutions for a certain set of dissipative non-linear dynamical systems. It is based on a suitable modification of the Hori kernel in the double-dimensional embedding phase space, allowing the inclusion of the dissipative (or generalized) forces. By means of this redefined auxiliary system, the path-integrals of the method can be performed in a domain of the phase space with the same dimensionality as the original problem.     

Exact solutions to magnetogasdynamics using Lie point symmetries

In the present work, we find some exact solutions to the first order quasilinear hyperbolic system of partial differential equations (PDEs), governing the one dimensional unsteady flow of inviscid and perfectly conducting compressible fluid, subjected to a transverse magnetic field. For this, Lie group analysis is used to identify a finite number of generators that leave the given system of PDEs invariant. Out of these generators, two commuting generators are constructed involving some arbitrary constants. With the help of canonical variables associated with these two generators, the assigned system of PDEs is reduced to an autonomous system whose simple solutions provide nontrivial solutions of the original system. Using this exact solution, we discuss the evolutionary behavior of weak discontinuities.     

Responses of a non-linearly coupled pitch-roll ship model under harmonic excitation

Stationary responses of nonlinearly coupled pitch and roll ship modes are studied using a modified averaging method, along with two second order multiple time scale (MTS) methods for comparison. Stability of the solutions is also studied. In the case of harmonic excitation all the three methods give fairly accurate results to the original problem but the averaging method is the most efficient. Analytic solutions are obtained from the averaged equations, which can be used to predict stationary responses both for small and for large excitations. From the averaging method several qualitatively different phenomena which cannot be addressed by the first order theory have been obtained: (i) the saturation phenomenon is lost, (ii) the bifurcation points are altered and (iii) a drift term is present which, although small, appears to have a significant effect on the accuracy of the solutions.     

Exact stationary probability density functions for non-linear systems under Poisson white noise excitation

The paper presents exact stationary probability density functions for systems under Poisson white noise excitation. Two different solution methods are outlined. In the first one, a class of non-linear systems is determined whose state vector is a memoryless transformation of the state vector of a linear system. The second method considers the generalized Fokker-Planck (Kolmogorov-forward) equation. Non-linear system functions are identified such that the stationary solution of the system admits a prescribed stationary probability density function. Both methods make use of the stochastic integro-differential equations approach. This approach seems to have some computational advantages for the determination of exact stationary probability density functions when compared to the stochastic differential equations approach.     

Perturbation methods for nonlinear autonomous discrete-time dynamical systems

Two perturbation methods for nonlinear autonomous discrete-time dynamical systems are presented. They generalize the classical Lindstedt-Poincaré and multiple scale perturbation methods that are valid for continuous-time systems. The Lindstedt-Poincaré method allows determination of the periodic or almost-periodic orbits of the nonlinear system (limit cycles), while the multiple scale method also permits analysis of the transient state and the stability of the limit cycles. An application to the discrete Van der Pol equation is also presented, for which the asymptotic solution is shown to be in excellent agreement with the exact (numerical) solution. It is demonstrated that, when the sampling step tends to zero the asymptotic transient and steady-state discrete-time solutions correctly tend to the asymptotic continuous-time solutions.     

Effect of External Excitations on a Nonlinear System with Time Delay

The trivial equilibrium of a two-degree-of-freedom autonomous system may become unstable via a Hopf bifurcation of multiplicity two and give rise to oscillatory bifurcating solutions, due to presence of a time delay in the linear and nonlinear terms. The effect of external excitations on the dynamic behaviour of the corresponding non-autonomous system, after the Hopf bifurcation, is investigated based on the behaviour of solutions to the four-dimensional system of ordinary differential equations. The interaction between the Hopf bifurcating solutions and the high level excitations may induce a non-resonant or secondary resonance response, depending on the ratio of the frequency of bifurcating periodic motion to the frequency of external excitation. The first-order approximate periodic solutions for the non-resonant and super-harmonic resonance response are found to be in good agreement with those obtained by direct numerical integration of the delay differential equation. It is found that the non-resonant response may be either periodic or quasi-periodic. It is shown that the super-harmonic resonance response may exhibit periodic and quasi-periodic motions as well as a co-existence of two or three stable motions.     

Second order hyperbolic equations with small non-linearities in the case of internal resonance

The monofrequent solutions of certain autonomous second order hyperbolic differential equations with weak non-linearities are found in the case when some of the natural frequencies of the generating equation are in integral ratio. The approach use is a development of the KrylovBogoliubov-Mitropolskii method. The solution found is applied to the case of the longitudinal vibrations of a rod for which the stress-strain relation contains a small non-linear term and to the case of the vibrations of a rod with small inhomogeneities of density and elastic modulus and with small damping.     

A staggered grid‐based explicit jump immersed interface method for two‐dimensional Stokes flows

In this paper the explicit jump immersed interface method (EJIIM) is applied to stationary Stokes flows. The boundary value problem in a general, non‐grid aligned domain is reduced by the EJIIM to a sequence of problems in a rectangular domain, where staggered grid‐based finite differences for velocity and pressure variables are used. Each of these subproblems is solved by the fast Stokes solver, consisting of the pressure equation (known also as conjugate gradient Uzawa) method and a fast Fourier transform‐based Poisson solver. This results in an effective algorithm with second‐order convergence for the velocity and first order for the pressure. In contrast to the earlier versions of the EJIIM, the Dirichlét boundary value problem is solved very efficiently also in the case when the computational domain is not simply connected. Copyright © 2007 John Wiley & Sons, Ltd.     

On Double Crater-Like Probability Density Functions of a Duffing Oscillator Subjected to Harmonic and Stochastic Excitation

It is a well-known phenomenon of the Duffing oscillator under harmonic excitation,that there is a frequency range, where two stable and one unstable stationarysolution coexist. If the Duffing oscillator is harmonically excited in thisfrequency range and additionally excited, e.g. by white noise, a double crater-likeprobability density function can be observed, if the noise intensity is smallcompared to the harmonic excitation. The aim of this paper is to calculate thisprobability density function approximately using perturbation techniques. Thestationary solutions in the deterministic case are calculated using theperturbation technique for the resonance case. In a second step, the probabilitydensity function of the perturbation of each of those stationary solutions iscalculated using the perturbation technique for the nonresonance case. This resultsin two crater-like probability density functions which are superimposed by usingthe probability of realization of each of the stationary solutions in thedeterministic case. The probability is calculated using numerical integration orthe method of slowly changing phase and amplitude. Finally, probability densityfunctions obtained in this manner are compared to Monte Carlo simulations.     

Dynamic stability and response of fluttered beams subjected to random follower forces

This paper presents a study of non-linear response of a fluttered, cantilevered beam subjected to a random follower force at the free end. The random follower force is characterized as the sum of a post-critical static force and a stationary process with a zero mean. First, the Ritz-Galerkin method is applied to yield a set of discretized system equations. The system equations are then partially uncoupled by a special modal analysis based on normal modes of the corresponding linear, autonomous system at the onset of fluttering. Next, the stochastic averaging method is utilized to get Ito's differential equation governing the amplitude of the fluttered mode. Finally, the probability density function for the amplitude of the fluttered mode is obtained by solving the FPK equation. Numerical results show that the probability density function for the amplitude of the fluttered mode is determined by the sample behavior of the beam near the trivial equilibrium configuration.