弹性杆动力学方程的高精度数值计算及仿真.pdf

第9卷 第2期 2009年 1月 科 学 技 术 与 工 程 Vo1．9 No．2 Jan．2009 1671·1819(2009)2—0246一o4 ScienceTechnologyandEngineering ⑥ 2009 Sci．Tech．Engng． 力 学 弹性杆动力学方程的高精度数值计算及仿真 杨 斌 赵维加 黄健飞 贾美娟 (青岛大学数学科学学院，青岛266071) 摘 要 对于描述运动弹性杆的非线性偏微分／代数方程组，利用三角谱方法离散弧长变量s，将方程组离散为常微分／代数 方程组，然后利用谱延迟修正方法进行数值求解，使数值离散方法达到谱精度。还分析了弹性杆运动的数值仿真I4题并给出 了数值结果。 关键词 弹性杆 动力学方程 谱方法 数值计算 中图法分类号 O313．3 0242．2； 文献标志码 A 弹性杆是一种重要的力学模型，广泛应用于 由杆的中心线决定。在惯性坐标系O一 下，设 s 描述细长的结构如海底 电缆、钻杆等，有关 的研 和t分别为弧长和时间变量，杆的中心线