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.     

The Asymptotic Perturbation Method for Nonlinear Continuous Systems

We apply the asymptotic perturbation (AP) method to the study of the vibrations of Euler--Bernoulli beam resting on a nonlinear elastic foundation. An external periodic excitation is in primary resonance or in subharmonic resonance in the order of one-half with an nth mode frequency. The AP method uses two different procedures for the solutions: introducing an asymptotic temporal rescaling and balancing the harmonic terms with a simple iteration. We obtain amplitude and phase modulation equations and determine external force-response and frequency-response curves. The validity of the method is highlighted by comparing the approximate solutions with the results of the numerical integration and multiple-scale methods.     

非线性结构动力学瞬态分析的摄动法

在用直接积分法求解非线性结构的动力响应时,常常需要做迭代运算。本文引入摄动方法后,加快了收敛速度,提高了计算效益。     

一类弱非线性振动问题的插值摄动解法

用插值摄动法［１］求解一类有阻尼的弱非线性振动问题。算例表明，当小参数很小时，本文结果和多尺度法的一级近似结果十分接近。当小参数不是很小时（即接近于强非线性振动时），本文结果，仍然相当准确，并优于多尺度法的一级近似结果。     

Types of Bifurcations of FitzHugh–Nagumo Maps

The FitzHugh–Nagumo-like systems are of fundamental importance to the understanding of the qualitative nature of nerve impulse propagation. Our work provides a numerical investigation of bifurcations associated with a family of piecewise differentiable canonical maps for a planar FitzHugh–Nagumo system. We describe the bifurcation structure of the maps with the variation of the parameters.     

Normal Modes and Boundary Layers for a Slender Tensioned Beam on a Nonlinear Foundation

In this paper, the nonlinear normal modes (NNMs) of a thin beamresting on a nonlinear spring bed subjected to an axial tension isstudied. An energy-based method is used to obtain NNMs. In conjunction with amatched asymptotic expansion, we analyze, through simple formulas, thelocal effects that a small bending stiffness has on the dynamics, alongwith the secular effects caused by a symmetric nonlinearity. Nonlinearmode shapes are computed and compared with those of the unperturbedlinear system. A double asymptotic expansion is employed to compute theboundary layers in the nonlinear mode shape due to the small bendingstiffness. Satisfactory agreement between the theoretical and numericalbackbone curves of the system in the frequency domain is observed.     

Perturbation Analysis for Wetting Fronts in Richard's Equation

Perturbation methods are used to study the interaction of wetting fronts with impervious boundaries in layered soils. Solutions of Richards' equation for horizontal and vertical infiltration problems are considered. Asymptotically accurate solutions are constructed from outer solutions and boundary-layer corrections. Results are compared with numerical simulations to demonstrate a high level of accuracy.     

Perturbation theory for nonlinear klein-gordon equation

ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｗｅｌｌ＿ｋｎｏｗｎｔｈａｔｔｈｅｎｏｎｌｉｎｅａｒＫｌｅｉｎ＿Ｇｏｒｄｏｎｅｑｕａｔｉｏｎｐｌａｙｓａｖｅｒｙｉｍｐｏｒｔａｎｔｒｏｌｅｉｎｎｏｎｌｉｎｅａｒｐｈｙｓｉｃｓｂｅｃａｕｓｅｏｆｉｔｓｗｉｄｅｓｐｒｅａｄａｐｐｌｉｃａｔｉｏｎｉｎｍａｎｙｆｉｅｌｄｓｏｆｐｈｙｓｉｃｓ[1,2 ].ＩｔｓｇｅｎｅｒａｌｆｏｒｍｉｓＵＴＴ-ＵＸＸ+ｍ2 Ｕ+λＵ3=0 ,( 1 )ｗｈｅｒｅｔｈｅｓｕｂｓｃｒｉｐｔｓｓｔａｎｄｆｏｒｐａｒｔｉａｌｄｉｆｆｅｒｅｎｔｉａｔｉｏｎｗｉｔｈｒｅｓｐ…     

求解非线性方程组Newton－PCG并行算法

基于向量机（ＶｅｃｔｏｒＣｏｍｐｕｔｅｒ），应用一次迭代法和不完全的Ｃｈｏｌｅｓｋｉ不塞入分解原理，建立了求解非线性方程组的Ｎｅｗｔｏｎ－ＰＣＧ并行算法。     

线性哈密顿系统的本征摄动法

