具任意脱层复合材料梁的非线性谐波响应

研究了具任意脱层复合材料梁的非线性谐波响应问题。基于弹性理论，建立了考虑剪切变形时的复合材料梁脱层的基本方程式。在空间上采用B样条函数和Galerkin积分法，在时间上采用增量谐波平衡法进行计算。通过实例计算，得出了简谐力作用下的非线性动力响应曲线。认为基谐波振动仍是非线性振动的主要部分。     

基于增量谐波平衡的参激系统非线性识别法

将增量谐波平衡法应用到非线性系统的建模和参数识别中,针对Mathieu-Duffing方程,推导了利用增量谐波平衡原理识别参数激励非线性系统参数的方法. 该方法改进了增量谐波平衡方法的推导过程,通过数值模拟对比研究了谐波平衡非线性识别(harmonic balance nonlinearity identification, HBNID)和增量谐波平衡非线性识别(incremental harmonic balance nonlinearity identification, IHBNID)的效果,验证了增量谐波平衡非线性识别的有效性. 结果表明,增量谐波平衡非线性识别的计算效率较高,计算精度和抗噪能力都优于谐波平衡非线性识别.      

直升机桨叶液压减摆器等效阻尼计算研究

直升机桨叶液压减摆器具有强非线性、大阻尼的特点。为了正确评估减摆器的等效阻尼，就旋转桨叶在水平面内的自由摆振运动，采用基于富里叶分析的移动矩形窗、基于希尔贝特（Hilbert）变换和能量平衡的方法对减摆器的等效阻尼进行了计算研究，并就阻尼强弱、衰减对数曲线的拟合等因素对计算结果的影响进行了分析。结果表明，对于线性阻尼和非线性弱阻尼振动系统，包络衰减法和能量平衡法计算等效阻尼能获得较一致的结果，对于非线性阻尼振动系统来说结果有分散性，应用时须谨慎。     

一种单元谐波平衡法

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

非惯性参考系中弹性壳的非线性振动分析

给出了弹性壳处于非惯性参系中的运动描述，基于Ｈａｍｉｌｔｏｎ原理建立了中厚壳在非惯性参考系中的非线性运动控制方程，应用多尺度法及谐波平衡法具体地分析了圆柱壳的非线性振动问题     

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

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

粘弹性结构动力稳定性分析的谐波平衡法

本文分析了线粘弹性结构的长期动力稳定特性。设材料具积分型本构关系，且其松弛模量能用Prony级数描述，将微分-积分型控制方程化成微分型方程，应用谐波平衡法确定动力稳定区域，着重讨论了材料参数及系统振动频率对动力稳定区域的影响发现该类粘弹性结构具有与一般阻尼系统不同的动力稳定特性。文中也将系统平衡法直接应用于微分-积分型控制方程，忽略卷积积分运算后产生的随时间衰减的非谐波项来得到决定动力稳定边界的特征方程，并对两种应用所得结果进行了比较。     

含孔von Kármán板中的高次谐波散射现象

基于vonKarman板大挠度弯曲理论，利用谐波平衡法及小参数振动法，研究了含孔vonKarman板的非线性波散射问题．通过分析发现：由于弯曲应力与中面力的非线性耦合，vonKarman板中会出现高次谐波散射现象．     

非惯性参考系中板的非线性振动分析

给出了弹性变形体处于非惯性参考系中的运动描述，基于Ｈａｍｉｌｔｏｎ原理建立了中厚矩形板在空间运动中的动力学一般方程，应用多尺度法及谐波平衡法具体地分析了简支矩形板在空间运动中的非线性振动问题。     

二次非线性粘弹性圆板的2/1+3/1超谐解

计及材料的非线性弹性和粘性性质，研究了圆板在简谐载荷作用下的2／1＋3／1超谐解，导出了相应的非线性动力方程。提出一类强非线性动力系统的叠加－叠代谐波平衡法。将描述动力系统的二阶常微分方程，化为基本解为未知函数的基本微分方程；及分岔解为未知函数的增量微分方程。通过叠加－迭代谐波平衡法得出了圆板的2／1＋3／1超谐解。对叠加迭代谐波平衡法和数值积分法进行了比较，两者结果吻合很高。并且讨论了2／1＋3／1超谐解的渐近稳定性。     

Brownian dynamics simulation of rodlike polymers under shear flow

The highly nonlinear behaviors of rodlike polymers in nematic phase under shear flow are studied with Brownian dynamics simulation. The LebwohlLasher nematogen model is taken as the prototype of the simulation and the mean-field approximation is avoided. By considering the nearest-neighbor intermolecular interaction, the spatial orientational correlation is introduced and therefore the spatial inhomogeneity such as the multiple-domain effect can automatically be incorporated. The transient order parameters, birefringence axes, shear stresses and first normal stress differences are calculated. The important finding of this work is that the director wagging and damped oscillation share the same molecular origin as director tumbling. The only difference is that the system is split into micro-domains which tumble with different phase angles in the wagging and damped oscillation regimes. The tumbling of the director of the whole system is suppressed due to the spatial inhomogeneity of director fields and then the damped oscillation of macroscopic stresses becomes predominant. The negative first normal stress difference exists at moderate shear rates, where both elasticity and viscosity play important role. Our simulation results including some dimensionless scaling parameters find good agreement with experimental observations in literature.     

强非线性自激振动问题的瞬态解

