基于祖冲之类方法具有保辛性

分析微分-代数方程的祖冲之类算法可以在时间格点上精确满足约束方程,计算格式简单,易于编程计算,因此十分实用。本文证明,基于祖冲之类方法不仅保约束,同时还保辛。     

基于辛本征空间的线性阻尼振动系统动力学分析

对于考虑阻尼项和陀螺项的一般线性动力学振动系统,建立基于辛本征空间展开求解的一般方法.基于Rayleigh商本征值的模态展开方法被广泛应用于复杂结构动力系统振动分析,但对于很多机械系统,由于其不能有效考虑陀螺效应的影响,其适用性却受到很大限制.该文首先讨论了无阻尼系统Rayleigh商本征值问题与辛本征值问题的对应关系,表明前者实际可由后者的一种退化形式给出（也即忽略陀螺效应）,而后者更具有一般性.在此基础上,进一步基于辛本征空间本征向量展开,推导了同时考虑阻尼和陀螺系统的一般线性动力学系统的有效求解方法.数值算例选取不考虑陀螺效应及考虑陀螺效应的两种线性阻尼振动系统对所提出的方法进行了验证,分析结果表明了该文所建立方法的正确性和有效性.     

求解周期性分段线性系统动态响应的高效数值方法

基于参变量变分原理,提出了一种求解具有大量间隙弹簧的周期性分段线性系统动态响应的高效率数值方法.通过参变量变分原理来描述间隙弹簧,将复杂的非线性动力问题转化为线性互补问题求解,避免了求解过程中的迭代和刚度阵更新,该算法能准确判断间隙弹簧的压缩和松弛状态.基于结构的周期性和能量传播速度的有限性,提出了一种求解系统动态响应的高效率精细积分方法.该算法指出周期结构的矩阵指数中存在大量的相同元素和零元素,从而不需要重复计算和存储这部分元素,节省了计算量并降低了计算机存储要求.分析了一个五自由度分段线性系统在简谐荷载作用下的动力学行为,包括稳定的周期运动、准周期运动和混沌运动.通过与Runge-Kutta方法的比较,该文方法的正确性和高效率得到了验证.     

一维离散结构能带结构与表面态的辛分析方法

本文基于辛几何方法推导了一维离散周期结构、半无穷周期结构和含杂质半无穷周期结构的本征方程,力求建立一个完整的辛分析体系。通过辛分析,将一维离散半无限周期结构转化到一个元胞上求解,大大简化了计算量。对于含杂质半无穷周期结构,结合辛分析和W-W算法,给出求解含杂质半无穷周期结构本征值问题的精确、稳定和高效算法。数值算例说明了本文算法的有效性。     

输入受限LQ控制的参变量变分原理和算法

在最优控制理论中根据模拟理论思想发展了塑性力学和接触力学中的参变量变分原理, 并建立了控制输入受限的线性二次(linear quadratic, LQ)最优控制问题的求解新方程---耦合的Hamilton正则方程与线性互补方程. 通过将连续时间离散成一系列等间距时间区段, 在离散时域内采用参数二次规划方法给出数值求解输入受限的LQ最优控制问题的新算法. 数值仿真验证了该算法在求解控制输入受限的LQ最优控制问题中的有效性, 并且该算法具有较快的收敛性, 在大步长下具有较高的计算精度.      

Improved precise integration method for differential Riccati equation

An improved precise integration method(IPIM) for solving the differential Riccati equation(DRE) is presented.The solution to the DRE is connected with the exponential of a Hamiltonian matrix,and the precise integration method(PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix.The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE.Based on the error analysis,the criterion for choosing two parameters of the PIM is given.Three kinds of IPIMs for solving the DRE are proposed.The numerical examples show that the IPIM is stable and gives the machine accuracy solutions.     

基于哈密顿动力系统新变分原理的保辛算法之二：算法保辛性质证明

文献[1]给出了哈密顿系统的一个新的变分原理,并基于此变分原理,通过选择一个时间步长两端不同广义位移或广义动量为独立变量,给出了四种不同类型的求解哈密顿动力系统的数值方法。本文将分别证明这四类数值方法都是保辛的数值方法。     

无网格模型的子结构分析

提出一种无网格模型的子结构分析方法。基于界带分析方法,处理无网格模型子结构之间的连接,能够保证各个子结构之间无缝连接,从而使无网格子结构模型与原无网格模型等价,同时还能够降低无网格刚度方程的计算规模。数值算例表明,本文方法是正确可靠的,可以保证各个子结构之间无缝连接。     

线性时变系统二次最优控制问题的保辛近似求解

状态空间的最优控制体系是保守的,其近似算法应当保辛．提出了基于分段常值精细积分方法的保辛摄动近似方法,在同一框架下求解了线性时变LQ最优控制中的计算问题,即变系数矩阵Riccati方程和状态反馈方程．该算法是保辛的,具有很好的数值稳定性和精度．算例验证了算法的有效性．     

On the stability analysis of plates and shells using a quadrilateral, 16-degrees of freedom plat shell element DKQ16

The linear buckling problems of plates and shells were analysed using a recently developped quadrilateral, 16-degrees of freedom flat shell element ( called DKQ16 ). The geometrical stiffness matrix was established. Comparison of the numerical results for several typical problems shows that the DKQ16 element has a very good precision for the linear buckling problems of plates and shells.