奇异Hamilton系统矩阵的精细积分法

常规位移有限元的结构振动方程是n个二阶常微分方程组.采用一般交分原理推导,将结构振动问题引入Hamiltoil体系,将得到2n个一阶常微分方程组.精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果.对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异.本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间.利用Hamiltoil体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解.数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点.

Precise integration method for singular Hamilton matrix

There are n second-order ordinary differential equations(ODE) for structural dynamics of finite element method.It's obtained 2n first-order ODEs when the dynamic system is introduced into Hamilton system through the principle of general variation.The precise integration method is good for solving ODEs.It can give precise numerical results approaching to the exact solution at the integration points when it is applied to linear time-invariant dynamic system.The system matrix will be singular when the structur...