Separatrices and limit cycles of strongly nonlinear oscillators by the perturbation-incremental method

The perturbation-incremental method is applied to determine the separatrices and limit cycles of strongly nonlinear oscillators. Conditions are derived under which a limit cycle is created or destroyed. The latter case may give rise to a homoclinic orbit or a pair of heteroclinic orbits. The limit cycles and the separatrices can be calculated to any desired degree of accuracy. Stability and bifurcations of limit cycles will also be discussed.     

A non-linear seales method for strongly non-linear oscillators

A non-linear seales method is presented for the analysis of strongly non-linear oseillators of the form % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaamXvP5wqonvsaeHbfv3ySLgzaGqbdiqb-Hha4zaadaGaey4kaSIa% am4zaiaacIcacqWF4baEcaGGPaGae8xpa0JaeqyTduMaamOzaiaacI% cacqWF4baEcqWFSaalcuWF4baEgaGaaiaabMcaaaa!4FEC!\[\ddot x + g(x) = \varepsilon f(x,\dot x{\text{)}}\], where g(x) is an arbitrary non-linear function of the displacement x. We assumed that % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaamXvP5wqonvsaeHbfv3ySLgzaGqbdiab-Hha4jaacIcacqWF0baD% cqWFSaalcqaH1oqzcaGGPaGaeyypa0Jae8hEaG3aaSbaaSqaaiaaic% daaeqaaOGaaiikaiabe67a4jaacYcacqaH3oaAcaGGPaGaey4kaSYa% aabmaeaacqaH1oqzdaahaaWcbeqaaiaad6gaaaaabaGaamOBaiabg2% da9iaaigdaaeaacaWGTbGaeyOeI0IaaGymaaqdcqGHris5aOGae8hE% aG3aaSbaaSqaaiab-5gaUbqabaGccaGGOaGaeqOVdGNaaiykaiabgU% caRiaad+eacaGGOaGaeqyTdu2aaWbaaSqabeaacaWGTbaaaOGaaiyk% aaaa!67B9!\[x(t,\varepsilon ) = x_0 (\xi ,\eta ) + \sum\nolimits_{n = 1}^{m - 1} {\varepsilon ^n } x_n (\xi ) + O(\varepsilon ^m )\], where % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaabsgacqaH+oaEcaGGVaGaaeizaiaadshacqGH9aqpdaaeWaqa% aiabew7aLnaaCaaaleqabaGaamOBaaaaaeaacaWGUbGaeyypa0JaaG% ymaaqaaiaad2gaa0GaeyyeIuoakiaadkfadaWgaaWcbaGaamOBaaqa% baGccaGGOaGaeqOVdGNaaiykaaaa!4FFC!\[{\text{d}}\xi /{\text{d}}t = \sum\nolimits_{n = 1}^m {\varepsilon ^n } R_n (\xi )\], % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiaabsgacqaH3oaAcaGGVaGaaeizaiaadshacqGH9aqpdaaeWaqa% aiabew7aLnaaCaaaleqabaGaamOBaaaaaeaacaWGUbGaeyypa0JaaG% imaaqaaiaad2gaa0GaeyyeIuoakiaadofadaWgaaWcbaGaamOBaaqa% baGccaGGOaGaeqOVdGNaaiilaiabeE7aOjaacMcaaaa!5241!\[{\text{d}}\eta /{\text{d}}t = \sum\nolimits_{n = 0}^m {\varepsilon ^n } S_n (\xi ,\eta )\], and R _n,S _nare to be determined in the course of the analysis. This method is suitable for the systems with even non-linearities as well as with odd non-linearities. It can be viewed as a generalization of the two-variable expansion procedure. Using the present method we obtained a modified Krylov-Bogoliubov method. Four numerical examples are presented which served to demonstrate the effectiveness of the present method.     

Periodic solutions of strongly non-linear Mathieu oscillators

