非线性方程的保辛近似算法

保守体系的微分方程可用Hamilton体系的方法描述,其特点是保辛。两个辛矩阵之和不能保辛,两个辛矩阵的乘积仍是辛矩阵。最常用的小参数摄动法用的是加法,因此对辛矩阵不能保辛。从保辛的角度,要用正则变换。本文针对非线性微分方程,运用自变量坐标变换,对原系统进行变换。由此推导出变换后系统的变分原理。引入Hamilton对偶变量,通过数学变换,得到变系数非线性方程。针对该方程,本文提出了保辛摄动算法。通过数值算例,对不同步长下,保辛摄动法、多尺度摄动法、龙格库塔法和精确解的结果做了比较。数值例题表明,对于非线性方程,本文提出的保辛摄动算法有良好的精度。在步长增大的情况下,保辛摄动保持了良好的稳定性。     

矩阵摄动法在汽车车架模态分析中的应用

本文采用以模态迭加原理为基础的实模态分析技术及初参数优选法对汽车车架的模态参数进行了识别,讨论了振动特征值问题中关于非重特征值和重特征值的矩阵摄动法,提出了利用有弹性元件悬挂的结构振动测试数据来得到自由——自由结构的模态参数的摄动修正方法.文中还给出了一些数值例子来说明此方法的应用,同时得到了一些重要结论.     

基于非局部效应和表面效应的输流碳纳米管稳定性分析简

应用非局部黏弹性夹层梁模型分析双参数弹性介质中输送脉动流碳纳米管的稳定性. 新模型中同时考虑了由管道内、外壁上的薄表面层引起的表面弹性效应和表面残余应力，经典的欧拉梁模型因此通过引入非局部参数和表面参数得到了改进. 用平均法对其控制方程进行求解，得到了管道稳定性区域. 数值算例揭示了纳米材料的非局部效应、表面效应及两个弹性介质参数对管道固有频率、临界流速和动态稳定性的复杂影响，结论可为纳米流体机械的结构设计和振动分析提供理论基础.     

不确定参数闭环控制系统特征值的区间分析

运用区间理论,讨论了区间参数的结构振动控制问题,给出了求解闭环系统区间特征值的一种方法．基于区间参数导出了区间刚度矩阵和质量矩阵,然后利用矩阵摄动理论和区间扩张理论,推导了复区间特征值上下界估计的算法．这些结果是从二阶系统的左右特征向量出发得到的．将该文方法应用到悬臂梁的控制问题,数值结果表明它是有效的．     

黏塑性本构计算的稳定性分析

提出本构方程计算方法的稳定性问题,针对黏塑性本构计算的显式精确算法的稳定性进行分析,发现该算法并非无条件稳定,使用小扰动方法给出了其计算稳定的必要条件,稳定性条件对数值计算中的时间步长提出限制要求。通过有限元算例验证了分析的正确性,计算结果也表明理论推导得到的稳定性公式能够准确预测满足计算稳定性条件要求的最大时间步长与各参数之间关系。     

基于非局部效应和表面效应的输流碳纳米管稳定性分析

应用非局部黏弹性夹层梁模型分析双参数弹性介质中输送脉动流碳纳米管的稳定性.新模型中同时考虑了由管道内、外壁上的薄表面层引起的表面弹性效应和表面残余应力,经典的欧拉梁模型因此通过引入非局部参数和表面参数得到了改进.用平均法对其控制方程进行求解,得到了管道稳定性区域.数值算例揭示了纳米材料的非局部效应、表面效应及两个弹性介质参数对管道固有频率、临界流速和动态稳定性的复杂影响,结论可为纳米流体机械的结构设计和振动分析提供理论基础.     

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

多维磁浮柔性转子控制系统分岔与控制器设计

讨论了多维悬浮柔性转子控制系统局部及全局分岔问题,首先建立了该复杂系统动力学模型,应用中心流形和求规范形综合方法,得到此系统非半简双零特征值问题的规范形及其普适开折,并进一步讨论了此控制系统的分岔 行为（余维二分岔）及稳定性;给出了为实现稳定控制,控制器参数、转子系统结构参数的相互关系及稳定控制域,即给出分岔 参数条件、分岔曲线、转迁集,最后,给出此柔性转子控制系统的数值仿真结果。     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

功能梯度材料扁锥壳的热弹性大变形分析

研究了功能梯度材料扁薄锥壳在横向非均匀升温场中的几何非线性大变形问题.基于von Kármán几何非线性理论推导出了以中面位移为基本未知量的功能梯度扁薄锥壳在横向非均匀热载荷作用下的轴对称大挠度控制方程.采用打靶法数值求解所得非线性常微分方程边值问题,得到了锥壳的大挠度弯曲变形数值解.给出了锥壳的变形与其形状参数、载荷和材料参数等变化的特征关系曲线,分析和讨论了温度参数和材料梯度变化参数对变形的影响.     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

Stability and bifurcation for a flexible beam under a large linear motion with a combination parametric resonance

Stability and bifurcation behaviors for a model of a flexible beam undergoing a large linear motion with a combination parametric resonance are studied by means of a combination of analytical and numerical methods. Three types of critical points for the bifurcation equations near the combination resonance in the presence of internal resonance are considered, which are characterized by a double zero and two negative eigenvalues, a double zero and a pair of purely imaginary eigenvalues, and two pairs of purely imaginary eigenvalues in nonresonant case, respectively. The stability regions of the initial equilibrium solution and the critical bifurcation curves are obtained in terms of the system parameters. Especially, for the third case, the explicit expressions of the critical bifurcation curves leading to incipient and secondary bifurcations are obtained with the aid of normal form theory. Bifurcations leading to Hopf bifurcations and 2-D tori and their stability conditions are also investigated. Some new dynamical behaviors are presented for this system. A time integration scheme is used to find the numerical solutions for these bifurcation cases, and numerical results agree with the analytic ones.     

