微尺度悬臂管颤振的有限维研究

基于修正的偶应力理论并考虑Lagrange应变张量所给出的几何非线性,运用Hamilton原理建立了微尺度悬臂管平面振动的积分-微分方程.通过Galerkin方法将原积分-微分方程离散成常微分方程组,研究了临界流速-质量比曲线的不同阶Galerkin近似解与精确解的符合程度以及它们对材料长度尺寸参数的依赖性.对不同的模态截断数,运用基于中心流形-范式理论的投影法计算了临界流速处系统的第一Lyapunov(李雅谱诺夫)系数和临界特征值关于流速的变化率,以此为基础分析了系统的分岔模式,探讨了模态截断数对系统动力学性质的影响.临界流速-质量比曲线的滞后部分及交点处的动力学性质表明,系统存在不同的分岔方向,用6个模态的Galerkin离散化方程作分岔图对此进行了验证,并通过理论分析及数值方法分别计算了颤振的固有频率.     

超声速蒙皮颤振微分积分方程的特征值问题

根据二维线化理论讨论超声速薄钣的动力稳定性,导致一类新颖的数学物理问题:非自共轭Volterra型四阶微分积分方程的复特征值问题.求得这一气动弹性系统的严格解.与其它近似分析对比,本法的临界曲线与实验数据符合良好,在低超声速范围不存在发散问题.此外,在数学物理实质方面,发现:(1)颤振频谱与固有频谱有互为间隔现象;(2)临界Mach数有简并现象.指出本法可以推广应用于三维机翼模型和燃气轮中叶栅的超声速颤振问题.     

Duffing系统在双参数平面上的分岔演化过程

给出了参数空间上最大Lyapunov指数的计算方法,数值计算了Duffing系统在双参数平面上的最大Lyapunov指数.结合单参数最大Lyapunov指数、分岔图、相图以及时间历程图,讨论了Duffing系统在双参数平面上的分岔以及随系统控制参数变化的分岔演化过程.结果发现在双参数平面上系统发生叉式分岔,出现具有缺边现象的两个不同区域,该区域内系统对初值有较强的敏感性,存在两吸引子共存现象;系统运动经过周期跳跃曲线时振动幅值突然减小;系统外激励频率较小时常引起颤振运动.此外,在两个具有缺边现象的区域内,随刚度系数的不断增加,系统出现了倍周期分岔曲线环,而且倍周期分岔曲线环内不断嵌套新的倍周期分岔曲线环,导致系统最终经倍周期分岔序列进入混沌状态,随着控制参数的变化,系统在双参数平面上的动力学特性变得非常复杂.     

超音速飞行器机翼颤振的时滞反馈控制

采用时滞反馈主动控制方法对超音速飞行器机翼颤振进行控制,以提高飞行器机翼系统的颤振临界速度.首先根据二元机翼的力学模型,制定时滞反馈控制策略并建立时滞反馈控制系统的数学模型;分别对无控、零时滞反馈控制和有时滞反馈控制系统进行稳定性分析,获得时滞反馈控制系统的颤振稳定性边界.利用MATLAB/SIMULINK进行时域数值模拟,验证理论稳定性分析结果的正确性.结果表明:通过调节时滞量,可有效提高飞行器机翼的颤振临界速度,且控制策略简单,效果较好.     

基于信息压缩矩阵特征值映射的时变系统UD分解辨识算法

在系统辨识领域遗忘因子UD分解算法(一种通过对系统数据矩阵进行UD分解的在线辨识算法)具有对时变系统阶次和参数同步估计的优异性能,但传统的遗忘策略不能从根本上解决信息压缩矩阵数据过饱和问题,为了拓展现有UD分解算法在时变系统的适用范围,同时针对数据空间分布不均匀性,提出一种基于信息压缩矩阵特征值映射的UD分解辨识算法.从理论上分析辨识算法跟踪能力与参数估计矩阵有界性的对应关系,从而构造出一种基于信息压缩矩阵特征值映射的有界函数,特征值映射函数能够根据系统数据传递过程中信息量的大小动态调整遗忘因子,解决了参数辨识过程中数据过饱和及数据分布不均匀问题.仿真结果表明,相比于常规时变遗忘因子策略,带有特征值映射的UD分解算法能够更加准确跟踪系统参数的变化,且能够保证系统不是2N阶持续激励信号的情况下,也能对时变系统参数进行跟踪.     

多体气动弹性系统超声速颤振分析

对一类多体气动弹性系统超声速颤振问题进行研究.分别采用多体动力学理论、活塞理论建立了弹性结构系统的动力学模型与超声速非定常气动力模型,得到了由微分代数方程表示的多体系统气动弹性动力学方程.通过数值求解微分代数方程的特征值问题,研究了多体系统在平衡位置小扰动运动的稳定性,完成了多体气动弹性系统超声速颤振分析.应用该方法研究了板状翼面及含操纵面翼面的超声速颤振问题,并得到了操纵面处于不同位置时翼面的颤振速度.结果表明,所发展的多体气动弹性系统超声速颤振分析方法,适用于由多个部件组成的工程结构颤振分析.     