The purpose of this paper is to continue our investigation into periodic solutions of strongly non-linear Mathieu oscillators. The modified version of the generalized averaging method which we developed recently is applied to derive highly accurate analytical expressions for these periodic solutions. These analytical results are used, together with the perturbation methods of multiple time scaling, to obtain second order expressions for the stability regions of these periodic solutions. The analytical research results are verified with numerical computations. Very good agreement is found, which shows the applicability of the modified version of the generalized averaging method to this kind of strongly non-linear oscillators. These oscillators may be used to model the beam-beam interaction in particle accelerators.     

A Perturbation-Incremental Method for Strongly Nonlinear Autonomous Oscillators with Many Degrees of Freedom

The perturbation-incremental method is extended to determine thebifurcations and limit cycles of strongly nonlinear autonomousoscillators with many degrees of freedom. Coupled van der Poloscillators and coupled Rayleigh oscillators are taken as numericalexamples. Limit cycles of the oscillators can be calculated to anydesired degree of accuracy. The stabilities of limit cycles are alsodiscussed.     

Response probability density functions of strongly non-linear systems by the path integration method

The paper discusses challenges in numerical analysis and numerical/analytical results for strongly non-linear systems—systems with “signum”-type non-linearities. Such non-linearities are implemented for instantaneous variations of the systems’ parameters, to reduce their mean energy response when subjected to random excitations. Numerical results for displacement and velocity response probability density functions (PDFs), energy response PDFs and various order moments are obtained by the path integration technique. Attention is also given to evaluation of mean upcrossing rate, related to the system's half period, via Rice's formula informally applied to discontinuous response PDFs.     

含高阶余项的非线性动力系统数值计算法

建立了一种求解非线性动力系统高精度数值计算的新方法,重构了等价的非线性动力系统方程,该方程考虑了非线性函数的任意高阶项,并给出了该方程的Duhamel积分表达式,在时间步长内用Newton-Raphson法进行数值迭代求解,该方法能连续满足微分方程而不只是在离散的步长端点满足方程,从而打破了传统的Euler型有限差分法。计算实例表明,该方法计算精度高于传统的Runge-Kutta,Newmark-β和Wilson-θ等方法。     

An equivalent nonlinearization method for strongly nonlinear oscillations

An equivalent nonlinearization method is proposed for the study of certain kinds of strongly nonlinear oscillators. This method is to express the nonlinear restored force of an oscillatory system by a polynomial of degree two or three such that the asymptotic solutions can be derived in terms of elliptic functions. The least squares method is used to determine the coefficients of approximate polynomials. The advantage of present method is that it is valid for relatively large oscillations. As an application, a strongly nonlinear oscillator with slowly varying parameters resulted from free-electron laser is studied in detail. Comparisons are made with other methods to assess the accuracy of the present method.     

A new method for detecting non-linear damping and restoring forces in non-linear oscillation systems from transient data

The purpose of this study is to recover the functional form of both non-linear damping and non-linear restoring forces in the non-linear oscillatory motions of an autonomous system. Using two sets of measured motion response data of the system, an inverse problem is formulated for recovering (or identification): the differential equation of motion is transformed into an equivalent integral equation of motion. The identification, which is non-linear, is shown to be one-to-one. However, the inverse problem formulated herein is concerned with the Volterra-type of non-linear integral equation of the first kind. This leads to numerical instability: solutions of the inverse problem lack stability properties. In order to overcome the difficulty, a regularization method is applied to the identification process. In addition, an L-curve criterion, combined with regularization, is introduced to find an optimal choice for the regularization parameter (i.e., the number of iterations), in the presence of noisy data. The workability of the identification is investigated for simultaneously recovering the functional form of the non-linear damping and the non-linear restoring forces through a numerical experiment.     

A method for finding the invariants and exact solutions of coupled non-linear differential equations with applications to dynamical systems

The polynomial invariants of (a set) non-linear differential equations are found by using a direct approach. The integrability of these invariants deserves the integrability of the given set of coupled differential equations. As applications, the Lorenz and Rikitake sets, among others, are studied. New invariants are obtained.     

A new method for the non-linear deflection analysis of an infinite beam resting on a non-linear elastic foundation

