Uniform Expansions of Periodic Solutions for the Third Superharmonic Resonance

A new method of uniform expansions of periodic solutions to ordinary differential equations has recently been proposed to study quasi-harmonic processes in non-linear dynamical systems, in particular, when a small parameter of non-linearity is absent. The main idea of the method consists in using the ratio of the amplitudes of higher harmonics to the amplitude of the first harmonic of a periodic solution as a small formal parameter that appears due to descending the amplitudes of harmonics monotonically with increasing their number (this is the condition that the term quasi-harmonic implies). In this paper, the method is generalized for the third superharmonic resonance (when the first and the third harmonics become of the same magnitude) in a harmonically forced oscillator with arbitrary odd polynomial non-linearity.     

Symbolic Computation of Local Stability and Bifurcation Surfaces for Nonlinear Time-Periodic Systems

In this paper we present a spectral technique for building asymptotic expansions which describe periodic processes in conservative and self-excited systems without assuming the oscillations to be weakly nonlinear. The small parameter of the expansion is connected with the ratio of the amplitudes of higher than the first harmonics in contrast to the traditional parameter connected with weak nonlinearity. In the case of an oscillator with power nonlinearity the frequency of the main harmonic and the complex amplitudes of higher harmonics are computed as the expansions of either integer (for weakly nonlinear oscillations) or algebraic (for strong nonlinearity) functions of the complex amplitude of the first harmonic depending on the character of the initial conditions and the maximum power of the nonlinear term in the equation. In the simplest case of weakly nonlinear oscillations the complete asymptotic expansion is shown to be valid in the whole domain of the periodic motions of definite type until the separatrix is reached. The expressions for the first terms of the expansion for concrete examples coincide with the expressions obtained both with the use of other methods and by expanding the exact solutions. For some special cases of the strongly nonlinear oscillations the comparison of the results with known exact solutions is carried out as well as the criteria of convergence of the expansions are determined.     

Synthesis of oscillators for exact desired periodic solutions having a finite number of harmonics

The synthesis of autonomous oscillators with exact desired periodic steady-state solution is described in this contribution. The vector field of the oscillator differential equation is built up with a conservative and a dissipative part. Both parts are synthesized using an algebraic function describing the desired limit cycle. The desired periodic motion is restricted by a finite numbers of harmonics, whereby the amplitude and the phase shift of every harmonic can be freely chosen, depending on the specific application. Afterwards the synthesis of a periodically driven oscillator with an exact desired periodic response is described. For this purpose, the differential equation of the autonomous oscillator is extended by a state-dependent compensation term that equals the excitation at the steady-state solution. Here the freely definable amplitudes and phase angles of the oscillator motion are restricted by the existence and stability conditions for synchronization.     

Accurate Higher-Order Analytical Approximate Solutions to Large-Amplitude Oscillating Systems with a General Non-Rational Restoring Force

A new approximate analytical approach for accurate higher-order nonlinear solutions of oscillations with large amplitude is presented in this paper. The oscillatory system is subjected to a non-rational restoring force. This approach is built upon linearization of the governing dynamic equation associated with the method of harmonic balance. Unlike the classical harmonic balance method, simple linear algebraic equations instead of nonlinear algebraic equations are obtained upon linearization prior to harmonic balancing. This approach also explores large parameter regions beyond the classical perturbation methods which in principle are confined to problems with small parameters. It has significant contribution as there exist many nonlinear problems without small parameters. Through some examples in this paper, we establish the general approximate analytical formulas for the exact period and periodic solution which are valid for small as well as large amplitudes of oscillation.     

高阶谐波平衡方法中非物理解来源分析及改进方法研究

对于周期性非定常问题,高阶谐波平衡（High-order Harmonic Balance, HOHB）方法将非定常方程的解用Fourier 级数展开至一定阶次,从而消除其中的时间导数项,大大降低了计算消耗. 本文以达芬振子方程为例,探讨了HOHB 方法中非物理解的来源,分析结果表明:非物理解出现的原因是在推导过程中非线性项的简化处理导致方程左右两边并不严格相等. 根据非线性项的特点,在其处理过程中扩充子时间层上的时域解,并将非线性项中出现的更高阶谐波截断,使方程左右两边严格相等. 通过对达芬振子方程进行数值模拟发现:改进方法在消除非物理解的同时,也显著减少了计算所需谐波数. 对比参考文献发现,同阶改进方法的精度和原始谐波平衡方法基本相当,证明了本方法的可行性. 最后将本方法应用于具有立方刚度非线性的气动弹性系统中,验证本方法的工程适用性. 但是,当方程中非线性项较多时,本方法所需要的计算消耗会有所增加.      

