轴向移动局部浸液单向板的1:3内共振分析

考虑单向板的轴向速度、轴向张力、流固耦合作用以及阻尼等因素, 基于由 von Kármán薄板大挠度方程得到的轴向移动局部浸液单向板的非线性振动方程, 研究了外激励作用下单向板在1:3内共振情况时的非线性振动特性. 首先利用Galerkin法对非线性振动方程离散化, 然后分别应用数值法和近似解析法对离散后模态方程组进行求解, 获得了系统内共振情况下复杂的幅频特性曲线, 并讨论了周期解的稳定性. 最后研究了1:3内共振系统平均方程组的运动分岔现象.     

基于平均法的旋转薄壁圆柱壳非线性行波振动研究

应用Donnell's简化壳理论,在考虑阻尼和几何非线性的情况下,基于平均法对旋转的薄壁悬臂圆柱壳在法向激振力作用下的非线性行波振动进行了研究.在分析过程中,首先,引入考虑阻尼及几何非线性的薄壁圆柱壳非线性波动方程,进行降阶处理后,得到模态坐标下的振动方程;其次,对模态方程进行平均化处理,确定转换矩阵,进行变量的幅值相角化,从而得到自治的标准化方程组;最后,由系统谐波共振周期解对应平均方程稳态解的原理,得到幅频特性方程.根据上述所得结果,进行了系统参数振动及稳定性研究,并进一步将结果与谐波平衡法及数值解作了比较.     

轴向运动复合材料圆柱壳的非线性振动研究

采用Runge–Kutta法和多尺度法对轴向运动分层复合材料薄壁圆柱壳的非线性振动特性进行了研究。首先根据层合壳理论建立轴向运动分层复合材料薄壁圆柱壳的波动方程,利用Galerkin法对方程进行离散,得到相互耦合模态方程组。然后应用Runge –Kutta法分析了不同参数条件下的幅频特性曲线,得到了系统由于固有频率接近所导致的内共振现象,以及系统呈现软特性等非线性特性。最后采用多尺度法进行了系统1:1内共振时的近似解析分析,对系统在不同参数下的振动研究表明,激振力幅值、阻尼、速度等参数对位移响应幅值、共振区间、模态间的耦合度及系统软特性程度均有影响,其结论与数值计算结果一致,并同时对解的稳定性进行了研究。     

悬索在外激励作用下的1:3内共振分析(Ⅱ):数值结果

本研究的第一部分已经推导了悬索在第一阶面内对称模态主共振和第三阶面内对称模态主共振下的平均方程,其中考虑了这两阶模态之间1∶3内共振.本文对平均方程的稳态解,周期解以及混沌解进行了研究.利用 Newton-Naphson 方法和拟弧长的延拓算法确定了主共振情况下的幅频响应曲线,通过利用 Jacobian 矩阵的特征值判断幅频响应曲线中解的稳定性.在这些幅频响应曲线中.都存在超临界 Hopf 分叉,导致平均方程的周期解.以这些超临界 Hopf 分叉为起点.利用打靶法和拟弧长的延拓算法确定了两种主共振情况下的周期解分支,同时通过利用 Floquet 理论判断这些周期解的稳定性.然后利用数值结果研究了两种主共振情况下的厨期解经过倍周期分叉通向混沌的过程.最后利用 Runge-Kutta 法研究了悬索两自由度离散模型的非线性响应.     

内共振条件下轴向运动梁磁弹性超谐共振研究

以载流导线激发的磁场中轴向运动梁为研究对象,同时考虑外激励力作用,推导出梁的磁弹性非线性振动方程．通过位移函数的设定和伽辽金积分法,将非线性振动方程离散为常微分方程组．采用多尺度法得到系统的近似解析解．应用Matlab 和Mathematica 软件求解幅频响应方程,并对稳态解进行稳定性判定．通过具体算例得到前两阶假设模态的响应幅值随不同参数的变化规律．结果发现：系统在内共振条件下发生超谐波共振时,二阶假设模态幅值明显小于一阶;随着外激励的增大,多值解区域范围明显缩小;随着电流强度增加,振动幅值减小,表明载流导线能够起到控制共振的作用．     

基于多尺度方法的平动圆柱贮箱航天器刚–液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚–液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚–液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励–幅值响应曲线.结果表明,随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚–液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚–液耦合非线性动力学特性具有重要参考价值.     

悬索在外激励作用下的1∶3内共振分析(II)∶数值结果

本研究的第一部分已经推导了悬索在第一阶面内对称模态主共振和第三阶面内对称模态主共振下的平均方程,其中考虑了这两阶模态之间1∶3内共振。本文对平均方程的稳态解、周期解以及混沌解进行了研究。利用Newton-Naphson方法和拟弧长的延拓算法确定了主共振情况下的幅频响应曲线,通过利用Jacobian矩阵的特征值判断幅频响应曲线中解的稳定性。在这些幅频响应曲线中,都存在超临界Hopf分叉,导致平均方程的周期解。以这些超临界Hopf分叉为起点,利用打靶法和拟弧长的延拓算法确定了两种主共振情况下的周期解分支,同时通过利用Floquet理论判断这些周期解的稳定性。然后利用数值结果研究了两种主共振情况下的周期解经过倍周期分叉通向混沌的过程。最后利用Runge-Kutta法研究了悬索两自由度离散模型的非线性响应。     

