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

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

强迫Van der Pol振子的动力学特性

采用增量谐波平衡方法导出强迫Van der Pol振子稳态周期响应的IHB计算格式．以外激励频率为参数进行跟踪延续获得了系统主共振时的幅频响应特性，并作出了特定系统参数下的周期响应极限环．其结果与Runge—Kutta方法进行了对比，结果表明该算法精度可以灵活控制，且收敛速度快，结果可靠，是非线性电路系统等工程应用中强非线性问题动力学特性分析的有效方法．     

一类非自治滞后-自激系统的主共振与锁模现象

研究一类受外激励作用的滞后－自激系统（van der Pol系统）的振动。对于主共振情况，用平均法求出了稳态响应方程，揭示了响应与系统参数的关系。结果表明，外激励的影响占主导地位，滞后因素对非线性共振频率有影响，而自激因素使得周期解在偏离非线性共振频率后发生Hopf分岔，转变为概周期运动。通过用映射法求得旋转数，揭示了主共振前、后的概周期运动中存在的锁模现象，其阶数分布符合Farey树规律。研究表明，锁模振动因滞后非线性因素而产生。随着滞后参数值（代表滞后程度）增大，各阶共振特别是亚谐振动的存在能力增强，对系统结构会带来较大影响。     

高速重载人字齿轮传动非线性动力学分析

人字齿轮承载能力强,重合度大,可靠性高,多于高速、重载工况下使用.探究高速重载人字齿轮传动系统非线性动力学特性,可为其设计提供参考.首先,计算齿轮副时变啮合刚度;引入齿侧间隙、间隙非线性函数和综合传动误差,计算时变啮合力;引入轴承游隙,计算轴承受力.随后,建立高速重载人字齿轮传动系统非线性动力学方程,使用4阶Runge-Kutta法对方程求解.最后,探究不同因素对系统动态响应影响.保持系统其他参数不变,分别改变输入转矩、啮合阻尼、齿侧间隙、啮合刚度及激励频率,绘制系统时间-位移图像、时间-速度图像、空间相图、空间频谱图及分岔图,观察系统非线性动力学响应变化趋势,判别系统运动状态.结果表明:在一定范围内,系统稳定性与啮合阻尼、啮合刚度呈正相关关系,与齿侧间隙、输入转矩呈负相关关系;逐渐增大外部激励频率时,系统运动从单周期运动逐渐变为混沌运动,随后又回归单周期运动.为保证系统平稳运行,应合理选取外部激励频率.     

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

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

压电复合材料层合梁的分岔、混沌动力学与控制

研究了简支压电复合材料层合梁在轴向、横向载荷共同作用下的非线性动力学、分岔和混沌动力学响应. 基于vonKarman理论和Reddy高阶剪切变形理论,推导出了压电复合层合梁的动力学方程. 利用Galerkin法离散偏微分方程,得到两个自由度非线性控制方程,并且利用多尺度法得到了平均方程. 基于平均方程,研究了压电层合梁系统的动态分岔,分析了系统各种参数对倍周期分岔的影响及变化规律. 结果表明,压电复合材料层合梁周期运动的稳定性和混沌运动对外激励的变化非常敏感,通过控制压电激励,可以控制压电复合材料层合梁的振动,保持系统的稳定性,即控制系统产生倍周期分岔解,从而阻止系统通过倍周期分岔进入混沌运动,并给出了控制分岔图.      

磁场中轴向运动条形板强非线性超谐共振及混沌运动

针对磁场环境中周期外载作用下轴向运动导电条形板的非线性振动及混沌运动问题进行研究。应用改进多尺度法对横向磁场中条形板的强非线性振动问题进行求解,得到超谐波共振下系统的分岔响应方程。根据奇异性理论对非线性动力学系统的普适开折进行分析,求得含两个开折参数的转迁集及对应区域的拓扑结构分岔图。通过数值算例,分别得到以磁感应强度、轴向拉力、激励力幅值和激励频率为分岔控制参数的分岔图和最大李雅普诺夫指数图,以及反映不同运动行为区域的动力学响应图形,讨论分岔参数对系统呈现的倍周期和混沌运动的影响。结果表明,可通过相应参数的改变实现对系统复杂动力学行为的控制。     

