含立方刚度和滞回的隔振系统动力学特性分析 1)

包含立方刚度和Bouc-Wen 型滞回的隔振系统具有复杂的非线性动力学特性。系统无阻尼响应模型可基于无滞回恢复力建立,利用谐波平衡法和泰勒展开求得近似解析解。系统有阻尼响应模型可利用解析/数值联合方法求解,该方法基于谐波平衡法和Levenberg-Marquardt 迭代算法,对于滞回产生的多值非光滑函数项,先计算时域响应再通过快速傅里叶变换求解谐波项系数。上述方法在含水平绞制梁的非线性隔振系统分析中得到有效应用。分析表明,在Bouc-Wen 型滞回和立方刚度的综合影响下,隔振系统呈现渐软–渐硬特性,滞回阻尼和线性阻尼都可以有效抑制共振,但前者高频隔振效果优于后者。     

DYNAMICAL ANALYSIS ON A KIND OF SEMI-ACTIVE VIBRATION ISOLATION SYSTEMS WITH DAMPING CONTROL 1)

研究了一类基于相对速度反馈的含立方刚度的单自由度非线性半主动隔振系统.通过平均法得到了系统分别在基于加速度-相对速度反馈的加速度驱动阻尼控制策略、速度-相对速度反馈的天棚阻尼控制策略和位移-相对速度反馈的地棚阻尼控制策略下主共振响应的近似解析解,并利用数值解验证了近似解析解的准确性.通过 Lyapunov 理论对不同控制策略下系统的稳定性进行了分析,讨论了系统参数对控制效果的影响.分析结果表明,对 3 种基于相对速度反馈的控制策略进行解析研究时,切换条件中的控制参数具有相同的表达式;在抑制共振响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略在低频时的减振效果最好,而基于加速度-相对速度反馈的加速度驱动阻尼控制策略在高频时的减振效果最优;在抑制瞬态响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略的减振效果最好.此类解析研究方法可应用到其他半主动开关控制策略中,为半主动隔振系统的控制策略研究提供了有效的方法和手段.     

一类阻尼控制半主动隔振系统的解析研究

研究了一类基于相对速度反馈的含立方刚度的单自由度非线性半主动隔振系统.通过平均法得到了系统分别在基于加速度-相对速度反馈的加速度驱动阻尼控制策略、速度-相对速度反馈的天棚阻尼控制策略和位移-相对速度反馈的地棚阻尼控制策略下主共振响应的近似解析解,并利用数值解验证了近似解析解的准确性.通过 Lyapunov 理论对不同控制策略下系统的稳定性进行了分析,讨论了系统参数对控制效果的影响.分析结果表明,对 3 种基于相对速度反馈的控制策略进行解析研究时,切换条件中的控制参数具有相同的表达式;在抑制共振响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略在低频时的减振效果最好,而基于加速度-相对速度反馈的加速度驱动阻尼控制策略在高频时的减振效果最优;在抑制瞬态响应振幅方面,基于速度-相对速度反馈的天棚阻尼控制策略的减振效果最好.此类解析研究方法可应用到其他半主动开关控制策略中,为半主动隔振系统的控制策略研究提供了有效的方法和手段.      

Nonlinear integral resonant controller for vibration reduction in nonlinear systems

A new nonlinear integral resonant controller (NIRC) is introduced in this paper to suppress vibration in nonlinear oscillatory smart structures. The NIRC consists of a first-order resonant integrator that provides additional damping in a closed-loop system response to reduce high-amplitude nonlinear vibration around the fundamental reso-nance frequency. The method of multiple scales is used to obtain an approximate solution for the closed-loop system. Then closed-loop system stability is investigated using the resulting modulation equation. Finally, the effects of different control system parameters are illustrated and an approximate solution response is verified via numerical simulation results. The advantages and disadvantages of the proposed controller are presented and extensively discussed in the results. The controlled system via the NIRC shows no high-amplitude peaks in the neighboring frequencies of the resonant mode, unlike conventional second-order compensation methods. This makes the NIRC controlled system robust to excitation frequency variations.     

多尺度法的推广及在非线性黏弹性系统的应用

黏弹性材料在航空、机械、土木等领域具有广阔的应用前景,而具有1.5自由度的非线性Zener模型能更好地描述其特性.因此,研究多尺度法的推广和应用具有重要意义.在传统多尺度法的基础上,推广并利用多尺度法对非线性奇数阶微分方程进行研究,解决非线性奇数阶系统的动力学求解问题.以非线性Zener模型为例,首先通过推广的多尺度法...     

