ANALYSIS OF FRICTION INDUCED VIBRATION UNDER THE ACTIVE DISTURBANCE REJECTION CONTROL FRAMEWORK 1)

自抗扰控制(active disturbance rejection control, ADRC)是一种具有两自由度控制结构的工程化方法, 由于其能够直观有效地处理多种扰动, 近些年来在许多机电系统上得到了成功应用. 当采用ADRC对带有摩擦力的机电系统进行调节时, 可能会产生极限环振动. 目前, 还没有ADRC框架下摩擦力振动精确分析的相关工作. 因此, 本文采用非线性动力学系统的分析工具对这一问题进行研究. 首先, 考虑两种典型摩擦力模型, 静态切换模型和动态LuGre 模型, 对一类二阶运动系统设计不同阶次的ADRC, 得到控制器的等效形式, 并揭示出与比例积分微分(proportional-integral-derivative, PID)控制之间的联系. 然后, 采用打靶法结合拟弧长延拓方法求解系统中的极限环, 并根据Floquet理论判断极限环的稳定性、可能出现的分岔以及分岔类型. 此外, 通过雅克比矩阵和近似数值方法对系统平衡点集的局部稳定性进行了分析. 最后, 通过数值计算研究了摩擦力模型和参数、ADRC阶次和参数对极限环和平衡点集的影响. 计算结果表明, 决定摩擦力Stribeck效应负斜率的参数$\beta$作用较大. 当$\beta>1$时, 两种摩擦力模型下的闭环系统呈现出相同的特性, 极限环会出现环面折叠分岔(cyclic fold bifurcation, CFB)且平衡点集是局部稳定的. 然而当$\beta<1$时, 两种闭环系统呈现出完全不同的特性. 此外, 不同阶次的ADRC在极限环的存在性和稳定性、平衡点集的稳定性上面的结论是相同的, 而低阶次的ADRC能够更好地解决摩擦力补偿和稳定鲁棒性之间的矛盾问题. 这些结论对实际现象的理解、ADRC阶次的选择以及参数整定提供了一定指导.     

热环境下壁板非线性颤振分析

基于一阶活塞气动力理论,采用Von Karman大变形应变-位移关系建立了无限展长壁板热环境下颤振方程,采用伽辽金方法对方程进行离散处理.取温度为分叉参数,研究壁板颤振时的分叉及混沌等复杂动力学特性.结果表明:温度载荷降低了系统的颤振临界动压,改变了颤振特性.在整个分岔参数范围内,系统呈现出较为复杂的变化,包括衰减振动、极限环振动、拟周期振动和混沌型振动.当考虑材料热效应时,系统的颤振动压将进一步降低,其响应也表现出更为丰富的非线性动态力学行为.     

一类混沌系统中的簇发振荡及其延迟叉形分岔行为

由于多时间尺度问题在实际工程系统中广泛存在,关于其复杂动力学行为及其产生机制的研究已成为当前国内外的热点课题之一.簇发振荡是多时间尺度系统复杂动力学行为的典型代表,而分岔延迟又是簇发振荡中的常见现象.本文为探讨非线性系统中分岔延迟所引发的簇发振荡的分岔机制,在一个三维混沌系统中引入参数激励,当激励频率远小于系统的固有频率时,系统产生了两时间尺度簇发振荡.将整个激励项看做慢变参数,激励系统转化为广义自治系统也即快子系统,分析快子系统平衡点的稳定性以及分岔条件,并运用快慢分析法和转换相图揭示了簇发振荡的动力学机理.文中考察了4组参数条件下系统的动力学行为,研究发现当慢变激励项周期性地通过分岔点时,系统产生了明显的超临界叉形分岔延迟行为,随着参数激励振幅的增大,分岔延迟的时间也逐渐延长,当这种延迟的动态行为终止于不同的参数区域时,导致系统轨线围绕不同稳定吸引子(平衡点,极限环)运动,从而得到了不同的簇发振荡行为.      

凸肩叶片的非线性振动特性与运动分岔

以凸肩叶片作为研究模型, 建立了考虑凸肩摩擦力, 几何大变形与阻尼的非线性振动方程.采用Galerkin法对振动方程离散化, 应用平均法对离散后模态方程组的非线性响应进行解析分析, 得到了非线性幅频特性曲线, 与数值解比较验证了解析解, 并讨论了系统周期解的稳定性. 用非线性振动理论详细研究了平均方程组的运动分岔现象, 揭示了平均方程组周期解的变化过程及其具有的非线性动力学性质. 解析结果表明, 凸肩之间的摩擦对系统第二阶非线性振动特性影响很大. 由于凸肩之间摩擦力方向的不断改变, 系统第二阶非线性幅频特性曲线不连续, 出现两个共振频域. 随着时间的推移, 系统振动的幅值会以$T/ 4$为周期在两个频域的幅频曲线上来回跳动, 这会使叶片的振动响应大幅降低.      

