圆弧状凹陷地形在平面SH波入射下的动力响应高频解答.pdfVIP

第 4卷 第 2期 震灾防御技术 、，01．4．NO．2 2009年6月 TechnologyforEarthquakeDisasterPrevention June．，2009 杨彩红，2009．圆弧状凹陷地形在平面SH波入射下的动力响应：高频解答．震灾防御技术，4(2)：158—166 圆弧状凹陷地形在平面 SH波入 射下的动力响应：高频解答 杨彩红 (中国地震灾害防御中心 ，北京 100029) 摘要 利用波函数的Fourier-Bessel级数展开法和外域型Graf加法公式，推导了圆弧状凹陷地形 对平面SH波二维散射问题的解析解。与现有其他解析解不同的是，该解析解的获得无须求解线 性方程组，从而避免了高频波入射下方程组系数矩阵可能出现的病态问题。因此，既使在入射波 频率很高的情况下，该解析解也能够收敛。利用该解，在一个非常宽的频带范围内讨论了凹陷地 形的几何形状对地面运动的影响。 关键词：平面 SH波 凹陷地形 散射 解析解 场地效应 引言 规则的凹陷地形在平面 SH波入射下，其动力响应解析解的研究是一个经典的工程波动 理论课题，该问题的研究对于揭示该类地形的波动效应具有理论与实际意义。Triftmac(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．tom 2期 杨彩红等：圆弧状凹陷地形在平面 SH波入射下的动力响应：高频解答 159 致方程组系数矩阵的病态，从而限制了上述解析解适用的频带范围。本文沿用 Yuan和 Liao (1994)的求解思路，但是在求解过程中，通过选择不同的求解极坐标系，避免了线性方程组 的求解，散射波函数的待定系数具有明确的解析表达式。利用该方法，本文在一个非常宽的 频带范围内给出了所研究问题的解，探讨 了级数解的收敛性，并分析了入射波与凹陷的特性 对地面运动的影响。 0 1 理论公式的推导 凹陷地形的模型如图1所示。半空间介质假 定为均匀、线弹性、各向同性，其剪切波速和密 度为 、P。凹陷表面为一圆弧，其圆心在 O】点， 半径为b，凹陷的半宽为 a，深度为 h，圆心 O1 与地表之问的垂直距离为 d。入射 SH波的幅值 为 1，圆频率为∞，入射角为 ，其位移在坐标系 图 1凹陷地形的模型 一 O2-Z2中具有如下表达式：