动力吸振器复合非线性能量阱对线性镗杆系统的振动控制

研究了线性动力吸振器复合非线性能量阱对线性镗杆在外部简谐激励下的振动控制. 忽略镗杆系统中的非线性因素, 建立了附加线性动力吸振器和非线性能量阱的镗杆系统的三自由度运动方程, 研究了附加复合式动力吸振器的镗杆系统的受迫振动. 通过平均法得到了附加复合式动力吸振器的镗杆系统的近似解析解, 并利用数值解验证了近似解析解的准确性, 两者具有很好的一致性. 利用近似解析解详细分析了线性动力吸振器和非线性能量阱的参数对镗杆振动抑制性能的影响. 对给定质量的复合式动力吸振器进行了参数优化, 其中线性动力吸振器参数采用H_∞优化方法的近似解析解进行了优化, 非线性能量阱的阻尼利用系统的近似解析解进行了优化. 分析结果表明, 线性动力吸振器与非线性能量阱组合可以有效抑制线性镗杆系统的振动, 而且采用参数优化后的复合式动力吸振器可以获得更好的减振效果. 通过附加非线性能量阱, 不但可以提高线性动力吸振器的振动抑制效果, 而且还可以提高振动控制系统的鲁棒性.      

On a Rayleigh Wave Equation with Boundary Damping

In this paper an initial-boundary value problem for a weakly nonlinear string (or wave) equation with non-classical boundary conditions is considered. One end of the string is assumed to be fixed and the other end of the string is attached to a dashpot system, where the damping generated by thedashpot is assumed to be small. This problem can be regarded as a simple model describing oscillations of flexible structures such as overhead transmission lines in a windfield. An asymptotic theory for a class ofinitial-boundary value problems for nonlinear wave equations is presented. Itwill be shown that the problems considered are well-posed for all time t. A multiple time-scales perturbation method incombination with the method of characteristics will be used to construct asymptotic approximations of the solution. It will also be shown that all solutions tend to zero for a sufficiently large value of the damping parameter. For smaller values of the damping parameter it will be shown how the string-system eventually will oscillate. Some numerical results are alsopresented in this paper.     

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

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

Research on linear/nonlinear viscous damping and hysteretic damping in nonlinear vibration isolation systems

A nonlinear vibration isolation system is promising to provide a high-efficient broadband isolation performance. In this paper, a generalized vibration isolation system is established with nonlinear stiffness, nonlinear viscous damping, and Bouc-Wen(BW)hysteretic damping. An approximate analytical analysis is performed based on a harmonic balance method(HBM) and an alternating frequency/time(AFT) domain technique.To evaluate the damping effect, a generalized equivalent damping ratio is defined w...     

多体系统碰撞动力学中接触力模型的研究进展

接触碰撞行为作为大自然与多体系统中的常见现象,其接触力模型对于多体系统的碰撞行为机理研究与性能预测至关重要.静态弹塑性接触模型与考虑能量耗散的连续接触力模型是研究接触碰撞行为的两类不同方法,在多体系统碰撞动力学中存在诸多共性与差异.本文分别从上述两类接触模型的发展历程入手,详细介绍了两类模型的区别与联系.首先,根据阻尼项分母中是否含有初始碰撞速度将连续接触力模型分为黏性接触力模型与迟滞接触力模型,讨论了能量指数与Hertz接触刚度之间的关系,阐述了现有连续接触力模型在计算弹塑性材料接触碰撞行为时存在的问题.其次,着重介绍了分段连续的准静态弹塑性接触力模型(可连续从完全弹性转换到完全塑性接触阶段),分析了利用此类弹塑性接触力模型计算碰撞行为的技术特点.同时,以恢复系数为桥梁和借助线性化的弹塑性接触刚度,避免了Hertz刚度对弹塑性接触刚度的计算误差,根据碰撞前后多体系统的能量与动能守恒推导了弹塑性接触模型等效的迟滞阻尼因子.探索了连续接触力模型与准静态弹塑性接触力模型之间的内在联系,数值计算结果定量说明了人为阻尼项代表的能量耗散与弹塑性接触力模型中加卸载路径代表的能量耗散具有等效性.另外...     

Rheological–dynamical analogy: Prediction of damping parameters of hysteresis damper

