层状场地中凹陷地形Rayleigh波斜入射下三维地震响应分析.docVIP

层状场地中凹陷地形Rayleigh波斜入射下三维地震响应分析 　　摘要： 采用以移动斜线荷载动力格林函数为基本解的2.5维间接边界元方法求解了层状场地中凹陷地形对斜入射Rayleigh波的三维散射问题。方法将总波场分解为自由波场和散射波场，自由波场是通过将层状场地三维精确动力刚度矩阵行列式取为零形成频散方程，采用搜索迭代方法求得频散曲线，进而求得自由波场。散射波场则通过施加在凹陷表面各离散斜线单元上的虚拟移动均布斜线荷载产生的动力响应来模拟，虚拟荷载的密度可通过凹陷表面的零应力边界条件建立方程求得。通过与二维散射结果的比较验证了方法正确性，并以均匀场地和基岩上单一土层场地中凹陷地形为模型，分别在频率和时域内进行了计算分析，研究了Rayleigh波斜入射下凹陷地形周围三维散射的基本规律。研究表明凹陷地形对Rayleigh波的三维散射取决于凹陷的界面形状、斜入射角度、入射频率以及层状场地的频散和多模态特性，层状场地中存在Rayleigh的高阶模态在某些频率处的位移幅值会显著大于其他模态及均匀场地位移幅值的现象。 　　关键词： 凹陷地形； Rayleigh波； 三维散射； 层状场地； 格林函数 　　中图分类号： P315.9； TU 311.3 文献标志码： A 文章编号： 1004-4523（2015）05-0809-13 　　引 言 　　局部地形对地震动的幅值及其空间分布有着显著影响，这已在多次的强震观测和震害调查中得以证实。同时近年来在中国长大桥梁、超高大坝及高速铁路等大型重要工程不断涌现，其中许多工程位于场地条件复杂的高烈度地震区，这些工程的规划建设及抗震设防亟需精确可靠的设计地震动参数，因此研究局部地形对地震动的影响有着显著的理论意义和工程参考价值。 　　凹陷地形作为常见的局部地形，其对地震波的散射问题是一直是国内外学者关注的热点问题，自Trifunac[1]开创性地给出了凹陷地形对SH波的散射解析解后，诸多学者针对该问题进行了研究。现有研究包括凹陷地形对SH波的二维平面外散射[14]，凹陷地形对P，SV和Rayleigh波的二维平面内散射[514]，二维凹陷地形对斜入射地地震波的三维散射[1519]和三维凹陷地形对地震波的散射[2024]。 　　地震波包括在地球内部传播的体波（SH，P和SV）和沿地表传播的面波。对于远场，相对于体波，面波振幅较大，波长较长，携带能量较多，是造成震害的主要原因。然而分析以上文献发现，关于凹陷地形对地震波散射的研究主要是针对体波，关于面波散射的研究非常有限，尤其是层状场地Rayleigh波的传播存在频散和多模态特性，目前鲜有研究，仅文献[6]针对层状场地中凹陷地形对Rayleigh波的二维散射给出了少量结果，关于层状场地中凹陷地形对斜入射Rayleigh波的三维散射，则至今未有研究。 　　针对上述问题，本文考虑天然土体成层特性，同时考虑Rayleigh波斜入射下凹陷地形的三维地震反应特性，在文献[25]给出层状场地三维精确动力刚度矩阵的基础上，推导层状半空间中移动斜线荷载动力格林影响函数，并以该格林函数为基本解，建立一种新的2.5维间接边界元方法求解了层状场地中凹下地形对斜入射Rayleigh波的三维散射。该方法的特色在于荷载可直接加在边界上而不会引起奇异性，保证了该方法的精度及对复杂边界条件的适应性，同时采用层状半空间动力格林函数，求解当中无须离散自由地表和土层交界面，最大限度降低了求解的自由度。 　　1 计算模型和相应公式 　　如图1所示，截面形状沿y轴保持不变的无限长凹陷地形位于层状半空间中。层状半空间有任意层相互平行的水平土层和下卧弹性基岩组成。入射波为Rayleigh面波，相速度速为ca，圆频率为ω，入射方向与y轴成夹角为ψ（ψ=90°为垂直入射，对应二维散射情况）。在随时间简谐变化的Rayleigh波斜入射下，由于入射方向与凹陷地形轴线不垂直，凹陷地形周围的地震反应是三维的，但该三维地震反应在凹陷地形的任意两个截面上完全相同，仅因沿y轴的位置不同而相差一个相位，因此可仅取一个截面进行离散求解，得到计算截面产生的动力响应后，其余截面产生动力响应，可依据其位置和Rayleigh波视速度，通过将计算截面动力响应偏移相应相位得到，最后将所有截面结果求和，也即沿凹陷地形轴线将所有截面动力响应积分，求得最终问题的解。 　　3.3 时域结果 　　图12首先给出了均匀场地Rayleigh波入射角度为ψ=45°和90°时凹陷附近地表位移时程。u，v和w分别表示沿x，y和z方向的位移时程序，而u0，v0和w0分别表示自由场表位移时程。从图11中可以看出，由于凹陷地形对Rayleigh波的散射，使得凹陷左侧位移幅值相对于自由场有着显著的放大，这与频域中结果一致。但凹陷地