非自旋航天器混沌姿态运动及其参数开闭环控制

研究万有引力场中受大气阻力且存在结构内阻尼的非自旋航天器在椭圆轨道上平面天平动的混沌及其参数开闭环控制问题．在建立数学模型的基础上确定出现混沌的必要条件并数值验证混沌的存在性，提出非线性振动系统混沌运动的参数开闭环控制并应用于控制航天器的混沌姿态运动．     

Subharmonic Melnikov method of six-dimensional nonlinear systems and application to a laminated composite piezoelectric rectangular plate

In this paper, the bifurcations of subharmonic orbits are investigated for six-dimensional non-autonomous nonlinear systems using the improved subharmonic Melnikov method. The unperturbed system is composed of three independent planar Hamiltonian systems such that the unperturbed system has a family of periodic orbits. The key problem at hand is the determination of the sufficient conditions on some of the periodic orbits for the unperturbed system to generate the subharmonic orbits after the periodic perturbations. Using the periodic transformations and the Poincaré map, an improved subharmonic Melnikov method is presented. Two theorems are obtained and can be used to analyze the subharmonic dynamic responses of six-dimensional non-autonomous nonlinear systems. The subharmonic Melnikov method is directly utilized to investigate the subharmonic orbits of the six-dimensional non-autonomous nonlinear system for a laminated composite piezoelectric rectangular plate. Using the subharmonic Melnikov method, the bifurcation function of the subharmonic orbit is obtained. Numerical simulations are used to verify the analytical predictions. The results of the numerical simulation also indicate the existence of the subharmonic orbits for the laminated composite piezoelectric rectangular plate.     

Periodic and Homoclinic Motions in Forced,Coupled Oscillators

We study periodic and homoclinic motions in periodically forced, weakly coupled oscillators with a form of perturbations of two independent planar Hamiltonian systems. First, we extend the subharmonic Melnikov method, and give existence, stability and bifurcation theorems for periodic orbits. Second, we directly apply or modify a version of the homoclinic Melnikov method for orbits homoclinic to two types of periodic orbits. The first type of periodic orbit results from persistence of the unperturbed hyperbolic periodic orbit, and the second type is born out of resonances in the unperturbed invariant manifolds. So we see that some different types of homoclinic motions occur. The relationship between the subharmonic and homoclinic Melnikov theories is also discussed. We apply these theories to the weakly coupled Duffing oscillators.     

A method for classifying orbits near asteroids

A method for classifying orbits near asteroids under a polyhedral gravitational field is presented, and may serve as a valuable reference for spacecraft orbit design for asteroid exploration. The orbital dynamics near aster- oids are very complex. According to the variation in orbit characteristics after being affected by gravitational perturbation during the periapsis passage, orbits near an as- teroid can be classified into 9 categories： （1） surrounding- to-surrounding, （2） surrounding-to-surface, （3） surrounding- to-infinity, （4） infinity-to-infinity, （5） infinity-to-surface, （6） infinity-to-surrounding, （7） surface-to-surface, （8） surface- to-surrounding, and （9） surface-to- infinity. Assume that the orbital elements are constant near the periapsis, the gravitation potential is expanded into a harmonic series. Then, the influence of the gravitational perturbation on the orbit is studied analytically. The styles of orbits are dependent on the argument of periapsis, the periapsis radius, and the periapsis velocity. Given the argument of periapsis, the orbital energy before and after perturbation can be derived according to the periapsis radius and the periapsis velocity. Simulations have been performed for orbits in the gravitational field of 216 Kleopatra. The numerical results are well consistent with analytic predictions.     

Complex dynamics in the stretch-twist-fold flow

