高维非线性振动系统参数识别

将增量谐波平衡非线性识别推广到高维振动系统, 推导了基于增量谐波平衡的多自由度非线性系统的识别方程. 针对一个两自由度系统进行了数值模拟计算, 讨论了系统在单周期、倍周期和混沌运动状态下的参数识别, 以及噪声对识别结果的影响, 验证了增量谐波平衡非线性识别在多自由度系统中的有效性. 结果表明, 该方法具有较高的计算效率和识别精度, 以及良好的抗噪能力.      

黏弹性复合材料梁动力学实验建模的研究

以ABS树脂为基材, 填充1%~10%的金红石纳米二氧化钛制成纳米复合材料样本系列,搭建了参数激励非线性振动实验系统. 采用实验建模的方法, 基于非线性增量谐波平衡识别理论,建立了黏弹性复合材料屈曲梁的动力学控制方程. 通过数值模拟与实验结果的比较, 验证了理论模型和实验系统在定性定量分析上的一致性, 并且对一类不同配比成分的纳米复合材料也有很好的适用性.      

参强联合作用非线性结构动力学实验建模

搭建以L 型梁为实验研究对象的参强联合作用多自由度非线性振动实验系统, 将增量谐波平衡非线性识别理论运用到实验建模方法中, 建立了L 型梁的动力学控制方程. 通过对不同激励频率和不同响应情况下的数值模拟与实验数据的比较, 验证了基于增量谐波平衡识别的实验建模方法对多自由度参强联合作用非线性动力学结构的有效性, 以及动力学控制方程的普适性.     

应用于具有二次,三次非线性系统的增量谐波平衡法

本文导出了适用于具有二次、三次非线性的微分方程组的增量谐波平衡法，研究了扁拱的相加型和相减型的联合共振问题以及二自由度系统的强非线性振动问题，算例表明，增量谐波平衡法是一个求解多自由度系统强非线性振动的有效的半解析的数值方法。     

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

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

谐波平衡法与复平均法计算强非线性系统稳态响应的区别

近些年,很多学者致力于利用非线性增强振动响应减少的效果或者能量采集器的效率。因而非线性系统的响应值需要从理论计算方面更准确地预测。另外,根据学者已取得的研究成就,非线性能量汇(NES)中存在的立方刚度非线性可以将结构中宽频域的振动能量传递至非线性振子部分。文章将一种由NES和压电能量采集器组成的NES-piezo装置与两自由度主结构耦合连接,系统受谐和激励作用。文章采用谐波平衡法和复平均法分别推导了系统稳态响应,参照数值结果,对比两种近似解析方法在求解强非线性系统稳态响应时的异同。计算结果表明,系统体现较弱非线性时,二者计算结果差异很小;当系统体现强非线性时,复平均法不能准确地呈现系统高阶响应,提高阶数的谐波平衡法能更准确地表示系统响应值。基于谐波平衡法和数值算法,讨论NES-piezo装置对于系统宽频域减振的影响。与仅加入非线性能量汇情况对比,结果表明NES-piezo装置不会恶化宽频域减振效果,并且在第一阶共振频率附近,可以稍微提高结构减振效率。另外,计算结果也表明,采用恰当的NES-piezo装置可实现宽频域范围的结构减振和压电能量采集一体化。此项研究工作为研究不同情形强非线性系统的响应提供了理论方法的指导。另外,研究结果也为宽频域范围的结构减振和压电能量采集一体化提供了理论依据。     

非线性时滞反馈控制系统超谐共振分析

采用增量谐波平衡法求解了非线性时滞微分方程的超谐共振解,研究了时滞、反馈控制增益、激励幅值、非线性项系数等系统参数对系统超谐共振响应的影响,分析了超谐共振响应随系统参数变化的规律。结果表明:三次谐波与一次谐波振幅的比值随时滞量呈周期性变化;反馈控制增益对系统超谐共振的影响与非线性项系数和激励幅值有关;随着非线性项系数和激励幅值的不断增大,三次谐波项与一次谐波项振幅的比值都是先增大后减小,而且减小的趋势逐渐减弱;一次谐波成份在振幅中占主导地位。     

非线性转子-轴承系统的周期解及近似解析表达式