基于模态缩减和模态实验数据的结构连接参数识别方法

利用在结构系统可测自由度上获得的不完备模态参数和子结构的有限元模型，根据模态缩减理论，建立了识别子结构间连接子结构参数的优化模型，采用逐次二次规划法求解，改善了测试噪声和模态截断误差的影响。该方法识别精度高、收敛速度快、计算量小，便于工程应用。     

Runge–Kutta methods for a semi-analytical prediction of milling stability

On the basis of Runge–Kutta methods, this paper proposes two semi-analytical methods to predict the stability of milling processes taking a regenerative effect into account. The corresponding dynamics model is concluded as a coefficient-varying periodic differential equation with a single time delay. Floquet theory is adopted to predict the stability of machining operations by judging the eigenvalues of the state transition matrix. This paper firstly presents the classical fourth-order Runge–Kutta method (CRKM) to solve the differential equation. Through numerical tests and analysis, the convergence rate and the approximation order of the CRKM is not as high as expected, and only small discrete time step size could ensure high computation accuracy. In order to improve the performance of the CRKM, this paper then presents a generalized form of the Runge–Kutta method (GRKM) based on the Volterra integral equation of the second kind. The GRKM has higher convergence rate and computation accuracy, validated by comparisons with the semi-discretization method, etc. Stability lobes of a single degree of freedom (DOF) milling model and a two DOF milling model with the GRKM are provided in this paper.     

Dynamics of an inertial two-neuron system with time delay

In this paper, we considered a delayed differential equation modeling two-neuron system with both inertial terms and time delay. By analyzing the distribution of the eigenvalues of the corresponding transcendental characteristic equation of its linearized equation, local stability criteria are derived for various model parameters and time delay. By choosing time delay as a bifurcation parameter, the model is found to undergo a sequence of Hopf bifurcation. Furthermore, the direction and the stability of the bifurcating periodic solutions are determined by using the normal form theory and the center manifold theorem. Also, resonant codimension-two bifurcation is found to occur in this model. Some numerical examples are finally given for justifying the theoretical results. Chaotic behavior of this inertial two-neuron system with time delay is found also through numerical simulation, in which some phase plots, waveform plots, power spectra and Lyapunov exponent are computed and presented.     

Second-order sensitivity of eigenpairs in multiple parameter structures

This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.     

A Numerical Method for Wave Propagation in Viscoelastic Stratified Porous Media

This paper presents a numerical method, a transmission matrix method, for the wave propagation in viscoelastic stratified saturated porous media. The wave propagation in saturated media, based on Biot theory, is a coupled problem. In this stratified three-dimensional model we do the Laplace transform for the time variable and the Fourier transform for the horizontal space coordinate. The original problem is transformed into ordinary differential equations with six independent unknown variables, which are only the function of the coordinate of depth. Thus, we get a transmission matrix of the wave problem for each layer. In the process of solution we use numerical method to calculate the eigenvalues and the eigenvectors of the transmission matrices. In the first step of the solution process we can obtain the wave field in the transformed space. The fast Fourier transform (FFT) method is used to do the inverse Laplace and the inverse Fourier transforms to get the solution in the time space. The detailed formulae are derived and some numerical examples are given.     

Several solution methods for the generalized complex eigenvalue problem with bounded uncertainties

The aim of this paper is to evaluate the effects of uncertain-but-bounded parameters on the complex eigenvalues of the non-proportional damping structures. By combining the interval mathematics and the finite element analysis, the mass matrix, the damping matrix and the stiffness matrix were represented as the interval matrices. Firstly, with the help of the optimization theory, we presented an exact solution—the vertex solution theorem, for determining the exact upper bounds or maximum values and exact lower bounds or minimum values of complex eigenvalues of structures, where the extreme values are reached on the boundary of the interval mass, damping and stiffness matrices. Then, an interval perturbation method was proposed, which needs less computational efforts. A numerical example of a seven degree-of-freedom spring-damping-mass system was used to illustrate the computational aspects of the presented vertex solution theorem and the interval perturbation method in comparison with Deif’s method.     

构造满足特征值要求的杂交元应力子空间

通过引进适当参数给出了杂交元应力空间特征值与参数之间的关系,从而可以通过调节参数来构造满足特征值即模态刚度要求的杂交元应力子空间。本在位移元本征应力模式基础上引进调节参数,同时,利用矩阵H对角化方法计算杂交元应力子空间的本征应力模式,然后由此方便有效地计算特征值,从而大大提高了计算效率。本通过建立Q4杂交应力元特征值与参数之间关系说明了这一方法是确定可行的。     

采用黏弹性人工边界时显式算法稳定性条件

黏弹性人工边界是处理无限域波动问题常用的数值模拟方法。采用显式时域逐步积分算法进行计算时,受黏弹性人工边界的阻尼、刚度等影响,人工边界区的稳定性比内部计算域的更严格,尚无明确、实用的稳定性判别准则,这限制了黏弹性人工边界在显式动力分析中的应用。针对二维黏弹性人工边界,利用基于局部子系统的稳定性分析方法和基于传递矩阵谱半径的稳定性判别准则,给出了可代表整体模型局部特征的不同边界子系统的稳定性条件解析解。通过对比分析不同计算区域的稳定性条件及其影响因素,证明了整体模型的稳定性由角点子系统控制。在此基础上,获得了含黏弹性人工边界的整体模型在显示动力计算中的统一稳定性判别准则和简化实用计算方法。在实际应用中,令积分时间步长满足稳定性条件,即可顺利完成整体模型的动力计算。以上研究可为将黏弹性人工边界应用于显式动力计算时积分时间步长的合理选取提供参考。