A modified and compact form of Krylov-Bogoliubov-Mitropolskii unified method for solving an nth order non-linear differential equation

A modified and compact form of Krylov-Bogoliubov-Mitropolskii (KBM) (Introduction to Nonlinear Mechanics, Princeton University Press, Princeton, NJ, 1947; Asymptotic Methods in the Theory of Nonlinear Oscillations, Gordan and Breach, New York, 1961) unified method (J. Franklin Inst. 339 (2002) 239) is determined for obtaining the transient response of an nth order (n?2) differential equation with small non-linearities. The formula presented in (J. Franklin Inst. 339 (2002) 239) is a changed form of KBM method. For n=2,3,4, some previous formulas were found separately by several authors in terms of amplitude and phase variables; but the formula of Shamsul Alam, J. Franklin Inst. 339 (2002) 239) is derived in terms of some unusual variables instead of amplitudes and phases. The formula of Shamsul Alam, J. Franklin Inst. 339 (2002) 239) is a general form and used arbitrarily to obtain asymptotic solution for n=2,3,4,…. However, a solution obtained by formula Shamsul Alam, J. Franklin Inst. 339 (2002) 239) is transformed to a formal form replacing the unusual variables by amplitude and phase variables. In the present paper, the formula of Shamsul Alam, J. Franklin Inst. 339 (2002) 239) is itself transformed to a usual form (i.e. in terms of amplitude and phase variables). The later form of the formula is similar to most of the previous formulas found by several authors when n=2,3,4. This form of the formula is also generalized and it is easier than those obtained in all previous papers (extension) and identical to that initiated by original contributors (Introduction to Nonlinear Mechanics, Princeton University Press, Princeton, NJ, 1947; Asymptotic Methods in the Theory of Nonlinear Oscillations, Gordan and Breach, New York, 1961).     

On second-order non-linear systems with conservative limit-cycles

A class of non-linear systems whose limit cycles are continuously conservative is suggested and dealt with. Most of the discussion is carried out by considering the behaviour of certain members of the class. Some significant properties of such systems are demonstrated. These properties appear important in the development of electronic function generators of high precision. The analysis of certain pumps is also related to models of similar systems.     

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

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

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.     

A unified compatibility method for exact solutions of non-linear flow models of Newtonian and non-Newtonian fluids

Criteria are established for higher order ordinary differential equations to be compatible with lower order ordinary differential equations. Necessary and sufficient compatibility conditions are derived which can be used to construct exact solutions of higher order ordinary differential equations subject to lower order equations. We provide the connection to generalized groups through conditional symmetries. Using this approach of compatibility and generalized groups, new exact solutions of non-linear flow problems arising in the study of Newtonian and non-Newtonian fluids are derived. The ansatz approach for obtaining exact solutions for non-linear flow models of Newtonian and non-Newtonian fluids is unified with the application of the compatibility and generalized group criteria.     

KBM unified method for solving an nth order non-linear differential equation under some special conditions including the case of internal resonance

The Krylov-Bogoliubov-Mitropolskii (KBM) unified method is used for obtaining the approximate solution of an nth order (n?4) ordinary differential equation with small non-linearities when a pair of eigen-values of the unperturbed equation is multiple (approximately or perfectly) of the other pair or pairs. The general solution can be used arbitrarily for over-damped, damped and undamped cases. In a damped or undamped case, one of the natural frequencies of the unperturbed equation may be a multiple of the other. Thus, the solution also covers the case of internal resonance which is an interesting and important part of non-linear oscillation. The determination of the solution is very simple and easier than the existing procedures developed by several authors (both in methods of averaging and multiple time scales) especially to tackle the case of internal resonance. The method is illustrated by an example of a fourth-order differential equation. The solution shows a good agreement with numerical solution in all of the three cases, e.g. over-damped, damped and undamped.     

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 new multistage spectral relaxation method for solving chaotic initial value systems

In this paper, we present a new pseudospectral method application for solving nonlinear initial value problems (IVPs) with chaotic properties. The proposed method, called the multistage spectral relaxation method (MSRM) is based on a novel technique of extending Gauss–Siedel type relaxation ideas to systems of nonlinear differential equations and using the Chebyshev pseudo-spectral methods to solve the resulting system on a sequence of multiple intervals. In this new application, the MSRM is used to solve famous chaotic systems such as the such as Lorenz, Chen, Liu, Rikitake, Rössler, Genesio–Tesi, and Arneodo–Coullet chaotic systems. The accuracy and validity of the proposed method is tested against Runge–Kutta and Adams–Bashforth–Moulton based methods. The numerical results indicate that the MSRM is an accurate, efficient, and reliable method for solving very complex IVPs with chaotic behavior.     

Multiple scale perturbation for second-order non-linear differential equations

Asymptotic solutions of a class of second-order non-linear differential equations with variable coefficients are studied. For large values of the parameter, the differential equations are of the singular-perturbation type and approximations are constructed by the generalized method of multiple scales.     

A probabilistic linearization method for non-linear systems subjected to additive and multiplicative excitations

A general method to obtain approximate solutions for the random response of non-linear systems subjected to both additive and multiplicative Gaussian white noises is presented. Starting from the concept of linearization, the proposed method of “Probabilistic Linearization” (PL) is based on the replacement of the Fokker–Planck equation of the original non-linear system with an equivalent one relative to a linear system subjected to additive excitation only. By means of the general scheme of the weighted residuals, the unknown coefficients of the equivalent system are determined. Assuming a Gaussian probability density function of the response process and by choosing the weighting functions in a suitable way, the equivalence of the proposed method, called “Gaussian Probabilistic Linearization” (GPL), with the “Gaussian Stochastic Linearization” (GSL) applied to the coefficients of the Itô differential rule is evidenced. In addition, the generalization of the proposed method, called “Generalized Gaussian Probabilistic Linearization” (GGPL), is presented. Numerical applications show as, varying the choice of the weighting functions, it is possible to obtain different linearizations, with a variable degree of accuracy. For the two examples considered, different suitable combinations of the weighting functions lead to different equivalent linear systems, all characterized by the exact solution in terms of variance.     