The aim of this paper is to develop a new method of analyzing the non-linear deflection behavior of an infinite beam on a non-linear elastic foundation. Non-linear beam problems have traditionally been dealt with by semi-analytical approaches that involve small perturbations or by numerical methods, such as the non-linear finite element method. In this paper, in contrast, a transformed non-linear integral equation that governs non-linear beam deflection behavior is formulated to develop a new method for non-linear solutions. The proposed method requires an iteration to solve non-linear problems, but is fairly simple and straightforward to apply. It also converges quickly, whereas traditional non-linear solution procedures are generally quite complex in application. Mathematical analysis of the proposed method is performed. In addition, illustrative examples are presented to demonstrate the validity of the method developed in the present study.     

A novel method of solving the exact equations for capillary-gravity waves

A new method is presented for the computation of two-dimensional periodicprogressive surface waves propagating under the combined influence of gravity and surfacetension.The nonlinear surface is expressed by Fourier series with finite number of terms,after the computational domain is transformed into a unit circle.The dynamic boundaryequation is used in its exact nonlinear form and the coefficients of Fourier series are foundby the Nweton-Raphson method successively.This is a neat method,Yielding highprescision with little computational effort.     

Three-dimensional arbitrary Lagrangian–Eulerian numerical prediction method for non-linear free surface oscillation

A numerical prediction method has been proposed to predict non-linear free surface oscillation in an arbitrarily-shaped three-dimensional container. The liquid motions are described with Navier–Stokes equations rather than Laplace equations which are derived by assuming the velocity potential. The profile of a liquid surface is precisely represented with the three-dimensional curvilinear co-ordinates which are regenerated in each computational step on the basis of the arbitrary Lagrangian–Eulerian (ALE) formulation. In the transformed space, the governing equations are discretized on a Lagrangian scheme with sufficient numerical accuracy and the boundary conditions near the liquid surface are implemented in a complete manner. In order to confirm the applicability of the present computational technique, numerical simulations are conducted for the free oscillations of viscid and inviscid liquids and for highly non-linear oscillation. In addition, non-linear sloshing motions caused by horizontal and vertical excitations and a transition from non-linear sloshing to swirling are numerically predicted in three-dimensional cylindrical containers. Conclusively, it is shown that these sloshing motions associated with high non-linearity are reasonably predicted with the present numerical technique. © 1998 John Wiley & Sons, Ltd.     

A natural neighbour method based on Fraeijs de Veubeke variational principle for materially non-linear problems

Abstract The natural neighbour method can be considered as one of many variants of the meshless methods. In the present paper, a new approach based on the Fraeijs de Veubeke （FdV） functional, which is initially developed for linear elasticity, is extended to the case of geometrically linear but materially non-linear solids. The new approach provides an original treatment to two classical problems： the numerical evaluation of the integrals over the domain A and the enforcement of boundary conditions of the type ui = hi on Su. In the absence of body forces （Fi = 0）, it will be shown that the calculation of integrals of the type fA .dA can be avoided and that boundary conditions of the type ui = hi on Su can be imposed in the average sense in general and exactly if hi is linear between two contour nodes, which is obviously the case for tTi = O.     

A semi-analytical approach for the non-linear large deflection analysis of laminated rectangular plates under general out-of-plane loading

A semi-analytical approach for the geometrically non-linear analysis of rectangular laminated plates with general inplane and out-of-plane boundary conditions under a general distribution of out-of-plane loads is developed. The analysis is based on the elastic thin plate theory with geometrically non-linear von Kármán strains. The solution of the non-linear partial differential equations is reduced to an iterative sequential solution of non-linear ordinary differential equations using the multi-term extended Kantorovich method. The efficiency, accuracy, and convergence of the proposed method are examined through a comparison with other semi-analytical methods and with finite element analyses. The capabilities of the approach and its applicability to the non-linear large deflection analysis of plate structures are demonstrated through various numerical examples. Emphasis is placed on combinations of lamination, boundary, and loading conditions that cannot be analyzed using alternative semi-analytical methods.