Non-linear rheology of a face-centred cubic phase in a diblock copolymer gel

The viscoelastic behaviour of a poly(oxyethylene)-poly(oxybutylene) diblock copolymer in aqueous solution forming a face-centred cubic (fcc) micellar phase has been investigated using oscillatory shear rheometry. With increasing strain amplitude, the micellar solution was observed to undergo a transition from linear to non-linear behaviour, characterized by strong shear thinning. The non-linear behaviour observed in the stress response was analyzed by Fourier transformation of the waveform. Fourier analysis revealed that the high harmonic contributions to the shear stress response increased with strain amplitude and up to the 81st harmonic was observed for very large amplitudes. The onset of non-linear response as defined from the dependence of isochronal dynamic shear moduli on strain amplitude was found to be in good agreement with that defined by the appearance of a higher harmonic in the stress waveform. The amplitudes of the harmonic coefficients are compared to the predictions of a model for the nonlinear rheological response of a lyotropic cubic mesophase based on the stress response to a periodic lattice potential (Jones and McLeish 1995). It is found that the model is able to account for qualitative trends in the data such as the development of finite higher harmonics with increasing strain, but it does not describe the full frequency and strain dependence of these coefficients. Received: 31 May 2000 Accepted: 21 August 2000     

Modified Parker-Sochacki method for solving nonlinear oscillators

A new approach is presented for solving nonlinear oscillatory systems. Parker-Sochacki method (PSM) is combined with Laplace-Padé resummation method to obtain approximate periodic solutions for three nonlinear oscillators. The first one is Duffing oscillator with quintic nonlinearity which has odd nonlinearity. The second one is Helmholtz oscillator which has even nonlinearity. The last one is a strongly nonlinear oscillator, namely; relativistic harmonic oscillator which has a fractional order nonlinearity. Solutions are also obtained using Runge-Kutta numerical method (RKM) and Lindstedt-Poincare method (LPM). However, the LPM could not be used to solve the relativistic harmonic oscillator since it is a strongly nonlinear oscillator. The comparison between these solutions shows that the convergence zone for the Parker-Sochacki with Laplace-Padé method (PSLPM) is remarkably increased compared to PSM method. It also shows that the PSLPM solutions are in excellent agreement with LPM solutions for Duffing oscillator and are superior to LPM solutions in case of Helmholtz oscillator. The PSLPM succeeded to give an accurate periodic solution for the relativistic harmonic oscillator. For a wide range of solution domain, comparing PSLPM with RKM prove the correctness of the PSLPM method. Hence, the PSLPM method can be used with satisfied confidence to solve a broad class of nonlinear oscillators.     

准周期响应对称破缺分岔点的一种快速计算方法

稳态响应如周期及准周期解的分岔计算, 是非线性动力学研究的难点问题之一. 与计算方法及分析理论相对完善的周期响应相比, 准周期响应的求解只是在近些年才得到较大进展, 而且其分岔分析更加棘手, 仍需要更有效的理论和方法. 目前, 稳态响应尤其是准周期响应的分岔计算, 一般需采用数值方法, 通过调节参数反复试算得到. 为此, 本文基于增量谐波平衡IHB法提出一种快速方法, 可以高效地确定准周期响应的对称破缺分岔点. 方法的理论基础是在准周期解的广义谐波级数表达基础上, 当响应发生对称破缺分岔时, 其偶次(含零次)谐波系数将逐渐由0变为小量. 基于此性质, 将零次谐波系数预先设定为小量, 同时将分岔控制参数视为可变的迭代变量, 进而通过IHB法构造迭代格式. 作为算例, 研究不可约频率作用下的双频激励Duffing系统以及Duffing-van der Pol耦合系统. 结果表明, 只要迭代格式收敛, 随着预设小量减小, 控制参数将逐渐接近分岔近似值; 同时, 通过提高谐波截断数可显著提高近似分岔值的计算精度. 所提方法无需反复试算, 只要迭代过程收敛、便可实现分岔点直接快速计算.      

