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

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

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

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

黏弹性轴向运动梁非线性受迫振动稳态幅频响应

本文研究了黏弹性轴向运动梁横向受迫振动稳态幅频响应问题.在控制方程的推导中,对黏弹性本构关系采用物质导数.把多尺度法直接应用于梁横向振动的非线性控制方程,利用可解性条件消除长期项,得到系统稳态的幅频响应曲线.运用Lyapunov一次近似理论分析幅频响应曲线的稳定性.通过算例研究了黏性系数,外部激励幅值以及非线性项系数对稳态幅频响应曲线及其稳定性的影响.运用数值方法对两端固定边界下黏弹性轴向运动梁的控制方程直接数值解,分析梁横向非线性振动的稳态幅频响应,通过数值算例验证直接多尺度法的结论.     

多关节机械臂系统非线性动力学分析

针对工程机械臂的非线性振动，研究了系统因周期运动失稳而产生Hopf分岔问题．工程机械臂在建模时往往忽略非线性因素，然而在特殊场合下非线性的影响不可忽略，本文利用拉格朗日法建立了新的非线性两关节机械臂系统的两自由度动力学方程，采用多尺度法对机械臂系统模型进行了求解，对系统发生主共振时的动力学行为进行了分析讨论．结果表明所求得的近似解析解与数值解具有较好的吻合性，同时由于系统方程中非线性项的影响，随着激励幅值与激励频率的变化，系统的振幅会发生跳跃，出现动态分岔现象．本研究为实际工程应用提供了一些理论依据和指导．     

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题. 对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应 的精确解. 遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解. 事实上,其数值解法严重受限于方程维度,而解析求解 则仅适用于少数特定的系统,且多是稳态解. 因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径. 本文针 对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降 维. 针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程. 建议了构造等价漂移系数 的条件均值函数方法. 进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答. 结合单自由度Rayleigh 振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精 度和效率,验证了其有效性.      

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题.对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应的精确解.遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解.事实上,其数值解法严重受限于方程维度,而解析求解则仅适用于少数特定的系统,且多是稳态解.因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径.本文针对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降维.针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程.建议了构造等价漂移系数的条件均值函数方法.进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答.结合单自由度Rayleigh振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精度和效率,验证了其有效性.     

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

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

轮对非线性动力学系统蛇行运动的解析解

在对铁路车辆系统的极限环幅值和非线性临界速度进行分析时通常采用数值方法, 不便于研究其随系统参数的变化规律. 轮对系统保留了影响车辆系统动力学性能的几个关键要素: 如轮轨几何非线性约束、轮轨接触蠕滑关系和悬挂系统等, 可以反映铁路车辆系统蛇行运动的本质特性. 轮对系统自由度少、参数少, 可以采用解析方法进行分析. 本文选取合适的特征量把轮对非线性动力学方程无量纲化, 得到了带有小参数的两自由度微分方程; 采用多尺度方法对该方程进行了解析求解; 给出了轮对系统极限环幅值的解析表达式并对其稳定性进行了判定; 给出了轮对系统的分岔速度解析表达式, 并进而获得系统的非线性临界速度的解析表达式. 在对得到的解析解用数值结果进行验证后, 用得到的解析解进行了系统参数影响分析. 传统的分岔图计算方法(如降速法、路径跟踪法等)需对微分方程进行大量数值积分计算方可求解系统的非线性临界速度值, 而通过本文获得的解析表达式可直接给出系统的非线性临界速度值和极限环幅值, 便于研究轮对系统动力学特性随参数的变化规律,进行快速方案比对和筛选, 为转向架结构优化设计提供参考.      

Duffing系统的主--亚谐联合共振

以Duffing系统为研究对象,研究在多频激励下同时发生主共振和1/3次亚谐共振的动力学行为与稳定性.首先,通过多尺度法得到系统的近似解析解,利用数值方法检验近似程度,结果吻合良好,证明了求解过程和解析解的正确性.然后,从解析解中导出稳态响应的幅频方程和相频方程,从幅频曲线以及相频曲线中发现系统最多存在7个不同的周期解,这种多解现象可用于对系统状态进行切换.基于Lyapunov稳定性理论,得到联合共振定常解的稳定条件,利用该条件分析了系统的稳定性,并与Duffing系统的主共振和1/3次亚谐共振单独存在时比较.最后,通过数值方法分析了非线性项和外激励对系统动力学行为与稳定性的影响,发现了联合共振特有的现象:刚度软化时,非线性项不仅影响系统的响应幅值,同时还影响系统的多值性和稳定性;刚度硬化时,非线性项对系统的影响与单一频率下主共振和1/3次亚谐共振类似,仅影响系统的响应幅值.这些结果对Duffing系统动力学特性的研究具有重要意义.     

