Van der Pol时滞弱耦合系统多尺度分析

对于非线性耦合项中带有时滞的van der Pol系统,采用多尺度法对该系统进行定性以及定量的分析.研究结果表明,对于van der Pol时滞耦合系统,时滞的存在影响了系统的稳定性,使系统的周期解发生了静态分岔和Hopf分岔.研究还发现,对于耦合强度较弱的情形,利用多尺度法对系统进行定嚣分析是合理可靠的.我们取不同的耦合强度作用了系统的时间历程图,相图和分岔图,分析了解析解与数值解之间产生误差的原因.本文所研究的系统来源于耦合的激光振荡器.     

含外激励van der Pol-Mathieu方程的非线性动力学特性分析

本文针对含有自激励,参数激励和外激励等三种激励联合作用下van der Pol-Mathieu方程的周期响应和准周期运动进行分析,发现其准周期运动的频谱中含有均匀边频带这一新的特性.首先,采用传统的增量谐波平衡法(IHB法)分析了van der Pol-Mathieu方程的周期响应,得到了其非线性频率响应曲线;再利用F...     

非线性动力系统线性模型数值计算的Taylor变换法

将非线性动力系统化为连续变化的线性系统,并导出任意自治或非自治非线性动力系统的瞬时线性化方程,该线性方程的连续变化描述了系统的全部复杂动力行为.进一步采用Taylor变换法求解系统的线性化方程,得到一种非线性动力系统数值计算的新方法,避免了指数矩阵展开的乘积运算.计算实例表明该方法在不增加计算机时的前提下,精度高于传统的Houbolt,Wilson-θ及Newmark-β等方法.计算了Duffing方程和van　Pol方程的混沌及周期特性.     

小参数法在悬垂索共振中的应用

当面外横向振动和面内横向振动频率的比接近1:2时,悬索会出现面内和面外耦合共振现象。为了研究悬索这种复杂独特的非线性特性,利用多尺度法对谐波激励的悬索动力学方程进行求解,得到对应于不同阶小量的偏微分方程组,其中二阶小量偏微分方程中的久期项不为0;采用提出的小参数法可以得到由久期项引起的悬索振动形态,解决久期项频率与系统频率相同但不能直接求解的问题;为了证明小参数法的准确性,采用Galerkin方法离散悬垂索的运动方程,然后利用多尺度法求解离散的运动方程,得到采用基函数描述的由久期项引起的连续系统的振动形态,与小参数法结论一致。     

多自由度非线性系统的同伦分析方法

大部分工程实际问题可以用多自由度非线性系统来描述,这些系统的数学模型是许多个耦合的两阶常微分方程.一般地,要精确求解这些方程非常困难,因此可以考虑它们的解析近似解.同伦分析方法是解非线性系统响应的有用工具,本文将它应用于多自由度非线性系统的求解中.利用求两自由度耦合van del Pol振子周期解的实例,展示了同伦分析方法的有效性和巨大潜力.同时,把得到的解析近似解与系统的Runge-Kutta数值解作了比较,结果表明同伦分析方法是求解多自由度非线性系统的有效方法.     

两自由度耦合van der Pol振子的拟主振动解

本文运用非线性系统的模态方法研究了两自由度耦合van der Pol振子。从退化系统稳定的主振动解出发,得到了原系统的拟主振动解,并给出了系统周期运动的条件,讨论了系统周期解、概周期解的分叉。     

悬索的多重内共振研究

本文对谐波激励的悬索的非线性响应进行了研究,同时考虑了如下问题(1):面内第三阶对称模态的主共振:(2):面内第一阶、第三阶对称模态和面外第五阶模态之间的内共振.本方首先针对考虑大变形的悬索动力学方程,由线性理论求得各阶频率,考察可能出现的内共振.然后利用直接法对悬索的运动学方程和边界条件进行非线性求解.由多尺度法得到系统的平均方程和悬索响应的二阶近似解.随后利用Newton-Raphson 方法和弧长法对特定张拉索进行数值仿真计算,得到面内第一阶对称模态、面内第三阶对称模态和面外第五阶模态的稳态解,并分析了解的稳定性.绘制幅频响应曲线,发现了关于悬索响应的多种分叉现象,并且对各种分叉现象周期解、混沌解进行了讨论.     

轴向运动正交各向异性板的超谐波共振及稳定性

