层状弹性介质波.docVIP

博士生综合能力考试读书报告 题目： 层状介质中弹性波传播 问题的研究概述 姓 名： 马宇立 学 号： 院 系： 力学与工程科学系 专 业： 固体力学 导 师： 苏先樾 教授 一、研究背景 1.1 工程背景 近几十年以来，随着科学技术的进步、工业的发展，工程结构正在向大型化、复杂化、智能化方向飞速地发展。大量先进技术和新材料在航空和航天领域得到广泛应用，如碳纤维增强复合材料已经应用在空间站的外壳、天线等大型部位；2002年，洛克希德·马丁公司联合攻击战斗机美国空军、海军、海军陆战队通用的新型隐形、超音速、轻型多用途战斗攻击机Ewing，Worzel和 Pekeris[1] 1948年测量和分析了声波在美国大西洋海岸浅水中的传播。在其发表的论文中，Pekeris采用由不同厚度的均匀层组成的层状液体模型来模拟包括砂石底部且沿垂直方向分层的海水，并进而用稳态波函数关于径向波数和圆频率的二重积分表出了瞬态波的传播。该二重积分可以由两种不同的过程求出：一种称为模态法，一种称为射线法。 在此之后，这种对水中声波的传播研究所采用的观测方法、理论和分析立即就被应用于层状体中弹性波传播的研究。层状介质是由许多平行平面形成，或由一族同心圆柱面或球面形成。在相邻两层的界面上，边界条件将由连续性条件所代替，对液体将是两个，对轴对称或平面应变的弹性问题是四个，而对一般的弹性问题则是六个。研究层状介质中波的传播，就是要揭示波经过各层界面的反射或透射在整个层状介质中行进的复杂过程。 早期，人们通过边界条件和各界面的连续性条件导出传播函数矩阵，然后导出频率方程并求其特征根，进而求得问题的解，这就是所谓的模态法。对于有许多层的介质，频率方程的根的计算变得十分困难。对于地震研究和勘探，为模拟地壳或其他材料性质急剧变化的不均匀介质，大量的层数是需要的。对此，Thomson[2](1950)及Haskell[3](1953) 对以上方法做出了一个重大的改进。他们的方法实质上是把对于具有n层的介质的一个大的矩阵简化为n个对于每层的小矩阵的乘积，大大简化了层状介质频散曲线的计算工作。 另一种求解方法，射线法，后来发展为分析多层固体中的瞬态波的广义射线理论。它除了被用于声学与地球物理学领域外，也同时在弹性动力学中发展起来的。Stokes(1948)分析了无界弹性体中由集中力所产生的弹性瞬态波。Lamb (1904)，Lapwood(1949)，Pekeris(1955)和Chao(1960)对弹性半空间，Knopoff(1958)和Davids(1959)对板，Cagniard (1939)和Brekhovskikh(1960)对相接的两个固体等均采用射线法进行了弹性瞬态波的分析。在1957年McGraw-Hill出版的《层状介质中的弹性波》这一名著中，W.M.Ewing，W.S. Jardetzky和F. Press对这些研究状况作了完美的综述。在该书中，他们将柱波或球波的Weyl-Sommerfeld积分的级数展式视为一种可能的数学表示，但他们并没有期望有富有成效的应用。按照广义射线理论，在层状介质中经多次反射或透射的弹性波可以视为由一组射线积分表示的沿不同路径传播的弹性波。后来人们发现这组积分可以用Cagniard 方法严格求解。Mulle[4](1969,1970)据此求解了由许多平行层组成的介质中弹性瞬态波的传播问题。Gilbert和Helmberger(1972)，Chen(1977)，Pao和Ceranoglu(1977)相继对球和圆柱面界定的介质，Chapman(1974,1976)对非均匀介质的波动问题进行了分析。特别值得指出的是，Pao和Gajewski[5](1977)重新阐述了波在两个不同介质的弹性体间的界面上折射这一边值问题，并且由Laplace变换与Fourier变换的组合严格地推导出平面波的级数展式。由此，他们系统地构造了广义射线积分和识别各类点源的源发射函数，导出了在内部与表面观察点处对于位移、速度、加速度和应力的接收函数，沿发射源与接收点之间几何上确定的所有可能路径传播的平面波的相函数，以及所研究的位移势的反射系数，透射系数等等。对这些积分施加变换, 根据Laplacce逆变换将给出复幔度平面内的有限数值积分，在层状体中按波在给定的源-接收器之间传播所达到的时间顺序计算这些积分。人们提出了各种数值算法来求关于稳态波函数的双重积分即广义射线积分,Kennet[6](1983)对此有详细的讨论。 近年来，Pao 等又发展了一种求解弹性波传播的新方法——回传矩阵-射线法。这种方法最初是针对刚架结构中波的传播而提出的(