多频激励下悬索非线性共振特性受温度影响分析

推导了考虑温度变化影响的悬索非线性运动微分方程,利用Galerkin法得到离散后的多自由度方程;考虑一阶正对称模态,以悬索同时发生主共振和1/3阶次谐波共振为例,利用多尺度法求解幅频响应方程组,并判断稳态解的稳定性;选取三组垂跨比及两组温度变化,基于幅频响应曲线和调谐相位曲线,探究温度变化影响下的主/次谐波联合共振响应。数值算例结果表明:主/次谐波联合共振时,系统响应变得更加复杂,同时展现出主共振和次谐波共振响应特性;温度变化会定性和定量地改变联合共振特性,改变系统振动的软/硬弹簧特性及程度;联合共振响应的幅值大小、相位和共振区间与温度变化密切相关;相同温度变化对联合共振响应的幅值和相位影响有差异,通过研究联合共振响应的相位,可以区分系统的多个稳态解。     

基于多尺度方法的平动圆柱贮箱航天器刚——液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.      

基于分数阶磁流变液阻尼器模型的车辆悬架组合控制

磁流变液阻尼器的分数阶Bingham模型结构形式简单, 而且可以更好地描述系统的滞回特性. 建立了含有分数阶Bingham模型的单自由度1/4车辆悬架系统模型, 利用磁流变液阻尼器对在路面简谐激励下的非线性车辆悬架系统进行振动控制. 研究了含有分数阶Bingham模型的悬架系统在天棚阻尼半主动控制下的主共振响应, 利用平均法得到了系统的近似解析解. 求解了系统定常解的幅频响应方程, 并根据李雅普诺夫稳定性理论得到了悬架系统的稳定性条件. 通过绘制数值解和解析解的幅频响应曲线对比图, 验证了近似解析解的正确性. 利用簧载质量垂直方向的加速度均方根值分析了半主动控制对车辆乘坐舒适性的影响, 发现天棚阻尼半主动控制策略在低频激励区域反而会降低车辆的乘坐舒适性. 因此提出了一种被动控制与半主动控制相结合的组合控制策略, 并分析了半主动控制参数对振动控制效果的影响. 分析结果表明, 该组合控制策略不但能够提高车辆的乘坐舒适性, 而且能有效抑制悬架系统的主共振振动幅值.      

Nonlinear Random Response of Ocean Structures Using First- and Second-Order Stochastic Averaging

This paper deals with the dynamic response of nonlinear elastic structure subjected to random hydrodynamic forces and parametric excitation using a first- and second-order stochastic averaging method. The governing equation of motion is derived by using Hamilton's principle, taking into account inertia and curvature nonlinearities and work done due to hydrodynamic forces. Within the framework of first-order stochastic averaging, the system response statistics and stability boundaries are obtained. Unfortunately, the effects of nonlinear inertia and curvature are not reflected in the final results, since the contribution of these nonlinearities is lost during the averaging process. In the absence of hydrodynamic forces, the method fails to give bounded response statistics, and the analysis yields stability conditions. It is the second-order stochastic averaging which can capture the influence of stiffness and inertia nonlinearities that were lost in the first-order averaging process. The results of the second-order averaging are compared with those predicted by Gaussian and non-Gaussian closures and by Monte Carlo simulation. In the absence of parametric excitation, the non-Gaussian closure solutions are in good agreement with Monte Carlo simulation. On the other hand, in the absence of hydrodynamic forces, second-order averaging gives more reliable results in the neighborhood of stochastic bifurcation. However, under pure parametric random excitation, the stochastic averaging and Monte Carlo simulation predict the on-off intermittency phenomenon near bifurcation point, in addition to stochastic bifurcation in probability.     

Stochastic averaging of energy harvesting systems

A stochastic averaging method is proposed for nonlinear energy harvesters subjected to external white Gaussian noise and parametric excitations. The Fokker–Planck–Kolmogorov equation of the coupled electromechanical system of energy harvesting is a three variables nonlinear parabolic partial differential equation whose exact stationary solutions are generally hard to find. In order to overcome difficulties in solving higher dimensional nonlinear partial differential equations, a transformation scheme is applied to decouple the electromechanical equations. The averaged Itô equations are derived via the standard stochastic averaging method, then the FPK equations of the decoupled system are obtained. The exact stationary solution of the averaged FPK equation is used to determine the probability densities of the displacement, the velocity, the amplitude, the joint probability densities of the displacement and velocity, and the power of the stationary response. The effects of the system parameters on the output power are examined. The approximate analytical outcomes are qualitatively and quantitatively supported by the Monte Carlo simulations.     

功能梯度简支矩形板的非线性动力响应

