非线性弦振动方程的多辛算法

摘要

本文利用Hamiltonian空间体系的多辛理论研究了非线性弦微小横向振动问题的数值解法。基于Bridges意义下的多辛积分理论，首先构造了非线性弦振动方程的一阶多辛偏微分方程组及其多种守恒律（多辛守恒律、局部能量守恒律和局部动量守恒律），随后构造了一种等价于Box 多辛格式的新隐式多辛格式，最后，运用该多辛格式对非线性弦振动方程进行了数值模拟，并将模拟结果与吕克璞等人得到的解析解进行比较。数值实验结果显示利用本文构造的多辛格式得到的数值解与吕克璞等人得到的解析解非常接近，这说明该多辛格式能够较为精确地模拟非线性弦振动问题，同时数值结果也反映出了多辛方法的两大优点：精确的保持多种守恒律和良好的长时间数值行为。