用新的改进的L-P法并结合权残法求得了一般情况下(阻尼项可以不是小量)强非线性自激振动问题的一级近似瞬态解。近似解有很好的精度。     

A return toward equilibrium in a 2D Rayleigh–Taylor instability for compressible fluids with a multidomain adaptive Chebyshev method

We numerically simulate a single-mode Rayleigh–Taylor instability between compressible miscible fluids with a highly accurate self-adaptive pseudospectral Chebyshev multidomain method in two two-dimensional boxes at small aspect ratios. The simulations are started from rest and pursued until the return toward mechanical equilibrium of the mixing. Four regimes—linear and weakly nonlinear, nonlinear steady bubble rise, return toward equilibrium, and finally a system of acoustic waves—can be identified. We show that this one-dimensional system of stationary acoustic waves is damped by the physical viscosity. This provides a reference solution.      

Bifurcation analysis for nonlinear multi-degree-of-freedom rotor system with liquid-film lubricated bearings

The oil-film oscillation in a large rotating machinery is a complex high-dimensional nonlinear problem. In this paper, a high pressure rotor of an aero engine with a pair of liquid-film lubricated bearings is modeled as a twenty-two-degree-of-freedom nonlinear system by the Lagrange method. This high-dimensional nonlinear system can be reduced to a two-degree-of-freedom system preserving the oil-film oscillation property by introducing the modified proper orthogonal decomposition (POD) method. The efficiency of the method is shown by numerical simulations for both the original and reduced systems. The Chen-Longford (C-L) method is introduced to get the dynamical behaviors of the reduced system that reflect the natural property of the oil-film oscillation.     

Post-divergence and post-flutter analysis of sub-tangentially loaded and damped Beck column using a nonlinear FE procedure

Both post-divergence and post-flutter behaviors of damped Beck columns subjected to a sub-tangentially follower force are rigorously explored using an exact co-rotational frame element. First a linear stability theory of the damped Beck column is summarized using a stability map. A geometrically nonlinear frame element based on the co-rotational formulation is then formulated including mass matrix, Rayleigh damping matrix, and load-correction stiffness matrix due to circulatory forces. The dynamic FE analysis is performed using Newmark integration method. Finally a Beck column model is parametrically analyzed in order to investigate non-linear stability characteristics of the internally and externally damped non-conservative system. In particular, interesting nonlinear behaviors of Beck column quite different from those predicted by the linear theory are reported through static and dynamic nonlinear analysis with variation of sub-tangentiality of the follower force.     

STUDY ABUT NONLINEAR FACTOR INFLUENCE ON THE VIBRATION ISOLATION EFFECT OF DAMPER

本文在同时考虑刚度和阻尼非线性的基础上, 研究了非线性因素对隔振器隔振效果的影响;获得了同时考虑立方刚度非线性和椭圆阻尼非线性的非线性系统运动响应和传递率计算公式;讨论了非线性因素对系统运动响应和传递率的影响.     

Relaxation oscillations,subharmonic orbits and chaos in the dynamics of a linear lattice with a local essentially nonlinear attachment

We study the dynamic interactions between traveling waves propagating in a linear lattice and a lightweight, essentially nonlinear and damped local attachment. Correct to leading order, we reduce the dynamics to a strongly nonlinear damped oscillator forced by two harmonic terms. One of the excitation frequencies is characteristic of the traveling wave that impedes to the attachment, whereas the other accounts for local lattice dynamics. These two frequencies are energy-independent; a third energy-dependent frequency is present in the problem, characterizing the nonlinear oscillation of the attachment when forced by the traveling wave. We study this three-frequency strongly nonlinear problem through slow-fast partitions of the dynamics and resort to action-angle coordinates and Melnikov analysis. For damping below a critical threshold, we prove the existence of relaxation oscillations of the attachment; these oscillations are associated with enhanced targeted energy transfer from the traveling wave to the attachment. Moreover, in the limit of weak or no damping, we prove the existence of subharmonic oscillations of arbitrarily large periods, and of chaotic motions. The analytical results are supported by numerical simulations of the reduced order model.     

Connection between Volterra series and perturbation method in nonlinear systems analyses

This paper is concerned with the connection between the Volterra series and the regular perturbation method in nonlinear systems analyses. It is revealed for the first time that, for a forced polynomial nonlinear system, if its derived linear system is a damped dissipative system, the steady response obtained through the regular perturbation method is exactly identical to the response given by the Volterra series. On the other hand, if the derived linear system is an undamped conservative system, then the Volterra series is incapable of modeling the forced polynomial nonlinear system. Numerical examples are further presented to illustrate these points. The results provide a new criterion for quickly judging whether the Volterra series is applicable for modeling a given polynomial nonlinear system.     

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

粘弹性阻尼修改后结构系统动特性的预测

无阻尼或比例阻尼结构系统经粘弹阻尼修改后,可变为粘弹性阻尼系统。要获得其动特性,需求解复特征值问题。但是,随之带来了计算量大、费用高等问题。尤其是粘弹性材料特性随频率而变化,需求解高阶非线性复特征值问题,这对于一个自由度较大的结构,计算量太大,通常较难实现。本文在特征值修改方法的基础上,提出粘弹性阻尼局部动力修改方法,即仅需已知原结构系统的实模态参数,就可求出粘弹阻尼修改后系统的复模态参数。还发展了特征值和修改量同时迭代方法,有效地解决了粘弹材料复模量随频率变化引起的非线性复特征值问题。