悬索在考虑1:3内共振情况下的动力学行为

研究了悬索在受到外激励作用下考虑1∶3内共振情况下的两模态非线性响应.对于一定范围内悬索的弹性-几何参数而言,悬索的第三阶面内对称模态的固有频率接近于第一阶面内对称模态固有频率的三倍,从而导致1∶3内共振的存在.首先利用Galerkin方法把悬索的面内运动方程进行离散,然后利用多尺度法对离散的运动方程进行摄动得到主共振情况下的平均方程.接下来对平均方程的稳态解、周期解以及混沌解进行了研究.最后利用Runge-Kutta法研究了悬索两自由度离散模型的非线性响应.     

非线性动力系统的频率近似法

给出非线性动力系统周期振动的频率近似法，本法将描述动力系统的非线性微分方程，化为以相角为自变量，振动频率为未知函数的积分方程，将弹性恢复力表示为线性及非线性两部分，从而得到积分方程的近似解，即频率的近似表达式。     

多尺度法二义性的一种解释

多尺度法是为解决含小参数系统发展起来的应用最广泛的摄动法之一. 在求解高阶近 似方程时，多尺度法一般只求特解. 用多尺度法求解van der Pol 方程的三阶解时 将出现矛盾. 以van der Pol方程为例，证明了忽略一阶修正量中的一阶谐波 项使得混合偏导数不能交换顺序，从而导致了多尺度法的二义性和另一个数学矛盾. 在求解一阶修正量时采用含有一阶谐波项的全解，消除了二义性和该矛盾. 该 方法所求得的近似解与数值解进行了比较，结果非常吻合，验证了其合理性.     

Effect of External Excitations on a Nonlinear System with Time Delay

The trivial equilibrium of a two-degree-of-freedom autonomous system may become unstable via a Hopf bifurcation of multiplicity two and give rise to oscillatory bifurcating solutions, due to presence of a time delay in the linear and nonlinear terms. The effect of external excitations on the dynamic behaviour of the corresponding non-autonomous system, after the Hopf bifurcation, is investigated based on the behaviour of solutions to the four-dimensional system of ordinary differential equations. The interaction between the Hopf bifurcating solutions and the high level excitations may induce a non-resonant or secondary resonance response, depending on the ratio of the frequency of bifurcating periodic motion to the frequency of external excitation. The first-order approximate periodic solutions for the non-resonant and super-harmonic resonance response are found to be in good agreement with those obtained by direct numerical integration of the delay differential equation. It is found that the non-resonant response may be either periodic or quasi-periodic. It is shown that the super-harmonic resonance response may exhibit periodic and quasi-periodic motions as well as a co-existence of two or three stable motions.     

A novel approach for generating multi-direction multi-double-scroll attractors

Attention is focused in this work on quasiperiodic motion of nonlinear systems whose spectrum contains uniformly spaced sideband frequencies with a distance \(\omega _{d}\) apart, around a frequency \(\omega \) with \(\omega \gg \omega _{d}\) and its integer multiples, which are referred to as carrier frequencies. The ratio of the two frequencies \(\omega \) and \(\omega _{d}\) is an irrational number. A new method based on the traditional incremental harmonic balance (IHB) method with multiple timescales, referred to as Lau method, where two timescales, \(\tau _{1}=\omega t\) (a fast timescale) and \(\tau _{2}=\omega _{d}t\) (a slow timescale), are introduced, is presented to analyze quasiperiodic motion of nonlinear systems. An amplitude increment algorithm is adapted to deal with cases where the two frequencies \(\omega \) and \(\omega _{d}\) are    unknown a priori, in order to automatically trace frequency response of quasiperiodic motion of nonlinear systems and accurately calculate all frequency components and their corresponding amplitudes. Results of application of the present IHB method to quasiperiodic free vibration of a hinged–clamped beam with internal resonance between two transverse modes are shown and compared with previously published results with Lau method and those from numerical integration. While differences are noted between results predicted by the present IHB method and Lau method, excellent agreement is achieved between results from the present IHB method and numerical integration even in cases of strongly nonlinear vibration. The present IHB method is also used to analyze quasiperiodic free vibration of high-dimensional models of the hinged–clamped beam.     

