hamilton弹性动力学及其辛算法-中国学术会议在线.pdfVIP

中国力学学会学术大会2005 ( CCTAM2005 ) Hamilton 弹性动力学及其辛算法∗ 1 1,2 1,3 罗恩 朱慧坚 黄伟江 (1.中山大学应用力学与工程系，广东 广州 510275; 2.佛山科学技术学院土木工程系, 广东 佛山 528000; 3.广州市建筑科学研究院，广东 广州 510030) Lagrange 力学最核心部分是Lagrange 方程和Hamilton 变分原理, Lagrange 方程是二阶微分方程组，对理论分析与数 值计算十分不便，而Hamilton 变分原理只能反映动力学的时间两端值问题。为此，本文建立了能反映动力学初值问题的非传 统Hamilton 变分原理，其泛函式如下： T T ο ο ~ t1 ~T T T Π(q) ∫ L (q, q)dt - m q q − (q − q )mq + m(q − q )q 0 1 1 0 0 0 0 0 1 式中m 为质量矩阵，带顶标ο的量为限制变分量。 Hamilton 力学最核心部分就是在相空间上的Hamilton 正则方程和相空间Hamilton 变分原理。Hamilton 正则方程是一阶 微分方程组，形式简单并呈现出美妙的对称性，而且还具有自然辛结构。这些独特的优点为定性研究和求解提供了新途径。 但是，对其具有的优越计算功能，人们相当长时间内没有察觉到。直到20 世纪80 年代和90 年代，冯康和钟万勰才分别提出 基于Hamilton 正则方程的辛几何算法和精细时程积分法。这些都是很有特色的优质算法。 相空间Hamilton 变分原理早在 1892年由法国数学家Poincaré提出,至今已有 100多年,一直没有受到应有的重视。国内 外对其论述的文献相当少,国内 1985年出版的中国大百科全书(力学卷)中也没有这种变分原理的条目。而对其具有的更为优 越的计算性能，国内外一直没有人察觉到。这种变分原理具有自然辛结构，但它只反映动力学的时间两端值问题。为此，作 者最近建立了能反映动力学初值问题的非传统相空间Hamilton 变分原理，其泛函如下： T T ο ο t1 T ~ T ~T Π2 (p , q) ∫0 [p q − H(p , q)]dt − p 1 q1 − (q0 − q0 )p 0 + (p 0 − p 0 )q1 并且还开展了基于这种非传统相空间Hamilton 型变分原理的辛算法研究。 从几何上看，Lagrange 力学对应于Riemann 几何，因此Lagrange 力学可称为Riemann 几何的力学；而Hamilton 力学 对应于辛几何，因此Hamilton 力学可称为辛几何的力学。 动力学体系本质上应当是Hamilton 体系，因此动力学的基本理论和分析方法要完全在Hamilton 体系下来建立，才是最 合理和最完美的选择；但是，100多年来动力学的基本理论和分析方法都是在Lagrange