The system associated with fluid particle motions of the stretch-twist-fold (STF) flow has displayed rich and attractive dynamic properties. Detailed research on the system has been done in this work. By using a high-dimensional generalization of the Melnikov method, the explicit parametric conditions for the existence of periodic solutions in the system can be determined. Then, by using the new-KAM-like theorems for perturbations of a three-dimensional generalized Hamiltonian system, the criteria for the existence of invariant tori in the STF flow have been obtained. In addition, one new first integral is found. On the basis of it, nonexistence of chaos in the system at α=0 is rigorously proved. Nonexistence of homoclinic orbits is also proved in the system if some conditions hold. And an interesting phenomenon is found, where the unit circle in the (y,z)-plane is filled with heteroclinic orbits of the system at α=0. The system with α=0 is also successfully reduced to a generalized Hamiltonian system, and further transformed to slowly varying oscillators.     

Chaos in a pendulum with feedback control

We study chaotic dynamics of a pendulum subjected to linear feedback control with periodic desired motions. The pendulum is assumed to be driven by a servo-motor with small inductance, so that the feedback control system reduces to a periodic perturbation of a planar Hamiltonian system. This Hamiltonian system can possess multiple saddle points with non-transverse homoclinic and/or heteroclinic orbits. Using Melnikov's method, we obtain criteria for the existence of chaos in the pendulum motion. The computation of the Melnikov functions is performed by a numerical method. Several numerical examples are given and the theoretical predictions are compared with numerical simulation results for the behavior of invariant manifolds.     

万有引力场中陀螺体的混沌运动

研究万有引力场中沿圆轨道运行的非对称陀螺体的姿态运动,引入Ｄｅｐｒｉｔ正则变量建立系统的Ｈａｍｉｌｔｏｎ结构,利用Ｍｅｌｎｉｋｏｖ方法证明在万有引力短作用的昆体产生混沌运动的可能性。对Ｐｏｉｎｃａｒｅ截面的数值计算表明提高陀螺体的转子转速可对混沌起抑制作用。     

Analysis method of chaos and sub-harmonic resonance of nonlinear system without small parameters

The Melnikov method is important for detecting the presence of transverse homoclinic orbits and the occurrence of homoclinic bifurcations.Unfortunately,the traditional Melnikov methods strongly depend on small parameters,which do not exist in most practical systems.Those methods are limited in dealing with the systems with strong nonlinearities.This paper presents a procedure to study the chaos and sub-harmonic resonance of strongly nonlinear practical systems by employing a homotopy method that is used ...     

Homoclinic orbits in a shallow arch subjected to periodic excitation

Homoclinic orbits in a shallow arch subjected to periodic excitation are investigated in the presence of 1:1 internal resonance and external resonance. The method of multiple scales is used to obtain a set of near-integrable systems. The geometric singular perturbation method and Melnikov method are employed to show the existence of the one-bump and multi-bump homoclinic orbits that connect equilibria in a resonance band of the slow manifold. These orbits arise from singular homoclinic orbits and are composed of alternating slow and fast pieces. The result obtained imply the existence of the amplitude-modulated chaos for the Smale horseshoe sense in the class of shallow arch systems.     

时变小扰动Hamilton系统的Hopf分岔

运用Melnikov方法研究了时变小扰动Ｈamilton系统周期轨道发生Ｈopf分岔的条件,并将这些条件应用到一类三维时变小扰动非自治系统,数值结果验证了本文理论的正确性,进一步数值积分表明,所研究的系统还存在复杂而有规律的环面分岔行为。     

Chaotic pitch motion of an asymmetric non-rigid spacecraft with viscous drag in circular orbit

We study the pitch motion dynamics of an asymmetric spacecraft in circular orbit under the influence of a gravity gradient torque. The spacecraft is perturbed by a small aerodynamic drag torque proportional to the angular velocity of the body about its mass center. We also suppose that one of the moments of inertia of the spacecraft is a periodic function of time. Under both perturbations, we show that the system exhibits a transient chaotic behavior by means of the Melnikov method. This method gives us an analytical criterion for heteroclinic chaos in terms of the system parameters which is numerically contrasted. We also show that some periodic orbits survive for perturbation small enough.     

Singular perturbation theory for homoclinic orbits in a class of near-integrable Hamiltonian systems

