不可压饱和多孔弹性杆动力响应的多辛方法.PDFVIP

　 ， ３６ ３ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１５年３月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３６，Ｎｏ．３，Ｍａｒ．１５，２０１５ 文章编号：１０００⁃０８８７（２０１５）０３⁃０２４２⁃１０ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 不可压饱和多孔弹性杆动力 响应的多辛方法∗ １，２ １ １ 刘雪梅 ，　 邓子辰 ，　 胡伟鹏 （１． 西北工业大学 力学与土木建筑学院，西安 ７１００７２； （２． 长安大学 工程力学系，西安 ７１００６４） （我刊编委邓子辰来稿） 摘要：　 研究了不可压饱和多孔弹性杆的一维动力响应问题．基于多孔介质理论，在流相和固相微 观不可压、固相骨架小变形的假定下，建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数 学模型．利用Ｈａｍｉｌｔｏｎ空间体系的多辛理论，构造了不可压饱和多孔弹性杆轴向振动方程的多辛形 式及其多种局部守恒律．采用中点Ｂｏｘ离散方法得到轴向振动方程的多辛离散格式和局部能量守 恒律以及局部动量守恒律的离散格式；数值模拟了不可压饱和多孔弹性杆的轴向振动过程，记录 了每一时间步的局部能量数值误差和局部动量数值误差．结果表明，已构造的多辛离散格式具有很 高的精确性和较长时间的数值稳定性，这为解决饱和多孔介质的动力响应问题提供了新的途径． 关　 键　 词：　 饱和多孔弹性杆；　 多辛方法；　 动力响应；　 离散 中图分类号：　 Ｏ３４３　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１５．０３．００２ 引　 　 言 近年来，随着饱和多孔介质在地震工程学、土力学、生物力学等领域的广泛应用，其相关理 论和研究已经越来越多地引起学术界的重视和关注，尤其有关各类饱和多孔结构流固耦合动 ［１］ 静力响应分析的研究更是取得了很大的成果．Ｔｈｅｏｄｏｒａｋｏｐｏｕｌｏｓ等 研究了多孔弹性板的弯曲 ［２］ 振动问题；Ａｎｋｅ等 基于Ｂｉｏｔ理论，建立了饱和多孔弹性Ｍｉｎｄｌｉｎ板的动力学数学模型，并给 ［３］ 出以挠度、转角和孔隙应力为基本未知量的虚功原理；杨骁等 基于多孔介质理论，建立了饱 和多孔弹性杆件弯曲变形时动力响应的数学模型，并利用Ｌａｐｌａｃｅ变换，给出了饱和多孔悬臂 梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等物理量随时间的响应曲线．然而，对于饱和 多孔结构流固耦合动静力响应问题在数值求解过程中如何保持系统的局部能量守恒律和局部 动量守恒律等多种局部几何性质，学术界却很少讨论． 众所周知，我国数学家Ｆｅｎｇ（冯康）先生于上世纪八十年代提出了基于辛几何原理计算 ［４］ Ｈａｍｉｌｔｏｎ系统的Ｈａｍｉｌｔｏｎ辛算法 ，这种算法有一个优于其它算法的重要特点，就是它可以对 守恒量长时间数值模拟，但是当用辛算法对无限维 Ｈａｍｉｌｔｏｎ 系统进行离散时就具有局限性， ∗ 收稿日期：　 ２０１４⁃１２⁃０２；修订日期：　 ２０１４⁃１２⁃２６ 基金项目：　 国家自然科学基金１１１７２２３９；中央高校基金（２０１４Ｇ１１２１０９６） 作者简介：　 刘雪梅（１９８０—），女，陕西榆林人，讲师，硕士（通讯作者． Ｅ⁃ｍａｉｌ：ｌｉｕｍｅｉ１１２＠１６３．ｃｏｍ）． ２４２ 刘　 雪　 梅