利用哈密顿算子辛自共轭的特点讨论了保守哈密顿体系的摄动问题,给出了哈密顿矩阵的本征值与本征向量的二阶摄动分析方法。即当系统在哈密顿框架下进行较小修改时,不重复求解大型哈密顿矩阵的本征问题,只需在原系统的模态参数基础上进行模态分析即可,这种矩阵摄动法给出了修改后矩阵的二阶本征值和本征向量,为一般线性保守体系的本征摄动求解提出了一个新方法。     

结动振动摄动分析的新方法

提出了一种用于结构动力修改的设计灵敏度分析的新方法。它发展了Ｎｅｌｓｏｎ［１］陈［２］等人的方法，较好地解决了结构修改量大时计算精度低之间的矛盾。数值算例表明，本文新方法计算精度高，易于计算机实施。     

扰动法在结构分枝失稳分析中的应用

对结构进行平衡路径的跟踪分析，是全面了解该结构的受力性能所必须进行的一项工作。目前的结构非线性稳定分析技术一般仅对极值点失稳型问题较为有效，而对分枝点失稳型问题则困难较多。对于具有缺陷敏感性的结构，如拱结构、壳体等，在普通荷载作用下，其失稳路径常包含分枝点。文献[1]提出并认为位移扰动法和力扰动法在分析结构分枝失稳时具有很好的效果。本文采用多个不同类型的算例，对扰动法在结构分枝失稳问题中的应用进行了分析比较，表明该方法具有较强的跟踪能力。最后，就扰动法在结构分枝失稳问题中的应用提出几点建议。     

Experimental Nonlinear Identification of a Single Mode of a Transversely Excited Beam

A procedure is presented for using a primary resonance excitation in experimentally identifying the nonlinear parameters of a model approximating the response of a cantilevered beam by a single mode. The model accounts for cubic inertia and stiffness nonlinearities and quadratic damping. The method of multiple scales is used to determine the frequency-response function for the system. Experimental frequency- and amplitude-sweep data is compared with the prediction of the frequency-response function in a least-squares curve-fitting algorithm. The algorithm is improved by making use of experimentally known information about the location of the bifurcation points. The method is validated by using the extracted parameters to predict the force-response curves at other nearby frequencies.We then compare this technique with two other techniques that have been presented in the literature. In addition to the amplitude- and frequency-sweep technique presented, we apply a backbone curve- fitting technique and a time-domain technique to the second mode of a cantilevered beam. Differences in the parameter estimates are discussed. We conclude by discussing the limitations encountered for each technique. These include the inability to separate the nonlinear curvature and inertia effects and problems in estimating the coefficients of small terms with the time-domain technique.     

Nonlinear Normal Modes and Their Bifurcations for an Inertially Coupled Nonlinear Conservative System

This work concerns the nonlinear normal modes (NNMs) of a 2 degree-of-freedom autonomous conservative spring–mass–pendulum system, a system that exhibits inertial coupling between the two generalized coordinates and quadratic (even) nonlinearities. Several general methods introduced in the literature to calculate the NNMs of conservative systems are reviewed, and then applied to the spring–mass–pendulum system. These include the invariant manifold method, the multiple scales method, the asymptotic perturbation method and the method of harmonic balance. Then, an efficient numerical methodology is developed to calculate the exact NNMs, and this method is further used to analyze and follow the bifurcations of the NNMs as a function of linear frequency ratio p and total energy h. The bifurcations in NNMs, when near 1:2 and 1:1 resonances arise in the two linear modes, is investigated by perturbation techniques and the results are compared with those predicted by the exact numerical solutions. By using the method of multiple time scales (MTS), not only the bifurcation diagrams but also the low energy global dynamics of the system is obtained. The numerical method gives reliable results for the high-energy case. These bifurcation analyses provide a significant glimpse into the complex dynamics of the system. It is shown that when the total energy is sufficiently high, varying p, the ratio of the spring and the pendulum linear frequencies, results in the system undergoing an order–chaos–order sequence. This phenomenon is also presented and discussed.     

Nonlinear Response of a Parametrically Excited System Using Higher-Order Method of Multiple Scales

Two fundamentally different versions of the method of multiple scales (MMS) are currently in use in the study of nonlinear resonance phenomena. While the first version is the widely used reconstitution method, the second version is proposed by Rahman and Burton [1]. Both versions of the second-order MMS are applied to the differential equation obtained for a parametrically excited cantilever beam with a lumped mass at an arbitrary position. The bifurcation and stability of the obtained response show the difference between the two versions. While the Hopf bifurcation phenomena with no jump is found in the case of second-order MMS version I, both jump-up and jump-down phenomena are observed in second-order MMS version II, which closely agree with the experimental findings. The results are compared with those obtained by numerically integrating the original temporal equation.     

Nonlinear System Identification of Multi-Degree-of-Freedom Systems