This research aims to predict the damping parameters of hysteresis damper based on an analytical rheological–dynamical (RDA) visco-elasto-plastic solution of one-dimensional longitudinal continuous vibrations of a bar. A visco-elasto-plastic bar or damper is an energy dissipation device. An attempt is made to estimate quantitatively the influence of material physical parameters of materials on the damping ratio in both the linear visco-elastic analysis and the nonlinear visco-elasto-plastic analysis of damper subjected to external vibration forces. Two types of damping are considered: viscous damping in the case of linear analysis, defined as stiffness and/or mass proportional and, in the case of nonlinear analysis, hysteresis damping caused by inelastic deformations of damper. Owing to the visco-elastic nature of the materials of the damper and the frequency dependence of the viscous damping ratio ξ, it is useful to consider separately the situations arising when ξ is positive (the system is stable) and when it is negative. A negative damping ratio means that the complementary solution of the response would not die away (the system is unstable because of factor e^{ξ · ω · t}). In the case of nonlinear analysis, the force–displacement relation is nonlinear, so it is very difficult to predict the actual damping and stiffness coefficients, even if the force–displacement characteristic is simply perfect elasto-plastic. Using the RDA method, which takes into account the rate of release of visco-elasto-plastic energy of the dissipation devices; nonlinear behaviors are linearized, enabling to obtain the equivalent damping and stiffness coefficients and the effective period for the damper.     

Dynamic analysis of nonlinear systems by modal synthesis techniques

Different kinds of modal synthesis method have been used widely in dynamic analysis of linear structure systems, but, in general, they are not suitable for nonlinear systems.In this paper, a kind of modal synthesis techniques is extended to dynamic analysis of nonlinear systems. The procedure is based upon the method suggested in [20],[21], which is applicable to vibration analysis for complex structure systems with coupling attachments but with simplified forms of linear springs and dampers. In fact, these attachments have nonlinear characteristics as those generally known to the cases of nonlinear elasticity and nonlinear damping, e.g., piecewise-linear springs, softening or hardening springs. Coulomb damping,elas-ioplastic hysteresis damping, etc. So long as the components of structure are still linear systems, we can get a set of independent free-interface normal mode information hut only keep the lower-order for each component. This can be done by computations or experiments or both. The global equations of linear vibration are set up by assembling of the component equations of motion with nonlinear coupling forces of attachments. Then the problem is reduced to less degrees of freedom for solving nonlinear equations. Thus considerable saving in computer storage and execution time can be expected. In the case of a very high-order system, if sufficient degrees of freedom are reduced, then it may be possible for the problem to be solved by the aid of a computer of ordinary grade.As the general nonlinear vibration of multiple degrees of freedom systems is quite involved, in general, the exact solution of a nonlinear system equations is not easy to find, so the numerical method can be adopted for solving the reduced nonlinear equations to obtain the transient response of system for arbitrary excitations.     

Damping Identification Using Perturbation Techniques and Higher-Order Spectra

Perturbation techniques and spectral moments are combined tocharacterize and quantify the damping and nonlinear parameters of thefirst mode of a three-beam two-mass frame. The frame is excitedharmonically near twice its lowest natural frequency. The response ismodeled with a second-order nonlinear equation with quadratic and cubicterms and linear and quadratic damping terms. The method of multiplescales is used to obtain a second-order approximate solution for thismodel. Bispectral analysis is used to quantify the level of couplingbetween modes and measure their phase difference. The amplitudes andphase difference between the excitation and response mode with differentfrequencies are substituted into the approximate solution to determinethe damping and nonlinear parameters.     

Application of Jacobian Elliptic Functions to the Analysis of the Steady-State Solution of the Damped Duffing Equation with Driving Force of Elliptic Type

In this work, Jacobian elliptic functions are applied to obtain the approximate steady-state solution of the damped Duffing equation with driving force of elliptic type. We assume that the damping coefficient and the nonlinear term need not be small. Numerical comparison between the approximate solution and the numerical integration indicates that our proposed approximate steady-state solution captures well the amplitude-time vibrational response. We have also found that for a certain range of values of the driving frequency ω_f and for a stiff spring, there are one, three, or five possible motions for the elliptic excitation amplitude-frequency response curve.     

分数阶拟周期Mathieu方程的动力学分析