通过对普通打靶方法进行改造提出一种确定非线性系统周期轨道及周期的新型打靶算法。首先通过改变系统的时间尺度，将非线性系统周期轨道的周期显式地出现在非线性系统的系统方程中，然后对传统打靶法进行改造，将周期也作为一个参数一起参与打靶法的迭代过程，迭代过程包含对周期轨道和周期的求解，迭代过程中的增量通过优化方法选择，从而能迅速确定出系统的周期轨道及其周期。应用所求的结果结合谐波平衡方法求得了非线性系统的周期轨道的近似解析表达式，理论上通过增加谐波的阶数任何精度的周期解都可以得到。最后将该方法应用于非线性转子轴承系统，求出了在某些参数下转子的周期解及其近似解析表达式，通过与四阶Runge-Kutta数值积分结果比较，验证了方法的有效性，计算结果对于转子系统运动的定量控制有重要理论指导意义。     

一种单元谐波平衡法

基于有限元离散,对于工程中的非线性响应问题,提出一种单元谐波平衡法．与常规的谐波平衡法不同,本文将谐波平衡方程建立在有限元素上,从而兼顾了有限元素法和常规谐波平衡法两大优势．有限元技术的应用能使得求解问题的范围扩大到复杂工程结构,而谐波平衡概念的使用将使得含有复杂变形和复杂本构关系的动力学响应问题得到有效解决．所提方法能适用于工程结构中具有复杂非线性关系的动力学响应问题．由于谐波平衡法的实施依赖于谐波系数方程及其切线刚度矩阵的解析推导,尽管已经局限到有限元素上,但对于较为复杂一些的本构关系,推导仍非易事．为解决这些问题,放弃通常对于变形梯度和应变张量所作的向量假设,而是从连续介质力学中基本的几何关系入手,提出一种矩阵分解形式．通过利用张量的内蕴导数定义以及关于迹函数的有关性质,给出应力增量的一种新的表现形式．当它与变形梯度的矩阵分解相结合时,使得切线刚度矩阵的导出变得十分简单,而且所得计算形式也比通常紧凑和方便许多．     

用IHB法分析双谐波激励下铰接塔-油轮系统的非线性动力学特性

研究了考虑平方阻尼情况下,铰接塔-油轮系统在双谐波激励下的非线性动力学特性.将该系统简化为单自由度分段线性恢复力,含平方阻尼的运动学分析模型,建立了铰接装载塔系统的分段非线性动力学方程.采用增量谐波平衡法获得系统周期解,使用Floquet理论判断系统的运动稳定性,结合路径跟踪法跟踪系统响应曲线,获得了系统所有可能的亚谐、谐波、组合谐波共振运动.分析了不对称恢复刚度比值对系统亚谐、组合谐波共振和对系统运动倍周期分岔点的影响,比较了考虑平方阻尼和不考虑平方阻尼情况下系统非线性动力学特性,得到了系统的一些重要的非线性动力学特点.     

Supercritical as well as subcritical Hopf bifurcation in nonlinear flutter systems

The Hopf bifurcations of an airfoil flutter system with a cubic nonlinearity are investigated, with the flow speed as the bifurcation parameter. The center manifold theory and complex normal form method are Used to obtain the bifurcation equation. Interestingly, for a certain linear pitching stiffness the Hopf bifurcation is both supercritical and subcritical. It is found, mathematically, this is caused by the fact that one coefficient in the bifurcation equation does not contain the first power of the bifurcation parameter. The solutions of the bifurcation equation are validated by the equivalent linearization method and incremental harmonic balance method.     

Nonlinear dynamic behavior of a bio-inspired embedded X-shaped vibration isolation system

