中的传播问题.pdfVIP

第 45 卷 第 4 期 地 球 物 理 学 报 Vol. 45 , No. 4 　　　　　　 　　　　　　 　2002 年 7 月 CHINESE JOURNAL OF GEOPHYSICS 　 July , 20 0 2 ( ) [文章编号] 　0001 - 5733 2002 04 - 0569 - 06 [ 中图分类号] 　P631 Hamilton 体系及弹性波在层状介质 中的传播问题 孙　雁　刘正兴 (上海交通大学工程力学系 ,上海　200030) [摘 　要] 　利用结构力学与最优控制的模拟理论 ,研究弹性波在层状介质中传播的数值计算 ( ) 方法. 将弹性波传播问题导向哈密顿 Hamilton 体系 ,在哈密顿体系中 ,推导出一种新的半解 析单元 ,称之为动力 - 部分杂交元 , 由此导出一套哈密顿体系下的半解析数值计算方法. 本 文给出了该方法在层状正交各向异性材料介质的弹性波传播问题的数值算例 ,分析了一定频 率的弹性波在层状介质中传播时的位移、应力的模式. 计算结果展现了 Hamilton 体系和辛几 何在弹性波传播问题研究的应用前景. [ 关键词] 　分层介质 ,弹性波 ,动力 - 部分杂交元 ,哈密顿体系 ,数值计算. 1 　引　　言 层状介质中的波传播问题是地球物理领域的重要研究课题之一. 该方面的研究已积 累了许多经典的理论研究成果[1 —3 ] . 但由于问题本身的复杂性 ,解析求解几乎不可能. 一般只能寻求数值计算方法来解决. 目前已发展了多种数值计算方法, 如有限差分 法[4 ,5 ] 、有限元法[6 ,7 ]等 ,这些方法已广泛应用于地震波传播的计算研究中. 哈密顿体系理论最早出现在经典的分析动力学中. 随着近代科学理论的发展 , 由于 哈密顿体系能本质地反映许多物理规律 ,正逐步得到人们的进一步重视. 冯康教授首先 用辛几何方法对哈密顿动力体系进行了研究[8 ,9 ] ,并取得了重要成果. 钟万勰教授则将哈 密顿体系引入到弹性力学及波传播问题中 ,导出了一套横向哈密顿算子矩阵的本征向量 展开法[10 ] . 近年来 , 辛几何和哈密顿体系理论开始应用于地球物理研究领域 ,研究成 果[11 —13 ]正逐步显示出该理论体系的重要性. 本文将哈密顿体系导入到弹性波传播问题中. 在哈密顿体系中 ,将位移向量和应力 向量合在一起构成全状态向量空间. 由于哈密顿体系的导入和与之相应的辛几何的应 用 ,在理论体系上进入了一个杂交系统中 ,本文由此提出了一种新的半解析单元 ———动力 - 部分杂交元 ,该单元与系统自然吻合 ,有望得到较好的计算精度 ,并将其应用于层状介 质的弹性波传播问题的数值计算中. [ 收稿日期] 　2001 - 09 - 30 收到 ,2002 - 02 - 12 收到修定稿. ( ) ( ) [ 基金项目] 　国家自然科学基金重点项目 ,上海市科学技术委员会资助项目 00QA14013 . [ 作者简介] 　孙　雁 ,女 ,1965 年生 ,副教授 ,主要从事工程力学、计算力学研究. Email :sunyan07 @sh163c. sta. net. cn