一类混沌系统中的簇发振荡及其延迟叉形分岔行为

由于多时间尺度问题在实际工程系统中广泛存在,关于其复杂动力学行为及其产生机制的研究已成为当前国内外的热点课题之一.簇发振荡是多时间尺度系统复杂动力学行为的典型代表,而分岔延迟又是簇发振荡中的常见现象.本文为探讨非线性系统中分岔延迟所引发的簇发振荡的分岔机制,在一个三维混沌系统中引入参数激励,当激励频率远小于系统的固有频率时,系统产生了两时间尺度簇发振荡.将整个激励项看做慢变参数,激励系统转化为广义自治系统也即快子系统,分析快子系统平衡点的稳定性以及分岔条件,并运用快慢分析法和转换相图揭示了簇发振荡的动力学机理.文中考察了4组参数条件下系统的动力学行为,研究发现当慢变激励项周期性地通过分岔点时,系统产生了明显的超临界叉形分岔延迟行为,随着参数激励振幅的增大,分岔延迟的时间也逐渐延长,当这种延迟的动态行为终止于不同的参数区域时,导致系统轨线围绕不同稳定吸引子(平衡点,极限环)运动,从而得到了不同的簇发振荡行为.     

高超声速飞行器横侧向失稳非线性分岔分析

针对滑翔式高超声速飞行器大攻角横侧向失稳问题,采用延拓算法和分岔理论,求解并分析了以俯仰舵偏为连续参数的稳态平衡分岔图和以副翼舵偏为连续参数的横侧向机动稳态平衡分岔图,对平衡分支的稳定性和突变点进行了分析,并给出了特征根拓扑结构变化.研究表明,高超声速飞行器存在极限分岔点、Hopf分岔点以及叉型分岔点,且从叉型分岔点延伸出多个平衡分支,引起横侧向的自滚转失稳;从Hopf分岔点延伸出极限环分支,该分支对应较为复杂的极限环运动,其中还包含倍周期分岔、花环分岔、极限环极限点分岔等复杂的分岔现象;在横侧向机动飞行情况下,模型存在横向操作偏离失稳问题,且存在多个不稳定的平衡点.研究结果为实现高超声速飞行器的稳定飞行和控制器的设计提供了极其重要的动力学信息.      

电磁作用激发的水电机组转子轴系振动研究

水轮发电机组的振动问题是长期困扰广大科研工作者的一个难题,也一直是学术界研究的热点.本文针对某大型水轮发电机转子轴系,首先建立了轴系扭振的微分方程组,给出了轴系整体扭振零频的计算公式,并计算了零频数值.接着在考虑了转子静偏心、振动偏心和转动偏心的综合影响后,求得了气隙磁场能量和单边电磁拉力的表达式,并应用机电耦联系统的能量法,求得了系统由电磁力激发的组合强迫共振响应的一次近似解.最后并应用奇异性理论对分岔方程进行了分析,求得了分岔集、滞后集、双极限集和转迁集的参数区间,给出了普适开折平面和物理参数平面上的分岔图,求得了机组稳定运行区域的电磁参数.本文所得结论对水电机组的工程设计人员具有重要的指导意义.     

二阶自治Birkhoff系统的平衡点分岔

研究二阶自治Birkhoff系统的奇点、闭轨和极限环,以及与其相关的稳定性问题。给出奇点判据和闭轨判据。应用这些判据讨论了二阶自治Birkhoff系统的平衡点分岔。     

非线性转子——密封系统低频失稳机理研究：不平衡系统的亚谐共振

本文采用Ｍｕｓｚｙｎｓｋａ密封力模型分析单圆盘转子－－密封系统的低频自激振动。文（１）研究了平衡转子的稳定性和分岔，本文研究不平衡转子在临界平衡点附近自激振动（周期扰动Ｈｏｐｆ分岔）的亚谐共振，给出了不同参数条件下的振动性态，为识别转子的亚谐共振故障及预防提供了一些新理论依据。     

冲击渐进振动系统相邻基本振动的转迁规律