A nonlinear system identification methodology based on theprinciple of harmonic balance is extended tomulti-degree-of-freedom systems. The methodology, called HarmonicBalance Nonlinearity IDentification (HBNID), is then used toidentify two theoretical two-degree-of-freedom models and anexperimental single-degree-of freedom system. The three modelsand experiments deal with self-excited motions of afluid-structure system with a subcritical Hopf bifurcation. Theperformance of HBNID in capturing the stable and unstable limitcycles in the global bifurcation behavior of these systems is alsostudied. It is found that if the model structure is well known,HBNID performs well in capturing the unknown parameters. If themodel structure is not well known, however, HBNID captures thestable limit cycle but not the unstable limit cycle.     

Free vibrations of a thin elastica by normal modes

In this work we investigate the existence, stability and bifurcation of periodic motions in an unforced conservative two degree of freedom system. The system models the nonlinear vibrations of an elastic rod which can undergo both torsional and bending modes. Using a variety of perturbation techniques in conjunction with the computer algebra system MACSYMA, we obtain approximate expressions for a diversity of periodic motions, including nonlinear normal modes, elliptic orbits and non-local modes. The latter motions, which involve both bending and torsional motions in a 2:1 ratio, correspond to behavior previously observed in experiments by Cusumano.     

强非线性振动系统周期解的能量迭代法

对于完全强非线性系统：x^.. g(x) f(x,x^.)x^.=0,提出求周期近似解析解以及这些解的稳定性的新方法。式中，g(x)、f(x,x^.)x^.分别是x,x、x^.的非线性函数。方法是基于能量原理，求出其一次近似解析解，然后引进牛顿迭代思想，得到周期系统数微分方程，最后根据谐波平衡原理及最小二乘法求其高次近似解，高次近似解的表达式由计算机辅助推导。计算参考文献[2]和[3]中的例题，令其中ε=1，研究该完全强非线性系统的周期解及其稳定性，本文方法与龙格-库塔数值法算得的结果对照如图1-3所示，它们表明本文方法不仅有效而且精度较高。     

Versal Deformation and Local Bifurcation Analysis of Time-Periodic Nonlinear Systems

In this study a local semi-analytical method of quantitativebifurcation analysis for time-periodic nonlinear systems is presented.In the neighborhood of a local bifurcation point the system equationsare simplified via Lyapunov–Floquet transformation whichtransforms the linear part of the equation into a dynamically equivalenttime-invariant form. Then the time-periodic center manifoldreduction is used to separate the `critical' states and reduce thedimension of the system to a possible minimum. The center manifoldequations can be simplified further via time-dependent normal formtheory. For most codimension one cases these nonlinear normal forms arecompletely time-invariant. Versal deformation of thesetime-invariant normal forms can be found and the bifurcation phenomenoncan be studied in the neighborhood of the critical point. However, ingeneral, it is not a trivial task to find a quantitatively correctversal deformation for time-periodic systems. In order to do so, onemust find a relationship between the bifurcation parameter of theoriginal time-periodic system and the versal deformation parameter ofthe time-invariant normal form. Essentially one needs to find theeigenvalues of the fundamental solution matrix of the time-periodicproblem in terms of the system parameters, which, in general, cannot bedone due to computational difficulties. In this work two ideas areproposed to achieve this goal. The eigenvalues of the fundamentalsolution matrix can be related to the versal deformation parameter bysensitivity analysis and an approximation of any desired order can beobtained. This idea requires a symbolic computational procedure whichcan be very time consuming in some cases. An alternative method issuggested for faster results in which a second or higher order curvefitting technique is used to find the relationship. Once thisrelationship is established, closed form post-bifurcation steady-statesolutions can be obtained for flip, symmetry breaking, transcritical andsecondary Hopf bifurcations. Unlike averaging and perturbation methods,the proposed technique is applicable at any bifurcation point in theparameter space. As physical examples, a simple and a double pendulumsubjected to periodic parametric excitation are considered. A simple twodegrees of freedom model is also studied and the results are comparedwith those obtained from the traditional averaging method. All resultsare verified by numerical integration. It is observed that the proposedtechnique yields results which are very close to the numericalsolutions, unlike the averaging method.     

线性蠕变问题的摄动边界元-有限元方法

本文将摄动、边界元、有限元方法结合起来,提出一种求解线性蠕变问题的新方法。该方法不采用一般增量法中在一个时段内各物理量保持不变或作线性变化的假设,加大了计算步长提高了精度。文中构造了边界元摄动格式,构造了包含钢筋在内的边界元有限元耦合摄动格式,并给出了满意的数值结果。