研究了轴向运动正交各向异性条形薄板在线载荷作用下的超谐波共振问题。通过哈密顿原理导出了几何非线性下正交各向异性条形板的非线性振动方程。运用伽辽金积分法,推得了关于时间变量的量纲归一化非线性振动微分方程组。应用多尺度法求解三阶超谐波共振问题,得到了稳态运动下一阶、二阶、三阶共振形式的共振幅值响应方程。利用Liapunov方法推得不同共振形式稳态解的稳定性判据,并据此分析不同参数对系统稳定性的影响。绘制了振幅特性变化曲线图和与之对应的激发共振多解临界点曲线图,分析系统参数对共振的影响,并预测系统进入非线性共振区域的临界条件。得出激励在特定位置区间时可激发系统的超谐波共振,随着激励幅值的增加,上稳定解支减小,下稳定解支增加,且一阶模态振幅大于二阶、三阶振幅。     

强非线性振动系统求解的两种解析方法

本文给出了拟保守系统的渐近解,该解是在非线性解的基础上的摄动解,因而可求解强非线性振动系统。文中利用此解研究了具有强非对称恢复力项的Liénard方程的极限环问题,给出了各种特殊情况下该方程的极限环幅值的计算公式,并讨论了非线性恢复力项对极限环的影响。此外,本文提出了一种改进的谐波平衡法,该方法是谐波平衡法与伽辽金方法结合的产物。     

三类随机系统广义概率密度演化方程的解析解

近年来逐步发展的概率密度演化方法理论为随机动力系统的分析与控制研究提供了新的途径.过去若干年来,已经发展了一系列数值方法如有限差分法、无网格法用于求解广义概率密度演化方程.但是,针对典型随机系统,关于这一方程解析解尚比较缺乏.本文以李群方法为工具,研究给出了Van der Pol振子、Riccati方程和Helmholtz振子3类典型随机非线性系统的广义概率密度演化方程解析解.这些结果,不仅可以作为检验求解广义概率密度演化方程的数值方法结果正确性的判别依据,也为概率密度演化理论的进一步深入研究提供了若干分析实例.      

Response statistics of van der Pol oscillators excited by white noise

The joint probability density function of the state space vector of a white noise exoited van der Pol oscillator satisfics a Fokker-Planck-Kolmogorov (FPK) equation. The paper describes a numerical procedure for solving the transient FPK equation based on the path integral solution (PIS) technique. It is shown that by combining the PIS with a cubic B-spline interpolation method, numerical solution algorithms can be implemented giving solutions of the FPK equation that can be made accurate down to very low probability levels. The method is illustrated by application to two specific examples of a van der Pol oscillator.     

A Modified Perturbation Technique Depending Upon an Artificial Parameter

In this paper, a modified perturbation method is proposed to search for analytical solutions of nonlinear oscillators without possible small parameters. An artificial perturbation equation is carefully constructed by embedding an artificial parameter, which is used as expanding parameter. It reveals that various traditional perturbation techniques can be powerfully applied in this theory. Some examples, such as the Duffing equation and the van der Pol equation, are given here to illustrate its effectiveness and convenience. The results show that the obtained approximate solutions are uniformly valid on the whole solution domain, and they are suitable not only for weak nonlinear systems, but also for strongly nonlinear systems. In applying the new method, some special techniques have been emphasized for different problems.     

强非线性动力系统的频率增量法

提出一类强非线性动力系统的暧时频率增量法,将描述动力系统的二阶常微分方程,化为以相位为自变量、瞬廛频率为未知函数的积分方程;用谐波平衡原理,将求解瞬时频率的积分问题,归结为求解以频率增量的Fourier系数为独立变量的线性代数方程组;给出了若干例子。     

Orthogonal trajectories for nonlinear equations

Abdelkader developed a method for solving the van der Pol equation numerically. He observed that the orthogonal trajectories of the van der Pol equation satisfy a Riccati equation, the solution to which can be expressed in terms of Bessel functions. In the present paper, all nonlinear equations whose orthogonal trajectories satisfy Riccati or generalized Riccati equations are determined. Further on, another representation of the trajectories of the van der Pol equation is derived. This representation appears to be simpler than Abdelkader's.     

Non-linear oscillator limit cycle analysis using a time transformation approach

A method is presented for the analysis of limit cycle behavior of autonomous non-linear oscillators characterized by second order ordinary differential equations containing a small parameter. The method differs from the classical perturbation methods in that the dependent variable is not expanded in a power series in the small parameter. Rather, a new independent variable is sought such that in its domain the motion is simple harmonic. Use of this time transformation technique to generate limit cycle phase portrait, amplitude and period is presented. We show results of the application of the method to the van der Pol oscillator, to an oscillator with quadratic damping, and to a modified van der Pol oscillator which is statically unstable in the limit of small motion.     

