欢迎光临《计算力学学报》官方网站!
孙雁.奇异Hamilton系统矩阵的精细积分法[J].计算力学学报,2009,26(1):46~51
本文二维码信息
码上扫一扫!
奇异Hamilton系统矩阵的精细积分法
Precise integration method for singular Hamilton matrix
投稿时间:2007-05-12  
DOI:10.7511/jslx20091008
中文关键词:  精细积分  Hamilton奇异矩阵  共轭辛正交  非齐次方程
英文关键词:precise integration method  singular Hamilton matrix  conjugate sympletic orthogonal normalization  non-homogeneous dynamic systems  complete pivot gauss-jordan elimination
基金项目:国家自然科学重点基金(10632032);国家自然科学基金(10672100)资助项目.
作者单位
孙雁 上海交通大学 船舶海洋与建筑工程学院 工程力学系,上海 200030 
摘要点击次数: 1626
全文下载次数: 1519
中文摘要:
      常规位移有限元的结构振动方程是n个二阶常微分方程组。采用一般变分原理推导,将结构振动问题引入Hamilton体系,将得到2n个一阶常微分方程组。精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果。对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异。本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间。利用Hamilton体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解。数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点。
英文摘要:
      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 structure has rigid body displacement for non-homogeneous dynamic systems. A new method named complete pivot Gauss-Jordan elimination is proposed. It is used to derive to zero eigen-solutions for singular matrix. Based on this method, it is easy to separate the subspace corresponding to zero eigen-solutions from singular Hamilton matrix by using conjugate sympletic orthogonal normalization between Hamilton eigen-vectors. Then the singular portion can be excluded through projection. The singular solution is derived from analysis. The nonsingular solution is derived from the precise integration method. The numerical result demonstrates the validity and efficiency of the method.
查看全文  查看/发表评论  下载PDF阅读器
您是第13579085位访问者
版权所有:《计算力学学报》编辑部
本系统由 北京勤云科技发展有限公司设计