圆弧状凹陷地形在平面SH波入射下的动力响应高频-震灾防御技术.pdfVIP

第 4 卷 第 2 期 震灾防御技术 Vol. 4, No. 2 2009 年 6 月 Technology for Earthquake Disaster Prevention June., 2009 杨彩红，2009. 圆弧状凹陷地形在平面 SH 波入射下的动力响应：高频解答. 震灾防御技术，4 （2 ）：158—166. 圆弧状凹陷地形在平面 SH 波入 射下的动力响应：高频解答1 杨彩红 （中国地震灾害防御中心，北京 100029 ） 摘要 利用波函数的 Fourier-Bessel 级数展开法和外域型 Graf 加法公式，推导了圆弧状凹陷地形 对平面 SH 波二维散射问题的解析解。与现有其他解析解不同的是，该解析解的获得无须求解线 性方程组，从而避免了高频波入射下方程组系数矩阵可能出现的病态问题。因此，既使在入射波 频率很高的情况下，该解析解也能够收敛。利用该解，在一个非常宽的频带范围内讨论了凹陷地 形的几何形状对地面运动的影响。 关键词：平面 SH 波 凹陷地形 散射 解析解 场地效应 引言 规则的凹陷地形在平面 SH 波入射下，其动力响应解析解的研究是一个经典的工程波动 理论课题，该问题的研究对于揭示该类地形的波动效应具有理论与实际意义。Trifunac （1973） 最先利用波函数的 Fourier-Bessel 级数展开法，求解了半圆形凹陷对 SH 波的散射。之后， Wong 和 Trifunac （1974）利用波函数的Mathieu 级数展开法，求解了半椭圆形凹陷对 SH 波 的散射。由于相对半圆形和半椭圆形，具有不同深、宽比的浅圆弧形，能更好地模拟实际凹 陷地形的边界，因此针对浅圆弧形凹陷对 SH 波的散射问题，Cao 和 Lee （1989）提出了近似 的解析解，利用一半径非常大的圆弧近似模拟水平地表，将散射波函数展开成Fourier-Bessel 级数的形式，引入边界条件后，波函数中的待定系数通过求解线性方程组确定。针对同样的 问题，Yuan 和 Liao （1994）提出了精确的级数解，避免了上述大圆弧假定的引入。该解析解 首先构造满足地表零应力边界条件的散射波级数表达式；然后利用外域型 Graf 加法公式并引 入凹陷圆弧边界的零应力条件，将 SH 波散射问题简化为一无穷维代数方程组的求解问题。 针对上述浅圆弧形凹陷对 SH 波的散射问题，现有的解析解（Cao 等，1989；Yuan 等， 1994）均需求解线性方程组，以确定散射波函数的待定系数。由于入射波频率的增大将会导 1 基金项目 国家科技支撑计划（2006BAC13B02）、国家自然科学基金项目和地震科学联合基金项目（A07045 ）资助 [收稿日期] 2009-02-22 [作者简介] 杨彩红，女，生于 1974 年。2006 年毕业于天津大学土木工程系获硕士学位。主要从事地震工程研究。 E-mail: yangc163.com 2 期 杨彩红等：圆弧状凹陷地形在平面 SH 波入射下的动力响应：高频解答 159 致方程组系数矩阵的病态，从而限制了上述解析解适用的频带范围。本文沿用 Yuan 和 Liao （1994）的求解思路，但是在求解过程中，通过选择不同的求解极坐标系，避免了线性方程组 的求解，散射波函数的待定系数具有明确的解析表达式。利用该方法，本文在一个非常宽的 频带范围内给出了所研究问题的解，探讨了级数解的收敛性，并分析了入射波与凹陷的特性 对地面运动的影响。 1 理论公式的推导 凹陷地形的模型如图 1 所示。半空间介质假 定为均匀、线弹性、各向同性，其剪切波速和密 度为β、ρ。凹陷表面为一圆弧，其圆心在O1 点， 半径为 b ，凹陷的半宽为 a ，深度为 h ，圆心 O1 与地表之间的垂直距离为 d 。入射 SH 波的幅值 为 1，圆频率为 ω，入射角为 γ，其位移在坐标系