Primary resonance of fractional-order van der Pol oscillator

In this paper the primary resonance of van der Pol (VDP) oscillator with fractional-order derivative is studied analytically and numerically. At first the approximately analytical solution is obtained by the averaging method, and it is found that the fractional-order derivative could affect the dynamical properties of VDP oscillator, which is characterized by the equivalent damping coefficient and the equivalent stiffness coefficient. Moreover, the amplitude–frequency equation for steady-state solution is established, and the corresponding stability condition is also presented based on Lyapunov theory. Then, the comparisons of several different amplitude–frequency curves by the approximately analytical solution and the numerical integration are fulfilled, and the results certify the correctness and satisfactory precision of the approximately analytical solution. At last, the effects of the two fractional parameters, i.e., the fractional coefficient and the fractional order, on the amplitude–frequency curves are investigated for some typical excitation amplitudes, which are different from the traditional integer-order VDP oscillator.     

Constructing invariant tori for two weakly coupled van der Pol oscillators

An algorithm is developed for the construction of an invariant torus of a weakly coupled autonomous oscillator. The system is put into angular standard form. The determining equations are found by averaging and are solved for the approximate amplitudes of the torus. A perturbation series is then constructed about the approximate amplitudes with unknown coefficients as periodic functions of the angular variables. A sequence of solvable partial differential equations is developed for determining the coefficients. The algorithm is applied to a system of nonlinearly coupled van der Pol equations and the first order coefficients are generated in a straightforward manner. The approximation shows both good numerical accuracy and reproducibility of the periodicities of the van der Pol system. A comparitive analysis of integrating the van der Pol system with integrating the phase equations from the angular standard form on the approximate torus shows numerical errors of the order of the perturbation parameter =0.05 for integrations of up to 10,000 steps. Applying FFT to the numerical periodicities generated by integrating the van der Pol system near the tours reveals the same predominant frequencies found in the perturbation coefficients. Finally an expected rotation number is found by integrating the phase equations on the approximate torus.Contribution of the National Institute of Standards and Technology, a Federal agency.     

Remarks on the Perturbation Methods in Solving the Second-Order Delay Differential Equations

The paper presents a study on the validity of perturbation methods, suchas the method of multiple scales, the Lindstedt–Poincaré method and soon, in seeking for the periodic motions of the delayed dynamic systemsthrough an example of a Duffing oscillator with delayed velocityfeedback. An important observation in the paper is that the method ofmultiple scales, which has been widely used in nonlinear dynamics, worksonly for the approximate solutions of the first two orders, and givesrise to a paradox for the third-order approximate solutions of delaydifferential equations. The same problem appears when theLindstedt–Poincaré method is implemented to find the third-orderapproximation of periodic solutions for delay differential equations,though it is effective in seeking for any order approximation ofperiodic solutions for nonlinear ordinary differential equations. Apossible explanation to the paradox is given by the results obtained byusing the method of harmonic balance. The paper also indicates thatthese perturbation methods, despite of some shortcomings, are stilleffective in analyzing the dynamics of a delayed dynamic system sincethe approximate solutions of the first two orders already enable one togain an insight into the primary dynamics of the system.     

Constructing transient response probability density of non-linear system through complex fractional moments

The probability density function for transient response of non-linear stochastic system is investigated through the stochastic averaging and Mellin transform. The stochastic averaging based on the generalized harmonic functions is adopted to reduce the system dimension and derive the one-dimensional Itô stochastic differential equation with respect to amplitude response. To solve the Fokker–Plank–Kolmogorov equation governing the amplitude response probability density, the Mellin transform is first implemented to obtain the differential relation of complex fractional moments. Combining the expansion form of transient probability density with respect to complex fractional moments and the differential relations at different transform parameters yields a set of closed-form first-order ordinary differential equations. The complex fractional moments which are determined by the solution of the above equations can be used to directly construct the probability density function of system response. Numerical results for a van der Pol oscillator subject to stochastically external and parametric excitations are given to illustrate the application, the convergence and the precision of the proposed procedure.     

Approximate analytical solution in slow-fast system based on modified multi-scale method

A simple, yet accurate modified multi-scale method (MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method (MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simplified system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking (SP), spiking-quiescence (SP-QS), and quiescence (QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method (HBM). The results obtained by the present method are considerably better than those obtained by traditional methods, quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.