Influences of Attack Angles on Aerodynamic Derivatives and Flutter Characteristics of Flat Box Girders北大核心CSCD

以南京第四长江大桥扁平箱梁为研究对象,通过节段模型自由振动风洞试验详细测试了模型在不同风攻角下的颤振响应,探讨了系统非稳态及稳态临界振幅随风速的演化规律.首先,基于颤振响应振幅包络,结合Hilbert变换,识别了系统振幅依存的模态阻尼,并初步阐释了颤振形态随风攻角转变的机理.其次,提取了系统在不同风攻角下的模态参数,基于双模态耦合闭合解法,识别了断面在不同风攻角下的非线性颤振导数,研究了关键颤振导数振幅依存性随风攻角变化的规律及对断面颤振形态和特性的潜在影响.最后,通过逐项拆解模态阻尼,深入剖析了风攻角对非耦合及耦合气动阻尼的影响,并阐明了分项阻尼导致系统颤振性能差异性的动力学机理.     

旋转输液管动力稳定性理论分析北大核心CSCD

基于Lagrange原理和假设模态法建立了旋转输液管的动力学模型.通过降阶升维的方法求解系统的特征值问题,并分析了旋转输液管自由振动特性.得到了不同端部集中质量和转速下,系统特征值随流速升高的演变轨迹.揭示了临界流速随系统参数的变化规律.研究发现,内部流体的流动对旋转输液管动力学特性存在显著影响.在某些参数组合下,系统低阶模态能够形成不同形式的内共振关系.预示了旋转输液管模型蕴含丰富的动力学现象.     

一类新的(2n-1)点二重动态逼近细分

利用正弦函数构造了一类新的带有形状参数ω的(2n-1)点二重动态逼近细分格式.从理论上分析了随n值变化时这类细分格式的C~k连续性和支集长度;算法的一个特色是随着细分格式中参数ω的取值不同,相应生成的极限曲线的表现张力也有所不同,而且这一类算法所对应的静态算法涵盖了Chaikin,Hormann,Dyn,Daniel和Hassan的算法.文末附出大量数值实例,在给定相同的初始控制顶点,且极限曲线达到同一连续性的前提下和现有几种算法做了比较,数值实例表明这类算法生成的极限曲线更加饱满,表现力更强.     

修正AHP中判断矩阵的近似加速算法

作为系统工程中典型的定性定量综合集成方法,层次分析法(AHP)在各种复杂系统综合评价和多目标决策中具有广泛的应用价值.为加速修正AHP中判断矩阵的一致性,通过构造修正特征值与步长函数,得到最优修正步长的近似加速算法(AAA).理论分析和实例分析的初步结果说明:近似加速算法结果稳定,精度高,修正速度快,在系统工程中具有推广应用价值.     

A block-Lanczos method for large continuation problems

The computation of solution paths of large-scale continuation problems can be quite challenging because a large amount of computations have to be carried out in an interactive computing environment. The computations involve the solution of a sequence of large nonlinear problems, the detection of turning points and bifurcation points, as well as branch switching at bifurcation points. These tasks can be accomplished by computing the solution of a sequence of large linear systems of equations and by determining a few eigenvalues close to the origin, and associated eigenvectors, of the matrices of these systems. We describe an iterative method that simultaneously solves a linear system of equations and computes a few eigenpairs associated with eigenvalues of small magnitude of the matrix. The computation of the eigenvectors has the effect of preconditioning the linear system, and numerical examples show that the simultaneous computation of the solution and eigenpairs can be faster than only computing the solution. Our iterative method is based on the block-Lanczos algorithm and is applicable to continuation problems with symmetric Jacobian matrices. This revised version was published online in June 2006 with corrections to the Cover Date.     

Dissipation-induced subcritical flutter in the acoustics of friction

We consider a gyroscopic system under the action of small dissipative and non–conservative positional forces, which has its origin in the models of rotating elastic bodies of revolution in frictional contact such as the singing wine glass or the squealing disc/drum brakes. The spectrum of the unperturbed gyroscopic system forms a spectral mesh in the plane ‘frequency versus gyroscopic parameter’ with double semi–simple purely imaginary eigenvalues at the nodes. In the subcritical range of the gyroscopic parameter the eigenvalues involved into the crossings have the same Krein signature and thus their splitting due to changes in the stiffness matrix, which break the rotational symmetry of the body, cannot produce complex eigenvalues and, therefore, flutter. We establish that perturbation of the gyroscopic system by the dissipative forces with the indefinite matrix can lead to the subcritical flutter instability even if the rotational symmetry is destroyed. With the use of the perturbation theory of multiple eigenvalues we explicitly find the linear approximation to the domain of the subcritical flutter, which turns out to have a conical shape. The orientation of the cone in the three dimensional space of the parameters, corresponding to gyroscopic, damping, and potential forces, is determined by the sign of an explicit expression involving the entries of both the damping and potential matrices. With the use of a time–dependent coordinate transformation we demonstrate that the conical zones of flutter for the original autonomous system coincide with the zones of the subcritical parametric resonance of the rotationally symmetric flexible body with the load moving in the circumferential direction. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