A frequency-domain method for solving linear time delay systems with constant coefficients

In an active control system, time delay will occur due to processes such as signal acquisition and transmission, calculation, and actuation. Time delay systems are usually described by delay differential equations (DDEs). Since it is hard to obtain an analytical solution to a DDE, numerical solution is of necessity. This paper presents a frequency-domain method that uses a truncated transfer function to solve a class of DDEs. The theoretical transfer function is the sum of infinite items expressed in terms of poles and residues. The basic idea is to select the dominant poles and residues to truncate the transfer function, thus ensuring the validity of the solution while improving the efficiency of calculation. Meanwhile, the guideline of selecting these poles and residues is provided. Numerical simulations of both stable and unstable delayed systems are given to verify the proposed method, and the results are presented and analysed in detail.     

A field method for solving the equations of motion of nonholonomic systems

In this paper, the field method^[1] for solving the equations of motion of holonomic nonconservative systems is extended to nonholonomic systems with constant mass and with variable mass. Two examples are given to illustrate its application. The project supported by the National Natural Science Foundation of China.     

PRECONDITIONED GAUSS-SEIDEL TYPE ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS

The preconditioned Gauss-Seidel type iterative method for solving linear systems, with the proper choice of the preconditioner, is presented. Convergence of the preconditioned method applied to Z-matrices is discussed. Also the optimal parameter is presented. Numerical results show that the proper choice of the preconditioner can lead to effective by the preconditioned Gauss-Seidel type iterative methods for solving linear systems.     

A symplectic finite element method based on Galerkin discretization for solving linear systems

We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems. For the dynamic responses of continuous medium structures, the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation. Thus, a symplectic finite element method with energy conservation is constructed in this paper. A linear elastic system can be discretized into multiple elements, and a Hamiltonian system of each element ...     

A direct method for non-linear normal modes

This paper presents a direct method for locating normal modes in certain holonomic, scleronomous, conservative non-linear two degree of freedom dynamical systems. The method does not require that the system studied be close to a linear system.     

A perturbation method for certain non-linear oscillators

We present a perturbation method for the analysis of single degree of freedom non-linear oscillation phenomena governed by an equation of motion containing a parameter ? which need not be small. The approach is to define a new parameter α = α(?) in such a way that asymptotic solutions in power series in α converge more quickly than do the standard perturbation expansions in power series in ?. Phenomena considered are free vibration of strongly non-linear conservative oscillators and steady state response of strongly non-linear oscillators subject to weak harmonic excitation.     

Construction of the stationary probability density for a family of SDOF strongly non-linear stochastic second-order dynamical systems

In this paper, a new procedure is proposed to construct the stationary probability density for a family of the single-degree-of-freedom (SDOF) strongly non-linear stochastic second-order dynamical systems subjected to parametric and/or external Gaussian white noises. First of all, the Fokker-Planck-Kolmogorov (FPK) equation associated with the original Itô stochastic differential equation is replaced by the equivalent FPK equation by adding arbitrary anti-symmetric diffusion coefficient. Then, a family of invariant measures depending on the arbitrary anti-symmetric diffusion coefficient and another arbitrary function is constructed by vanishing the probability flows in two directions. Finally, the drift vector associated with a family of Itô stochastic differential equations is deduced by giving, a priori, these two arbitrary functions. It is shown that the known invariant measures dependent on energy are only the special cases of invariant measures presented in this paper, while some other classes of invariant measures are independent of energy. Thus, the invariant measures constructed in this paper are those belonging to the most general class of the SDOF strongly non-linear stochastic second-order dynamical systems so far.     

求解非线性方程组的混沌分形方法

混沌分形是动力系统普遍出现的一种现象,牛顿-拉夫森NR(Newton-Raphson)方法是重要的一维及多维迭代技术,其迭代本身对初始点非常敏感,该敏感区是牛顿-拉夫森法所构成的非线性离散动力系统Julia集,在Julia集中迭代函数会呈现出混沌分形现象,提出了一种寻找牛顿-拉夫森函数的Julia点的求解方法,利用非线性离散动力系统在其Julia集出现混沌分形现象的特点,提出了一种基于牛顿-拉夫森法的非线性方程组求解的新方法,计算实例表明了该方法的有效性和正确性.     

A unified symmetry of lagrangian systems

In this paper, the definition and the criterion of a unified symmetry are presented. A new conserved quantity, as well as the Noether conserved quantity and the Hojman conserved quantity, deduced from the unified symmetry, are obtained. Some examples are given to illustrate the application of the results. The project supported by the National Natural Science Foundation of China (10272021)     

A METHOD FOR SOLVING THE DYNAMICS OF MULTIBODY SYSTEMS WITH RHEONOMIC AND NONHOLONOMIC CONSTRAINTS

ＡＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＴＨＥＤＹＮＡＭＩＣＳＯＦＭＵＬＴＩＢＯＤＹＳＹＳＴＥＭＳＷＩＴＨＲＨＥＯＮＯＭＩＣＡＮＤＮＯＮＨＯＬＯＮＯＭＩＣＣＯＮＳＴＲＡＩＮＴＳ￥ＳｈｕｉＸｉａｏｐｉｎｇ（水小平）ＺｈａｎｇＹｏｎｇｆａ（张永发）（Ｄｅｐａｒ...