层状介质中波传播的研究的进展.pdfVIP

苏先樾：层状介质中波传播研究的进展 185 层状介质中波传播研究的进展* 苏先樾 北京大学力学与工程科学系，北京 100871 , 摘 要 本文简要综述了层状介质波传播研究的历史 并着重介绍了新近发展起来的迴传矩 - 阵 射线法。 关键词 层状介质；波传播；传递矩阵；广义射线；迴传矩阵 中图分类号 O347.4 1 引言 从 dAlembert 以来，诸多研究兴趣持续地关注弹性体中的波动现象。层状介质中波传播的 研究由于有在光学, 声学与地球物理学等诸多领域, 乃至在复合材料与压电材料的层合板等结 构的瞬态响应及无损探测等工程应用的广泛背景也十分活跃。 层状介质是由许多平行平面形成,或由一族同心圆柱面或球面形成。在相邻两层的界面上, 边界条件将由连续性条件所代替, 对液体将是两个, 对轴对称或平面应变的弹性问题是四个, 而对一般的弹性问题则是六个。研究层状介质中波的传播,就是要揭示波经过各层界面的反射 或透射在整个层状介质中行进的复杂过程。 层状介质中波传播的研究可追溯到上个世纪 40 年代, Ewing, Worzel 和 Pekeris 1948 年测 量和分析了声波在美国大西洋海岸浅水中的传播。在其发表的论文中，Pekeris 采用由不同厚 度的均匀层组成的层状液体模型来模拟包括砂石底部且沿垂直方向分层的海水, 并进而用稳 态波函数关于径向波数和圆频率的二重积分表出了瞬态波的传播。该二重积分可由两种不同的 过程求出: 一种称为模态法,一种称为射线法。他们对水中声波的传播研究所采用的观测方法, 理论和分析立即就被应用于层状体中弹性波传播的研究, 详见 Ewing,Jardetzky 和 Press(1957), 及 Brekhovskikh( 1960,1980)分别在专著中的综述。后来 Pao 和 Gajew ski(1977), 以及 Kennet (1983)又分别对广义射线理论及射线积分的各种算法作了深入的阐述。 2 传递矩阵 早期,人们通过边界条件和各界面的连续性条件导出传播函数矩阵, 然后导出频率方程并 * 基金项目：国家自然科学基金项目 ;高等学校博士点专项科研基金项目。 186 苏先樾：层状介质中波传播研究的进展 求其特征根, 进而求得问题的解, 这就是所谓的模态法。对于有许多层的介质,频率方程的根 的计算变得十分困难。 对于地震研究和勘探,为模拟地壳,或为模拟其材料性质急剧变化的不 均匀介质, 大量的层数是需要的。对此, Thomson (1950)及 Haskell (1953) 对以上方法做出了 一个重大的改进。他们的方法实质上是把对于具有 n 层的介质的一个大的矩阵简化为 n 个对于 每层的小矩阵的乘积, 大大简化了层状介质频散曲线的计算工作。后来, Gilbert 和 Backus (1966)证明 Thomson–Haskell 法是 Volterra(1887)为解变系数常微分方程组而采用的传播函数 矩阵的一个特殊情形。 3 广义射线理论 另一种求解方法,射线法, 后来发展为分析多层固体中的瞬态波的广义射线理论。它除了 被用于声学与地球物理学领域外, 也同时在弹性动力学中发展起来的。Stokes(1948) 分析了无 界弹性体中由集中力所产生的弹性瞬态波。Lamb (1904),Lapwood(1949), Pekeris(1955) Chao (1960) 对弹性半空间, Knopoff(1958) 和 Davids (1959) 对板, Cagniard (1939) 和 Brekhovskikh (1960) 对相接的两个固体 等均采用射线法进行了弹性瞬态波的分析。在 1957 年 McGraw-