哈密顿系统的稳定性边界

通过采用摄动法对线性哈密顿参数系统的特征值和特征向量进行灵敏度分析,给出了此类系统的稳定性边界的判据,结果表明:具有约当链的系统重特征根对系统的稳定性起至关重要的作用。     

Smooth continuation of solutions and eigenvalues for variational inequalities based on the implicit function theorem

The implicit function theorem is applied in a nonstandard way to abstract variational inequalities depending on a (possibly infinite-dimensional) parameter. In this way, results on smooth continuation of solutions as well as of eigenvalues and eigenvectors are established under certain particular assumptions. The abstract results are applied to a linear second order elliptic eigenvalue problem with nonlocal unilateral boundary conditions (Schrödinger operator with the potential as the parameter).     

Characterization of critical eigenvalues of axial flow engine compressor PDE model

In this article, we study the spectrum of axial flow engine compressor PDE model and show that there are three types of critical eigenvalues which determine the stability of the steady solution of the system depending on the compressor geometry. All of these critical eigenvalues cross the imaginary axis with nonzero speeds as the throttle coefficient increases, leading to Hopf bifurcations of the system.     

DYNAMIC BIFURCATION OF NONLINEAR EVOLUTION EQUATIONS

§1. IntroductionA key problem in the study of problems in mathematical physics and mechanics is tounderstand and predict patterns and their transitions/evolutions. In ?uid mechanics, forinstance, it is important to study the periodic, quasi-periodic, ape…     

Predicting stochastic characteristics of generalized eigenvalues via a novel sensitivity-based probability density evolution method

This paper proposes a novel numerical method for predicting the probability density function of generalized eigenvalues in the mechanical vibration system with consideration of uncertainties in structural parameters. The eigenproblem of structural vibration is presented by first and the sensitivity of generalized eigenvalues with respect to structural parameters can be derived. The probability density evolution method is then developed to capture the probability density function of generalized eigenvalues considering uncertain material properties. Within the proposed method, the probability density evolution equation for the generalized eigenvalue problem is established accounting for the sensitivity of generalized eigenvalues with respect to structural parameters. A new variable which connects generalized eigenvalues to structural parameters is then introduced to simplify the original probability density evolution equation. Next, the simplified probability density evolution equation is solved by using the finite difference method with total variation diminishing schemes. Finally, the probability density function as well as the second-order statistical quantities of generalized eigenvalues can be predicted. Numerical examples demonstrate that the proposed method yields results consistent with Monte-Carlo simulation method within significantly less computation time and the coefficients of variation of uncertain parameters as well as the total number of them have remarkable effects on stochastic characteristics of generalized eigenvalues.     

Numerical stability analysis and computation of Hopf bifurcation points for delay differential equations

We present a numerical technique for the stability analysis and the computation of branches of Hopf bifurcation points in nonlinear systems of delay differential equations with several constant delays. The stability analysis of a steady-state solution is done by a numerical implementation of the argument principle, which allows to compute the number of eigenvalues with positive real part of the characteristic matrix. The technique is also used to detect bifurcations of higher singularity (Hopf and fold bifurcations) during the continuation of a branch of Hopf points. This allows to trace new branches of Hopf points and fold points.     

Evans Function for Lax Operators with Algebraically Decaying Potentials

We study the instability of algebraic solitons for integrable nonlinear equations in one spatial dimension that include modified KdV, focusing NLS, derivative NLS, and massive Thirring equations. We develop the analysis of the Evans function that defines eigenvalues in the corresponding Lax operators with algebraically decaying potentials. The standard Evans function generically has singularities in the essential spectrum, which may include embedded eigenvalues with algebraically decaying eigenfunctions. We construct a renormalized Evans function and study bifurcations of embedded eigenvalues, when an algebraically decaying potential is perturbed by a generic potential with a faster decay at infinity. We show that the bifurcation problem for embedded eigenvalues can be reduced to cubic or quadratic equations, depending on whether the algebraic potential decays to zero or approaches a nonzero constant. Roots of the bifurcation equations define eigenvalues which correspond to nonlinear waves that are formed from unstable algebraic solitons. Our results provide precise information on the transformation of unstable algebraic solitons in the time-evolution problem associated with the integrable nonlinear equation. Algebraic solitons of the modified KdV equation are shown to transform to either travelling solitons or time-periodic breathers, depending on the sign of the perturbation. Algebraic solitons of the derivative NLS and massive Thirring equations are shown to transform to travelling and rotating solitons for either sign of the perturbation. Finally, algebraic homoclinic orbits of the focusing NLS equation are destroyed by the perturbation and evolve into time-periodic space-decaying solutions.