冲击振动现象广泛存在于动力机械系统中,使得系统表现出复杂的动力学响应.目前对冲击振动系统的p/1类基本振动的稳定性及分岔研究报道较少,而且已有的对冲击振动系统动力学的研究基本都是基于单参数分岔进行分析的.研究以小型振动冲击式打桩机为工程背景,建立了冲击渐进振动系统的力学模型.分析了激振器和缓冲垫发生碰撞的类型,以及滑块渐进运动的条件.给出了系统可能呈现的四种运动状态的判断条件和运动微分方程.通过二维参数分岔分析得到系统在(ω,l)参数平面内存在的各类周期振动的参数域和分布规律.详细分析了相邻p/1类基本振动的转迁规律.在5/1基本振动的参数域的右边区域,相邻p/1基本振动的参数域临界线上存在一个奇异点X_p,相邻p/1类基本振动的分岔特点以奇异点X_p为临界点.在l小于l_X_p的区域内,相邻p/1基本振动经实擦边分岔和鞍结分岔相互转迁,实擦边分岔线和鞍结分岔线之间存在迟滞域,迟滞域内,系统存在两个周期吸引子共存的现象.在l大于l_X_p的区域内,相邻p/l类基本振动的参数域之间存在一个中间过渡区域.中间过渡区域内,系统呈现(2p+2)/2和(2p+1)/2周期振动等.在5/1基本振动的参数域的左边区域,p/1基本振动经多重滑移分岔产生(P+1)/1基本振动.     

Non-linear control of friction-induced self-excited vibration

The present paper introduces a new method of controlling friction-driven self-excited vibration. The control law is primarily derived using the Lyapunov's second method. A single degree-of-freedom oscillator on a moving belt represents the primary model of the system. The control action is achieved by modulating the normal load at the frictional interface based on the state of the oscillatory system. The basic mechanism of the control action utilises subcritical Hopf bifurcation of the equilibrium followed by cyclic-fold bifurcation (of limit cycle oscillations) to globally stabilise the equilibrium. The basic mechanism is qualitatively independent of the exact model of friction. Different models of friction, like, algebraic model, LuGre model and switch model with time-dependent static friction are considered to substantiate the above claim. An approximate method for estimating the critical value of the control parameter that ensures global stability of the equilibrium is also proposed.     

Stability and Neimark-Sacker bifurcation analysis for a discrete single genetic negative feedback autoregulatory system with delay

In this paper, we deal with a discrete single genetic negative feedback autoregulatory system with delay by using Euler method. Choosing the delay $\tau $ as the bifurcation parameter and analyzing the associated characteristic equation corresponding to the unique positive fixed point, it is found that the stability of the positive equilibrium and Neimark-Sacker bifurcation may occur when $\tau $ crosses some critical values. Then the explicit formula which determines the stability, direction, and other properties of bifurcating periodic solution is derived by using the center manifold theorem and normal form theory. Finally, in order to illustrate our theoretical analysis, numerical simulations are also included in the end.     

Two-parameter bifurcation analysis of limit cycles of a simplified railway wheelset model

The effect of the nonlinear terms on bifurcation behaviors of limit cycles of a simplified railway wheelset model is investigated. At first, the stable equilibrium state loses its stability via a Hopf bifurcation. The bifurcation curve is divided into a supercritical branch and a subcritical one by a generalized Hopf point, which plays a key role in determining the occurrence of flange contact and derailment of high-speed railway vehicles, and the occurrence of this critical situation is an important decision-making criteria for design parameters. Secondly, bifurcations of limit cycles are discussed by comparing the bifurcation behavior of cycles for two different nonlinear parameters. Unlike local Hopf bifurcation analysis based on a single bifurcation parameter in most papers, global bifurcation analysis of limit cycles based on two bifurcation parameters is investigated, simultaneously. It is shown that changing nonlinear parameter terms can affect bifurcation types of cycles and division of parameter domains. In particular, near the branch points of cycles, two symmetrical limit cycles are created by a pitchfork bifurcation and then two symmetrical cycles both undergo a period-doubling bifurcation to form two stable period-two cycles. Around the resonant points, period orbits can make several turns, whose number of turns corresponds to the ratio of resonance. Thirdly, near the Neimark–Sacker bifurcation of cycles, a stable torus is created by a supercritical Neimark–Sacker bifurcation, which shows that the orbit of the model exhibits modulated oscillations with two frequencies near the limit cycle. These results demonstrate that nonlinear parameter terms can produce very complex global bifurcation phenomena and make obvious effects on possible hunting motions even though a simple railway wheelset model is concerned.

     

BIFURCATIONS OF AIRFOIL IN INCOMPRESSIBLE FLOW