This paper describes a new type of orbits homoclinic to resonance bands in a class of near-integrable Hamiltonian systems. It presents a constructive method for establishing whether small conservative perturbations of a family of heteroclinic orbits that connect pairs of points on a circle of equilibria will yield transverse homoclinic connections between periodic orbits in the resonance band resulting from the perturbation. In any given example, this method may be used to prove the existence of such transverse homoclinic orbits, as well as to determine their precise shape, their asymptotic behavior, and their possible bifurcations. The method is a combination of the Melnikov method and geometric singular perturbation theory for ordinary differential equations.     

Determining the nature of orbits in disk galaxies with non-spherical nuclei

We investigate the regular or chaotic nature of orbits of stars moving in the meridional plane (R,z) of an axially symmetric galactic model with a flat disk and a central, non-spherical and massive nucleus. In particular, we study the influence of the flattening parameter of the central nucleus on the nature of orbits, by computing in each case the percentage of chaotic orbits, as well as the percentages of orbits of the main regular families. In an attempt to maximize the accuracy of our results upon distinguishing between regular and chaotic motion, we use both the Fast Lyapunov Indicator (FLI) and the Smaller ALingment Index (SALI) methods to extensive samples of orbits obtained by integrating numerically the equations of motion as well as the variational equations. Moreover, a technique which is based mainly on the field of spectral dynamics that utilizes the Fourier transform of the time series of each coordinate is used for identifying the various families of regular orbits and also to recognize the secondary resonances that bifurcate from them. Varying the value of the flattening parameter, we study three different cases: (i) the case where we have a prolate nucleus, (ii) the case where the central nucleus is spherical, and (iii) the case where an oblate massive nucleus is present. Furthermore, we present some additional findings regarding the reliability of short time (fast) chaos indicators, as well as a new method to define the threshold between chaos and regularity for both FLI and SALI, by using them simultaneously. Comparison with early related work is also made.     

Hamiltonian structure for rigid body with flexible attachments in a circular orbit

Hamiltonian structure of a rigid body in a circular orbit is established in this paper. With the reduction technique, the Hamiltonian structure of a rigid body in a circular orbit is derived from Lie-Poisson structure of semidirect product, and Hamiltonian is derived from Jacobi's integral. The above method can be generalized to establish the Hamiltonian structure of a rigid body with a flexible attachment in a circular orbit. At last, an example of stability analysis is given. The project supported by National Natural Science Foundation of China and Aeronautic Science Foundation.     

Limit Cycles Near Homoclinic and Heteroclinic Loops

In the study of near-Hamiltonian systems, the first order Melnikov function plays an important role. It can be used to study Hopf, homoclinic and heteroclinic bifurcations, and the so-called weak Hilbert’s 16th problem as well. The form of expansion of the first order Melnikov function at the Hamiltonian value h ₀ such that the curve defined by the equation H(x, y) = h ₀ contains a homoclinic loop has been known together with the first three coefficients of the expansion. In this paper, our main purpose is to give an explicit formula to compute the first four coefficients appeared in the expansion of the first order Melnikov function at the Hamiltonian value h ₀ such that the curve defined by the equation H(x, y) = h ₀ contains a homoclinic or heteroclinic loop, where the formula for the fourth coefficient is new, and to give a way to find limit cycles near the loops by using these coefficients. As an application, we consider polynomial perturbations of degree 4 of quadratic Hamiltonian systems with a heteroclinic loop, and find 3 limit cycles near the loop. Dedicated to Professor Zhifen Zhang on the occasion of her 80th birthday     

The Fast Norm Vector Indicator (FNVI) method: a new dynamical parameter for detecting order and chaos in Hamiltonian systems

In the present article, we introduce and also deploy a new, simple, very fast, and efficient method, the Fast Norm Vector Indicator (FNVI) in order to distinguish rapidly and with certainty between ordered and chaotic motion in Hamiltonian systems. This distinction is based on the different behavior of the FNVI for the two cases: the indicator after a very short transient period of fluctuation displays a nearly constant value for regular orbits, while it continues to fluctuate significantly for chaotic orbits. In order to quantify the results obtained by the FNVI method, we establish the dFNVI, which is the quantified numerical version of the FNVI. A thorough study of the method??s ability to achieve an early and clear detection of an orbit??s behavior is presented both in two and three degrees of freedom (2D and 3D) Hamiltonians. Exploiting the advantages of the dFNVI method, we demonstrate how one can rapidly identify even tiny regions of order or chaos in the phase space of Hamiltonian systems. The new method can also be applied in order to follow the time evolution of sticky orbits. A detailed comparison between the new FNVI method and some other well-known dynamical methods of chaos detection reveals the great efficiency and the leading role of this new dynamical indicator.     