Subharmonic resonance of a clamped-clamped buckled beam with 1:1 internal resonance under base harmonic excitations

The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated. The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton's principle. A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method. A high-dimensional model of the buckled beam is derived, concerning nonlinear coupling. The incremental harmonic balance (IHB) method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve, and the Floquet theory is used to analyze the stability of the periodic solutions. Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited. Bifurcations including the saddle-node, Hopf, perioddoubling, and symmetry-breaking bifurcations are observed. Furthermore, quasi-periodic motion is observed by using the fourth-order Runge-Kutta method, which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.     

悬臂输流管道在基础激励与脉动内流联合作用下的参激振动

首先建立了悬臂输流管道在基础激励与脉动内流联合作用下的运动方程;然后基于Galerkin法研究了该系统的非线性动力学行为,分析了系统运动状态随激励频率和相位差的变化,以及混沌百分比随频率比(基础激励频率与脉动频率之比)和相位差的变化。结果表明,无论以激励频率还是以相位差为分岔参数,系统都具有多种形式的动态响应,包括周期运动、概周期运动和混沌运动,但进入和脱离混沌的途径不同。相位差和频率比对系统的混沌百分比有重要影响:相位差为π/2时系统混沌百分比最大;频率比为1时系统混沌百分比最小,频率比较小或较大时系统混沌百分比与只有基础激励时接近。     

Energy harvesting in a Mathieu–van der Pol–Duffing MEMS device using time delay

This paper investigates quasi-periodic vibration-based energy harvesting in a delayed nonlinear MEMS device consisting of a delayed Mathieu–van der Pol–Duffing type oscillator coupled to a delayed piezoelectric coupling mechanism. We use the multiple scales method to approximate the quasi-periodic response and the related power output near the principal parametric resonance. The effect of time delay on the energy harvesting performance is studied. It is shown that for appropriate combination of time delay parameters, there exists an optimum range of excitation frequency beyond the resonance where quasi-periodic vibration-based energy harvesting is maximum. Numerical simulations are performed to confirm the analytical predictions.     

Simultaneous combination and 1:3:5 internal resonances in a parametrically excited beam-mass system

The nonlinear dynamics of a base-excited slender beam carrying a lumped mass subjected to simultaneous combination parametric resonance of sum and difference type along with 1:3:5 internal resonances is investigated. Method of normal form is applied to the governing nonlinear temporal differential equation of motion to obtain a set of first-order differential equations which are used to obtain the steady-state, periodic, quasi-periodic and chaotic responses for different control parameters viz., amplitude and frequency of external excitation and damping. Frequency response, phase portraits, time spectra and bifurcation diagram are plotted to visualize the system behaviour with variation in the control parameters. Here, two distinct zones of trivial instability, blue sky catastrophe phenomena, jump down phenomena, simultaneous occurrence of periodic and chaotic orbits, period doubling of the mixed-mode periodic orbits leading to chaos, attractor merging crisis, boundary crisis, type II and on-off intermittencies are observed. Bifurcation diagram is plotted to facilitate the designer to choose a safe operating zone.     

Dynamics of two strongly coupled van der pol oscillators

A perturbation method is used to study the steady state behavior of two Van der Pol oscillators with strong linear diffusive coupling. It is shown that a bifurcation occurs which results in a transition from phase-locked periodic motions to quasi-periodic motions as the coupling is decreased or the detuning is increased. The analytical results are compared with a numerically generated solution.