分数阶微积分有着诸多优异的特点, 目前在动力学领域主要用来提高非线性系统振动特性研究的准确性. 本文在拟周期Mathieu方程的基础上, 引入分数阶微积分理论, 研究了分数阶微分项参数对方程稳定性的影响. 首先, 采用摄动法得到方程稳定区和非稳定区分界线(即过渡曲线)近似表达式, 利用数值方法验证了解析结果的准确性, 图像显示两者吻合较好. 随后, 通过归纳总结不同情况下的过渡曲线近似表达式, 发现在系统中分数阶微分项以等效线性刚度和等效线性阻尼的方式存在. 根据这一特点, 得到了系统等效线性阻尼和等效线性刚度的一般形式, 并且定义了非稳定区域厚度. 最后, 通过数值仿真直观地分析了分数阶微分项参数对方程稳定区域大小和过渡曲线位置的影响. 结果发现, 分数阶微分项不仅具有阻尼特性还具有刚度特性, 并且以等效线性刚度和等效线性阻尼的方式影响着方程稳定区域大小和过渡曲线位置. 合理选择分数阶微分项参数可以使其呈现不同程度的刚度特性或阻尼特性, 方程稳定区域的大小和过渡曲线的位置也因此产生了不同程度的变化.      

横向荷载作用下弹性地基上梁的主共振分析

基于Winkler地基模型及Euler-Bernoulli梁理论,建立了弹性地基上有限长梁的非线性运动方程.运用Galerkin方法对运动方程进行一阶模态截断,并利用多尺度法求得该系统主共振的一阶近似解.分析了长细比、地基刚度、外激励幅值和阻尼系数等参数对系统主共振幅频响应的影响,然后通过与非共振硬激励情况对比分析主共振对其动力响应的影响.结果表明:主共振幅频响应存在跳跃和滞后现象;阻尼对主共振响应有抑制作用;主共振显著增大系统稳态动力响应位移.     

基础谐波激励单侧约束简支梁系统的动力学特征分析 1)

为研究间隙碰撞对系统动力学响应的影响,以理想简支梁的振型函数为Rayleigh--Ritz基函数建立了单侧约束简支梁系统的非线性离散动力学方程组,应用数值方法研究了系统在基础谐波激励下的动力学响应特征及其对共振频率线性等效方法适用性的影响。研究表明：非线性动力学系统间隙产生的局部碰撞,使得系统振动能量在系统各阶模态之间转移,使得线性等效方法失效;即使进行非线性分析,也需要考虑系统固有频率远大于激励频段上限的模态。     

Nonlinear dispersion of long internal waves

The propagation of long weakly nonlinear waves in an atmospheric waveguide is considered. A model system of Kadomtsev-Petviashvili equations [1], which describes the propagation of such waves, is derived. In the case of one excited wave mode the system of model equations goes over into the Kadomtsev-Petviashvili equation, in which, however, the variables x and t are interchanged. The reasons for this are clarified. In the two-dimensional case an approximate solution of the model equations is constructed, and steady nonlinear waves and their interaction in a collision are considered. The results of a numerical verification of the stability of the approximate steady solutions and of the solution to the problem of decay of the wave into quasisolitons are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 151–157, May–June, 1988.     

Approximating piecewise nonlinearities in dynamic systems with sigmoid functions: advantages and limitations

In the industry field, the increasingly stringent requirements of lightweight structures are exposing the ultimately nonlinear nature of mechanical systems. This is extremely true for systems with moving parts and loose fixtures which show piecewise stiffness behaviours. Nevertheless, the numerical solution of systems with ideal piecewise mathematical characteristics is associated with time-consuming procedures and a high computational burden. Smoothing functions can conveniently simplify the mathematical form of such systems, but little research has been carried out to evaluate their effect on the mechanical response of multi-degree-of-freedom systems. To investigate this problem, a slightly damped mechanical two-degree-of-freedom system with soft piecewise constraints is studied via numerical continuation and numerical integration procedures. Sigmoid functions are adopted to approximate the constraints, and the effect of such approximation is explored by comparing the results of the approximate system with the ones of the ideal piecewise counter-part. The numerical results show that the sigmoid functions can correctly catch the very complex dynamics of the proposed system when both the above-mentioned techniques are adopted. Moreover, a reduction in the computational burden, as well as an increase in numerical robustness, is observed in the approximate case.

     

聚苯胺体系电流变流体的阻尼机理及阻尼模型研究

对聚苯胺体系电流变流体在正弦信号作用下的阻尼力响应及其频谱特性进行了理论和实验研究.提出电流变流体的非线性阻尼力的频谱是由基频和频率为基频奇数倍的高次谐波分量组成.建立了由粘性阻尼和迟滞阻尼组成的电流变流体阻尼模型,通过对模型的傅立叶变换及数值仿真,验证了理论分析及阻尼模型的正确性.