基于黎曼几何的非线性振子动力学分析递推解析算法

来源：计算力学学报杂志2019年第3期北京时间：

作者：杨喆,陈国海,杨迪雄

单位：大连理工大学 工程力学系 工业装备结构分析国家重点实验室, 大连, 116024,大连理工大学 工程力学系 工业装备结构分析国家重点实验室, 大连, 116024,大连理工大学 工程力学系 工业装备结构分析国家重点实验室, 大连, 116024

摘要：基于黎曼几何和变分原理，推导了黎曼流形上非线性耗散动力系统的二阶微分动力学方程，并运用流形收缩的概念将动力学方程离散化，进而建立了相应的递推求解格式。选取3个自治非线性阻尼振子系统，分别采用递推解析算法和龙格库塔法求解微分动力学方程，并比较分析了不同的时间步长下两种算法的计算耗时。结果表明，与龙格库塔法相比，基于黎曼几何的递推算法不仅能得到每一时步的解析表达式，而且计算耗时短，计算效率高。基于黎曼流形的动力学方程递推算法为非线性动力学系统的解析求解提供了新思路。

关键词：非线性振子,黎曼几何,递推解析算法,龙格库塔法,计算效率,

基金资助：国家自然科学基金（51478086；11772079）；土木工程防灾国家重点实验室开放基金（SLDRCE17-03）资助项目.