To improve the vibration isolation performance and bandwidth, loading capacity and supporting stability of passive vibration isolation system by utilizing nonlinearity, a bio-inspired embedded X-shaped vibration isolation (BIE-XVI) structure is proposed considering muscle/tendon contractile functions, joint rotational friction and connecting rod mass simultaneously. Furthermore, the dynamic model with pure linear elements and geometric relationship are established and the nonlinear variation properties are investigated. The effects of the key parameters of the BIE-XVI structure on frequency response characteristics and vibration isolation range are analyzed thoroughly by incremental harmonic balance method in various working conditions. From the parametric investigations, it can be found that the sensitivities of the nonlinear resonance properties are markedly different with respect to the different structure parameters. For longer rod length, larger assembly angle and higher stiffnesses, the hardening nonlinearity is weakened, but the resonance peak does not necessarily decrease. Besides, the softening nonlinearity and hardening nonlinearity can be interconverted with changing isolated mass and excitation amplitude. The BIE-XVI structure can widen the isolation frequency range and reduce the resonance peak to improve the vibration isolation properties by adjusting/designing the structural parameters, which could realize quasi-zero-stiffness property for vibration isolation.

     

Limit cycle oscillations of rectangular cantilever wings containing cubic nonlinearity in an incompressible flow

Limit cycle oscillations (LCO) as well as nonlinear aeroelastic analysis of rectangular cantilever wings with a cubic nonlinearity are investigated. Aeroelastic equations of a rectangular cantilever wing with two degrees of freedom in an incompressible potential flow are presented in the time domain. The harmonic balance method is modified to calculate the LCO frequency and amplitude for rectangular wings. In order to verify the derived formulation, flutter boundaries are obtained via a linear analysis of the derived system of equations for five different cases and compared with experimental data. Satisfactory results are gained through this comparison. The problem of finding the LCO frequency and amplitude is solved via applying the two methods discussed for two different cases with hardening cubic nonlinearities. The results from first-, third- and fifth-order harmonic balance methods are compared with the results of an exact numerical solution. A close agreement is obtained between these harmonic balance methods and the exact numerical solution of the governing aeroelastic equations. Finally, the nonlinear aeroelastic analysis of a rectangular cantilever wing with a softening nonlinearity is studied.     

INCREMENTAL HARMONIC BALANCE METHOD FOR AIRFOIL FLUTTER WITH MULTIPLE STRONG NONLINEARITIES

The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.     

谐波平衡法在动导数快速预测中的应用研究

谐波平衡法以傅里叶级数展开为基础,将周期性非定常流场的非定常求解过程转化为几个定常流场的耦合求解过程,并通过重建得到整个流场的非定常过程. 建立了基于谐波平衡法的动导数快速预测方法,数值模拟了超声速带翼导弹俯仰的动态流场,并通过积分法获取了俯仰动导数,与实验结果吻合很好;且在同等计算精度下,谐波平衡法的计算效率是双时间步方法的13 倍. 应用谐波平衡法研究了较大范围内减缩频率对俯仰动导数的影响规律. 研究发现,对于本外形,当减缩频率降低到一定值后,俯仰动导数的值迅速变化,甚至发生变号;对此现象产生的原因进行了深入分析,并通过对导弹自激俯仰运动的数值模拟验证了该结果. 此外,针对大攻角条件下动态流场非线性强的特点,开展了谐波平衡法在大攻角下的适用性研究. 结果表明,谐波平衡法在大攻角下也能取得很好的计算结果.      

Averaging Oscillations with Small Fractional Damping and Delayed Terms

We demonstrate the method of averaging for conservative oscillators which may be strongly nonlinear, under small perturbations including delayed and/or fractional derivative terms. The unperturbed systems studied here include a harmonic oscillator, a strongly nonlinear oscillator with a cubic nonlinearity, as well as one with a nonanalytic nonlinearity. For the latter two cases, we use an approximate realization of the asymptotic method of averaging, based on harmonic balance. The averaged dynamics closely match the full numerical solutions in all cases, verifying the validity of the averaging procedure as well as the harmonic balance approximations therein. Moreover, interesting dynamics is uncovered in the strongly nonlinear case with small delayed terms, where arbitrarily many stable and unstable limit cycles can coexist, and infinitely many simultaneous saddle-node bifurcations can occur.     

Improving convergence of incremental harmonic balance method using homotopy analysis method

We have deduced incremental harmonic balance an iteration scheme in the （IHB） method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.     

Dynamic stability of a nonlinear multiple-nanobeam system

Nonlinear Dynamics - We use the incremental harmonic balance (IHB) method to analyse the dynamic stability problem of a nonlinear multiple-nanobeam system (MNBS) within the framework of...