A stability in the large investigation for a non-linear second order two-degree-of-freedom system with periodic excitation

An extension to an algorithm due to Simpson has been developed for the analysis of a non-linear second order two-degree-of-freedom system with external periodic excitation. The form of equations considered arises from the study of mechanical systems with a single concentrated weak non-linearity and the method assumes a solution made up of harmonic terms whose amplitudes vary slowly in time. The system considered is such that in the absence of external excitation, it possesses a stable equilibrium point and an unstable limit cycle arising from a sub-critical Hopf bifurcation. When forcing is applied, the stable equilibrium point may then be replaced by a stable periodic attractor, and the limit cycle by an unstable multi-periodic attractor. The method has been applied to the problem of locating these attractors, and if they exist, of finding the stable attractor's basin of attraction in terms of initial conditions. The method reduces the problem from a search in four-dimensional phase space to a search for a boundary in a plane defined by amplitudes a₁ and a₂ in the assumed form of the solution.The method was applied to three non-linear systems in which the non-linearity was due to either a linear spring with a small amount of cubic hardening or a linear spring with freeplay. Agreement was shown to be good in those cases where the non-linearity was weak. However, the method would not be expected to give such accurate results if the non-linear effect was more significant. This was illustrated for a case involving the freeplay non-linearity.     

Developed steady wave motions of a liquid film falling along a vertical plane

The falling of a thin viscous fluid layer (film) along a vertical plane under the effect of gravity is accompanied by wave motions in which capillary forces play an essential part. An equation for the film thickness h(x, t) is used extensively in analyses of these motions. This equation, obtained from the Navier—Stokes equations and the boundary conditions under different assumptions, reduces to an ordinary third-order nonlinear differential equation [1–7] for steady plane motions. Periodic solutions of this equation were sought by the methods of asymptotic expansions in the amplitude or by Fourier series expansions [1–7], which assumes a sequential accounting of the nonlinearity as a small perturbation. This limits the validity of the results obtained to the domain of small amplitudes. The case of arbitrary amplitudes is considered in this paper. A solution of the problem, based on an asymptotic expansion in the parameter ε is constructed. In this expansion the equation for the first approximation remains nonlinear but admits of integration, which discloses the class of bounded periodic solutions. Moreover, strict integral relations (for any ε) are obtained, and a variational problem about seeking the lower bound of values of the mean film thickness and other characteristics of the ultimately developed optimal motions is formulated and solved on their basis. The results obtained agree with experiments.     

Primary resonance of multiple degree-of-freedom dynamic systems with strong non-linearity using the homotopy analysis method

A homotopy analysis method（HAM）is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existence of small parameters in the considered equation.The HAM provides a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter.Two examples are presented to show that the HAM solutions agree well with the results of the modified Linstedt-Poincar＇e method and the incremental harmonic balance method.     

Periodic and quasiperiodic motion for complex non-linear systems

A damped complex non-linear system corresponding to two coupled non-linear oscillators with a periodic damping force is investigated by an asymptotic perturbation method based on Fourier expansion and time rescaling. Four coupled equations for the amplitude and the phase of solutions are derived. Phase-locked solutions with period equal to the damping force period are possible only if the oscillators amplitudes are equal. On the contrary, if the oscillators amplitudes are different, periodic solutions exist only with a period different from the damping force period. These solutions are stable only for perturbations that conserve the phase difference and the square amplitude sum of the oscillators. Energy considerations are used in order to study existence and characteristics of quasiperiodic motion. We demonstrate that modulated motion can be also obtained for appropriate values of the detuning parameter and in this case an approximate analytic solution is easily constructed. If the detuning parameter decreases the modulation period increases and then diverges, an infinite-period bifurcation occurs and the resulting motion becomes unbounded. Analytic approximate solutions are checked by numerical integration.     

Bifurcations in a parametrically excited non-linear oscillator