研究了功能梯度简支矩形板在横向简谐激励作用下的非线性动力响应问题。采用幂律分布规律描述功能梯度材料的等效材料参数,基于Galerkin法建立了系统广义坐标的常微分控制方程。利用平均法得到了系统的幅频响应特性,分析了功能梯度矩形薄板的非线性主共振特性。数值算例验证了平均化方法的正确性,揭示了功能梯度平板主共振响应中的多值性和跳跃现象;同时分析发现初始条件会改变功能梯度平板主共振的响应幅值。最后讨论了功能梯度材料的梯度指数对系统幅值响应的影响。     

A modified averaging scheme with application to the secondary Hopf bifurcation of a delayed van der Pol oscillator

In this paper, a modified averaging scheme is presented for a class of time-delayed vibration systems with slow variables. The new scheme is a combination of the averaging techniques proposed by Hale and by Lehman and Weibel, respectively. The averaged equation obtained from the modified scheme is simple enough but it retains the required information for the local nonlinear dynamics around an equilibrium. As an application of the present method, the delay value for which a secondary Hopf bifurcation occurs is successfully located for a delayed van der Pol oscillator.     

A micromechanical analysis of the nonlinear elastic and viscoelastic constitutive relation of a polymer filled with rigid particles

Micromechanical theory is applied to study the nonlinear elastic and viscoelastic constitutive relations of polymeric matrix filled with high rigidity solid particles. It is shown that Eshelby's method can be extended to the case of nonlinear matrix and Eshelby's tensor still exists provided that Poisson's ratio of the nonlinear matrix assumes constant value in deforming process and the rigidity of elastic filling particles is much higher than that of the matrix. A new method for averaging process is proposed to overcome the difficulty that occured in applying the ordinary equivalent inclusion method or the self-consistant method to nonlinear matrices. A rather simple constitutive equation is obtained finally and the strengthening effect of solid particles to composites is investigated. The work supported by the LNM, Institute of Mechanics, Chinese Academy of Sciences and by the National Natural Science Foundation of China     

具有多个极限环非线性动力系统的解析近似

应用一种新的解析方法------同伦分析法，研究了一种具有多个 极限环的Rayleigh振子问题. 与所有其他传统方法不同，该方法不依赖于小参数， 且提供了一个简便的途径以确保级数解的收敛, 因此，特别适用于强非线性问题. 将同伦分析法与平均法以及四阶的龙格库塔方法(数值解)做了比较. 结果 表明，平均法在强非线性情况失效, 四阶的龙格库塔法不能找到非稳定的极限环，而同伦分析法不仅适用于强非线性情 况，而且给出了非稳定的极限环.     

Influences of random disturbances on simultaneous resonances of nonlinear coupling systems

The influences of random disturbances on simultaneous resonances of nonlinear coupling systems are dealt with in this paper. First, the approximate probability distribution of behavior of nonlinear system is presented through using the stochastic averaging method. Secondly, using the Monte Carlo numerical simulation method, we study the above mentioned system. Both conclusions are nearly same. It is confirmed that the stochastic averaging method is one of efficient methods for dealing with nonlinear random vibration problems.     

有限宽轴承-转子系统非线性行为研究

利用平均本征值法求得有限宽轴承的油膜力．以刚性Jeffcott转子为研究对象，借助数值积分并结合Poincare映射研究了轴承-转子系统的非线性动力学行为．所得结果为平均本征值法用于同类型的工程实际问题打下了一定的基础．     

随机平均规范形方法

计算随机规范形系数是应用随机规范形方法的关键．提出一种应用随机平均计算随机规范形系数的方法．为了说明方法的有效性,对白噪声激励的Duffing系统,经过变换,对于相应的平均方程,比较了精确解、规范形方法解和能量包络方法解的稳态概率密度,结果表明,当非线性项系数较小时,三者完全一致．当非线性项系数较大时,规范形方法所得解与精确解相差不大．     

Complexification-averaging method via the averaged Lagrangian

In this paper, the complexification-averaging (CX-A) method for multi-DOF nonlinear vibratory systems is rederived in a new way based upon the averaged Lagrangian. The complex variables are introduced to represent the original displacements and velocities, and then the fast–slow decomposition of the complex variables is made. The time averaging of the Lagrangian over the fast variables is performed. Two different expressions for the kinetic energy are presented, and this results in two schemes for deriving the governing equations of the slow variables. For the scheme I, through the order analysis of the derivatives of the slow variables, it is shown that the second-order terms appeared in the averaged Lagrangian can be omitted, and thus a reduced averaged Lagrangian is obtained. Via the reduced averaged Lagrangian, the corresponding Lagrangian equations are derived. For the scheme II, through time averaging, the averaged Lagrangian is obtained, and then the corresponding equations for the slow variables can be obtained. Finally, two nonlinear vibratory systems with two-DOF and four-DOF, respectively, are given as examples to illustrate the new procedure for the CX-A method. The loci of nonlinear normal modes on the potential surface are studied in the first example, and the frequency-energy plot is investigated in the second example.