Bifurcations of an airfoil with nonlinear pitching stiffness in incompressible flow are investigated. The pitching spring is regarded as a spring with cubic stiffness. The motion equations of the airfoil are written as the four dimensional one order differential equations. Taking air speed and the linear part of pitching stiffness as the parameters, the analytic solutions of the critical boundaries of pitchfork bifurcations and Hopf bifurcations are obtained in 2 dimensional parameter plane. The stabilities of the equilibrium points and the limit cycles in different regions of 2 dimensional parameter plane are analyzed. By means of harmonic balance method, the approximate critical boundaries of 2-multiple semi-stable limit cycle bifurcations are obtained, and the bifurcation points of supercritical or subcritical Hopf bifurcation are found. Some numerical simulation results are given.     

On the effect of nonsmooth Coulomb friction on Hopf bifurcations in a 1-DoF oscillator with self-excitation due to negative damping

This article deals with the question, to what extent damping due to nonsmooth Coulomb friction may affect the stability and bifurcation behavior of vibrational systems with self-excitation due to negative effective damping which??for the smooth case??is related to a Hopf bifurcation of the steady state. Without damping due to Coulomb friction, the stability of the trivial solution is controlled by the effective viscous damping of the system: as the damping becomes negative, the steady state loses stability at a Hopf point. Adding Coulomb friction changes the trivial solution into a set of equilibria, which??for oscillatory systems??is asymptotically stable for all values of effective viscous damping. The Hopf point vanishes and an unstable limit cycle appears which borders the basin of attraction of the equilibrium set. Moreover, the influence of nonlinear damping terms is discussed. The effect of Coulomb frictional damping may be seen as adding an imperfection to the classical smooth Hopf scenario: as the imperfection vanishes, the behavior of the smooth problem is recovered.     

On the effect of non-smooth Coulomb damping on flutter-type self-excitation in a non-gyroscopic circulatory 2-DoF-system

This article deals with self-excited vibrations, attractivity of stationary solutions, and the corresponding bifurcation behavior of two-dimensional differential inclusions of the type $\mathbf{M}\mathbf{q}'' + \mathbf{D}\mathbf{q}' + (\mathbf{K} + \bar{\mu}\mathbf{N})\mathbf{q} \in-\mathbf{R}\operatorname{Sign}(\mathbf{q}')$ . For the smooth case R=0, the equilibrium may become unstable due to non-conservative positional forces stemming from the circulatory matrix N. This type of instability is usually referred to as flutter instability and the loss of stability is related to a Hopf bifurcation of the steady state, which occurs for a critical parameter $\bar{\mu}= \bar{\mu}_{\mathrm{crit}}$ . For R≠0, the steady state is a set of equilibria, which turns out to be attractive for all values of the bifurcation parameter $\bar{\mu}$ . Depending on $\bar{\mu}$ , the basin of attraction of the equilibrium set can be infinite or finite. The transition from an infinite to a finite basin of attraction occurs at the stability threshold $\bar{\mu}_{\mathrm{crit}}$ of the underlying smooth problem. For the finite basin of attraction, its size is proportional to the Coulomb friction and inverse-proportional to $(\bar{\mu}- \bar{\mu}_{\mathrm{crit}})$ . By adding Coulomb damping the notion of steady state stability for the smooth problem is replaced by the question whether the basin of attraction of the steady state is infinite or finite. Simultaneously, the local Hopf-bifurcation is replaced by a global bifurcation. This implies that in the presence of Coulomb damping the occurrence of self-excited vibrations can only be investigated with regard to the perturbation level.     

Local Bifurcation Analysis of a Second Gradient Model for Deformations of a Rectangular Slab

In this paper we carry out a derivation of the equilibrium equations of nonlinear elasticity with an added second-gradient term proportional to a small parameter . These equations are given by a fourth order semilinear system of pdes. We discuss different types of possible boundary conditions for these equations. We then specialize the equations to a rectangular slab and study the linearized problem about a homogenous deformation. We show that these equations admit solutions representable as Fourier series in one of the independent variables. Furthermore, we obtain the characteristic equation for the eigenvalues (possible bifurcation points) for the linear problem and derive asymptotic representations for this equation for small . We used these expressions to show that in the limit as the characteristic equation for converges uniformly (in certain regions of the parameter space) to the corresponding characteristic equation for . When the base material () is that of a Blatz–Ko type, we get conditions for the existence of eigenvalues of the linear problem with and small. Our numerical results in this case indicate that the number of bifurcation points is finite when and that this number monotonically increases as . For the problem with we get conditions for the existence of local branches of non-trivial solutions.