粘弹性矩形板的混沌和超混沌行为

从薄板Karman理论的基本假设出发；利用线性粘弹性理论中的Boltzman叠加原理，建立了粘弹性薄板非线性动力学分析的初边值问题，其运动方程是一组非线性积分──微分方程．在空间域上利用Galerkin平均化法之后，得到了变型的非线性积分──微分型的Duffing方程．综合利用动力系统中的多种方法，揭示了粘弹性矩形板在横向周期激励下的丰富的动力学行为，如不动点、极限环、混沌、奇怪吸引子、超混沌等，其中，混沌和超混沌是交替出现的．     

Thermal-induced snap-through buckling of simply-supported functionally graded beams

The instability of functionally graded material (FGM) structures is one of the major threats to their service safety in engineering applications. This paper aims to clarify a long-standing controversy on the thermal instability type of simply-supported FGM beams. First, based on the Euler-Bernoulli beam theory and von Kármán geometric nonlinearity, a nonlinear governing equation of simply-supported FGM beams under uniform thermal loads by Zhang's two-variable method is formulated. Second, an approximate analytic solution to the nonlinear integro-differential boundary value problem under a thermal-induced inhomogeneous force boundary condition is obtained by using a semiinverse method when the coordinate axis is relocated to the bending axis (physical neutral plane), and then the analytical predictions are verified by the differential quadrature method (DQM). Finally, based on the free energy theorem, it is revealed that the symmetry breaking caused by the material inhomogeneity can make the simply-supported FGM beam under uniform thermal loads occur snap-through postbuckling only in odd modes; furthermore, the nonlinear critical load of thermal buckling varies non-monotonically with the functional gradient index due to the stretching-bending coupling effect. These results are expected to provide new ideas and references for the design and regulation of FGM structures.     

Swell-induced surface instability of confined hydrogel layers on substrates

In response to external stimuli, polymeric hydrogels can change volume and shape dramatically. Experimental studies have observed a variety of instability patterns of hydrogels, due to swelling or shrinking, many of which have not been well understood. The present paper considers swell-induced surface instability of a hydrogel layer on a rigid substrate. Based on a recently developed theoretical framework for neutral polymeric gels, a linear perturbation analysis is performed to predict the critical condition for the onset of the surface instability. Using a nonlinear finite element method, numerical simulations are presented to show the swelling process, with the evolution of initial surface perturbations followed by the formation of crease-like surface patterns. In contrast to previously suggested critical conditions for surface creasing, the present study suggests a material-specific condition that predicts a range of critical swelling ratios from about 2.5 to 3.4 and quantitatively relates the critical condition to material properties of the hydrogel system. A stability diagram is constructed with two distinct regions for stable and unstable hydrogels depending on two dimensionless material parameters.     

横向载荷作用下轴向运动梁的屈曲失稳以及非线性振动特性

梁的轴向运动会诱发其产生横向振动并可能导致屈曲失稳，对结构的安全性和可靠性产生重大的影响。本文重点研究了横向载荷作用下轴向运动梁的屈曲失稳及横向非线性振动特性。基于Hamilton变分原理，建立了横向载荷作用下轴向运动梁的动力学方程，获得了梁的后屈曲构型。使用截断Galerkin法，将控制方程改写成Duffing方程的形式。用同伦分析方法确定载荷作用下轴向运动梁的非线性受迫振动的封闭形式的表达式。结果表明，后屈曲构型对轴向速度和初始轴向应力有明显的依赖性。通过同伦分析法得出非线性基频的显式表达式，获得了初始轴向力会影响非线性频率随初始振幅和轴向速度的线性关系。另外，轴向外激励的方向也会改变系统固有频率。     

DYNAMIC PULL-IN INSTABILITY OF DOUBLE-SIDED ACTUATED NANO-TORSIONAL SWITCHES