一类非自治动力系统超次谐周期解的不存在性

在非线性动力系统的研究中， Melnikov函数被广泛地用来作为微扰哈密顿系统是否发生次谐或超次谐分岔乃至混沌的判 据. 但是在大多数情况下，经典的Melnikov方法往往只给出存在次谐周期解的结论. 产生 该结果的原因被归之为在经典的Melnikov方法中只采取了一阶近似，因而高阶Melnikov方 法被发展用来判断超次谐周期解的存在性. 本文对一类非自治微分动力系统进行了研究，证 明了在这样一类系统中如果存在周期解则只可能是次谐周期解，超次谐周期解不可能存在， 并进一步证明了在一类平面问题中所定义的旋转(R)型超次谐周期解同样不可能存在.作为 该结论的一个应用，文中考察了几个典型的算例，结果表明现有的二阶Melnikov方法判断 平面扰动系统是否存在超次谐周期解的结论是不恰当的，并提供了一个简单的几何上的解释.     

Global Bifurcations in Parametrically Excited Systems with Zero-to-One Internal Resonance

In this work we study the existence of Silnikov homoclinicorbits in the averaged equations representing the modal interactionsbetween two modes with zero-to-one internal resonance. The fast mode isparametrically excited near its resonance frequency by a periodicforcing. The slow mode is coupled to the fast mode when the amplitude ofthe fast mode reaches a critical value so that the equilibrium of theslow mode loses stability. Using the analytical solutions of anunperturbed integrable Hamiltonian system, we evaluate a generalizedMelnikov function which measures the separation of the stable and theunstable manifolds of an annulus containing the resonance band of thefast mode. This Melnikov function is used together with the informationof the resonances of the fast mode to predict the region of physicalparameters for the existence of Silnikov homoclinic orbits.     

Homoclinic Tubes in Discrete Nonlinear Schrödinger Equation under Hamiltonian Perturbations

In this paper, we study the discrete cubic nonlinear Schrödinger lattice under Hamiltonian perturbations. First we develop a complete isospectral theory relevant to the hyperbolic structures of the lattice without perturbations. In particular, Bäcklund–Darboux transformations are utilized to generate heteroclinic orbits and Melnikov vectors. Then we give coordinate-expressions for persistent invariant manifolds and Fenichel fibers for the perturbed lattice. Finally based upon the above machinery, existence of codimension 2 transversal homoclinic tubes is established through a Melnikov type calculation and an implicit function argument. We also discuss symbolic dynamics of invariant tubes each of which consists of a doubly infinite sequence of curve segments when the lattice is four dimensional. Structures inside the asymptotic manifolds of the transversal homoclinic tubes are studied, special orbits, in particular homoclinic orbits and heteroclinic orbits when the lattice is four dimensional, are studied.     

A simple feedback control system: Bifurcations of periodic orbits and chaos

We consider a pendulum subjected to linear feedback control with periodic desired motions. The pendulum is assumed to be driven by a servo-motor with small time constant, so that the feedback control system can be approximated by a periodically forced oscillator. It was previously shown by Melnikov's method that transverse homoclinic and heteroclinic orbits exist and chaos may occur in certain parameter regions. Here we study local bifurcations of harmonics and subharmonics using the second-order averaging method and Melnikov's method. The Melnikov analysis was performed by numerically computing the Melnikov functions. Numerical simulations and experimental measurements are also given and are compared with the previous and present theoretical predictions. Sustained chaotic motions which result from homoclinic and heteroclinic tangles for not only single but also multiple hyperbolic periodic orbits are observed. Fairly good agreement is found between numerical simulation and experimental results.