A non-linear parametrically excited oscillator, that includes van der Pol as well as Duffing type non-linearities, is studied for its small non-linear motions using the method of averaging. The averaged equations, which form a dynamical system on the plane and depend on the linear damping and the detuning, are analyzed for their constant and periodic solutions. Bendixon's criterion is used to deduce the existence and the non-existence of limit cycle solutions for various values of the parameters. Then, using local bifurcation theory for “saddle-node”, pitchfork and “Hopf” bifurcations and some results from one and two parameter unfoldings of degenerate singularities, a partial bifurcation set is constructed. Since constant and periodic solutions of the averaged system correspond, respectively, to the periodic solutions and almost periodic or amplitude modulated motions of the original oscillator, the bifurcation set indicates some ways in which periodic solutions can become “entrained” or can break the entrainment for almost periodic oscillations.     

Oscillators with a quasi-constant restoring force: approximations for motion

In this work approximate solutions to conservative single-degree of freedom oscillators with a restoring force close to the one with a constant magnitude are derived. Approximate solutions are assumed as a truncated Fourier series and harmonic balancing is applied. In addition, the assumption that the response of the oscillators considered is close to the response of the antisymmetric oscillator is introduced. It is suggested in a novel way how to modify the differential equation of motion with the assumed solution so as to derive explicit expressions for the frequency and the amplitudes of harmonics in the first, second and third approximation are presented. The comparison of the results obtained with numerical solutions as well as with some existing approximate analytical results from the literature is also carried out, showing excellent accuracy.     

1/3 subharmonic solution of elliptical sandwich plates

IntroductionInrecentyears,withtheessentialadvantageoflightweightandhighrigidity ,sandwichplatesandshellshavebeenusedasanimportantpatternofstructuralelementsinaeronautical,astronauticalandnavalengineering .However,nonlinearproblemsforsandwichplatesandshellsareonlyinvestigatedbyafewbecauseofthedifficultiesofnonlinearmathematicalproblems.LiuRen_huaiandXuJia_chu[1,2 ]andothershavemadesomeinvestigationsinthisfield .Bifurcationofnonlinearvibrationforsandwichplateshasnotyetbeeninvestigated .Inthisp…     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton's method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

Analytical approximation to large-amplitude oscillation of a non-linear conservative system

This paper deals with non-linear oscillation of a conservative system having inertia and static non-linearities. By combining the linearization of the governing equation with the method of harmonic balance, we establish analytical approximate solutions for the non-linear oscillations of the system. Unlike the classical harmonic balance method, linearization is performed prior to proceeding with harmonic balancing, thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations. Hence, we are able to establish analytical approximate formulas for the exact frequency and periodic solution. These analytical approximate formulas show excellent agreement with the exact solutions, and are valid for small as well as large amplitudes of oscillation.     

Nonlinear analysis of chatter vibration in a cylindrical transverse grinding process with two time delays using?a?nonlinear time transformation method

A nonlinear dynamic system of cylindrical transverse grinding process is studied in this paper. The system consists of a grinding wheel and a workpiece, which are connected to the base by spring-damper elements, interacting with nonlinear normal forces. This two DOF model includes two time delays originated from the regenerative effects of the workpiece and the grinding wheel. Bifurcation points are located using a numerical algorithm by which we can find all the eigenvalues in a given rectangular region on the complex plane for the delayed differential equations. Supercritical bifurcation has been found for some sets of system parameter values. The amplitudes of the limit cycles are predicted using a nonlinear time transformation method, which is similar to the harmonic balance approach in that a periodic solution is approximated by a Fourier series. However, the main difference is that a nonlinear time ? is introduced in the Fourier series rather than the physical time t. The analytical solutions of stable limit cycles up to the third harmonics are compared with numerical simulations for the retarded system. It is shown that the proposed method gives accurate approximate solutions.     

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.     

The Response of an Oscillator Moving along a Parabola to an External Excitation

The behavior of a mass point moving along a parabola under theeffect of an external periodic excitation in resonance with the naturalfrequency of the oscillator is studied. The asymptotic perturbationmethod based on temporal rescaling and balancing of the harmonic termswith a simple iteration is used in order to determine the nonlinearmodulation equations for the amplitude and the phase of the oscillation.External force-response curves are shown and moreover jump phenomena arealso observed. In certain cases a second low frequency appears inaddition to the forcing frequency and then stable two-periodquasi-periodic motions are present with amplitudes depending on theinitial conditions. The value of the low frequency depends on theamplitude of the external excitation. A higher order perturbationanalysis is developed and the validity of the method is highlighted bycomparing the leading order and the higher order approximate analyticsolutions to numerical results.