A nonlinear frequency-amplitude relation is developed to investigate the vibrational amplitude effect on the dynamic pull-in instability of double-sided-actuated nano-torsional switches. The governing equation of a nano-electro-mechanical system pre-deformed by an electric field contains the quintic nonlinear term. The influences of basic parameters on the pull-in instability and natural frequency are investigated using a powerful analytical approach called the homotopy perturbation method. It is demonstrated that two terms in series expansion are sufficient to produce an acceptable solution. The numerical results obtained have verified the soundness of the asymptotic procedure. The phase portraits of the double-sided nano-torsionalactuator exhibit periodic, homoclinic and heteroclinic orbits.     

Nonlinear Dynamics and Stability of a Machine Tool Traveling Joint

In this work, we investigate the nonlinear dynamics and stability of a machine tool traveling joint. The dynamical system considered includes contacting elements of a lathe joint and the cutting process where the onset of instability is governed by mode coupling. The equilibrium equations of the dynamical system yield a unique fixed point that can change its stability via a Hopf bifurcation. The unstable domain is primarily governed by the cutting tool location, the contact stiffness of the joint and the depth of material to be removed. Self excited vibrations due to a mode coupling instability evolve around the unstable fixed point and one or more limit cycles may coexist in the statically unstable domain. Stability and accuracy of the approximate analytical solutions are analyzed by applying Floquet analysis. Perturbation of the dynamical system with weak periodic excitation results with periodic and aperiodic solutions.     

Horizontal thermal convection in a porous medium

Linearised instability and nonlinear stability bounds for thermal convection in a fluid-filled porous finite box are derived. A nonuniform temperature field in the steady state is generated by maintaining the vertical walls at different temperatures. The linearised instability threshold is shown to be well above the global stability boundary, which is strongly suggested by the lack of symmetry in the perturbed system. The numerical results are evaluated utilising a newly developed Legendre-polynomial-based spectral method.     

Chaotic belt phenomena in nonlinear elastic beam

IntroductionThechaoticphenomenainsolidmechanicsfieldsbringmoreandmoreinterest.In 1 998,F .C .Moon[1]analyzedthechaoticbehaviorsofbeamsexperimentallyfirst.Thenhestudiedthedynamicsresponseoflinearelasticbeamsubjectedtransverseperiodicload .Thechaoticmotionsoflineardampingbeamshavebeenstudiedbymanyscholarsathomeandabroadinrecentyears[2 ,3].ThedynamicbehaviorsofnonlineardampingbeamssubjectedtotransverseloadP=δP0 (f+cosωt)sin(πx/l)arestudiedinthispaper.Thecriticconditionsthatchaosoccursinthes…     

Weakly nonlinear stability analysis of a falling film with countercurrent gas flow

Weakly nonlinear stability analysis of a falling film with countercurrent gas–liquid flow has been investigated. A normal mode approach and the method of multiple scales are employed to carry out the linear and nonlinear stability solutions for the film flow system. The results show that both supercritical stability and subcritical instability are possible for a film flow system when the gas flows in the countercurrent direction. The stability characteristics of the film flow system are strongly influenced by the effects of interfacial shear stress when the gas flows in the countercurrent direction. The effect of countercurrent gas flow in a falling film is to stabilize the film flow system.     

Material and structural instabilities of single-wall carbon nanotubes

The nonlinear atomistic interactions usually involve softening behavior. Instability resulting directly from this softening are called the material instability, while those unrelated to this softening are called the structural instability. We use the finite-deformation shell theory based on the interatomic potential to show that the tension instability of single-wall carbon nanotubes is the material instability, while the compression and torsion instabilities are structural instability.     

Instability analysis and simulation of water infiltration into an unsaturated elasto-viscoplastic material

A linear instability analysis was performed in order to investigate which variables have a significant effect on the onset of the instability of an unsaturated viscoplastic material subjected to water infiltration. It was found that the onset of the growing instability of the material system mainly depends on the specific moisture capacity, the suction, and the hardening parameter. Then, in order to simulate the water infiltration process of a one-dimensional unsaturated soil column, a multiphase coupled elasto-viscoplastic finite element analysis was performed based on the theory of porous media. The results of the numerical simulations are discussed with respect to the effect of the specific moisture capacity and the initial suction on the development of volumetric strain. We found that rapid transitions from unsaturated to saturated states and higher levels of initial suction lead to the contractive behavior of the material and instability. The instability detected by the numerical results is consistent with the theoretical results obtained through the linear instability analysis.