SH波传播的三角形格子法.pdfVIP

第27卷增刊 防灾减灾工程学报 V01．27．Suppl． 2007年4月 ofDisasterPreventionand Apr．2007 Journal MitigationEngineering SH波传播的三角形格子法’ 刘铁林，张 俏，钟 伟 (沈阳建筑大学土木工程学院，沈阳110168) 摘要：提出了一种弹性介质中SH波传播数值模拟的新方法——三角形格子法。该方法可以用于模拟具有任意形 状自由表面争内部介质交界面的不均匀介质中SH波的传播问题。本文算法可以按照连续介质的形状及内部分界 面任意割分网格，具有网格剖分的灵活性。通过对被研究体的研究，给出了算法的计算会式。与解析解比较，表明本 文算法有良好的精度。给出了具有起伏地表形状非均匀介质中瞬态SH波传播的数值算例，研究了水平地震波作用 下应力集中和位移放大区域的分布情况。本文算法具有程序实现简单、计算量小、计算机内存占用少的特点，可以 用于模拟更大区域范围内任意地表形状和复杂地质构造中SH波的传播问题，为土动力学中波动问题研究的一种 新方法。 关键词：SH波；三角形格子f半空间；非均匀介质；起伏形状自由表面 中图分类号：P315．3文献标识码。A 文章编号；1672—2132(2007)增刊一0011—04 匀地质构造中的SH波传播。与解析解进行了比较， 0 引言 验证了该方法的准确性与有效性。同时，本文给出了 地表具有任意起伏形状的非均匀介质中瞬态SH波 对于SH波传播问题的研究有实验方法、解析 传播的数值算例，研究了水平地震波作用下应力集 方法和数值方法。刘殿魁[1]用复变函数方法研究了 中和位移放大区域的分布情况。 各向异性介质中一圆孔周围SH波散射和动应力集 中问题，陆建飞[23用Fourier变换方法研究了双介质 1数值算法 材料中裂纹产生的SH波的散射问题，刘有军[31用波 函数展开法研究SH波的散射波与动应力集中。数 1．1控制方程 值方法中主要包括有限差分法、有限元法和边界元 一阶应力、位移混合变量弹性波方程： 法。Jean Virieux[43用有限差分方法研究SH波在不 均匀介质中的传播问题，Adams，BrianM．[51用二维 P碧=誓+誓可2蔷十蔷 ㈤‘1) 弹性有限元法研究了弱运动地震响应。近年来边界 元法得到了提高普及，与近似法领域相比，具有降 吒，：产罢 (2)【z) ％2户磊 维、无限远处条件的传播比较容易实现和结果精度 叮。，：∥窑 L3)(3) 高的优点。Semblat，J．F．嘲应用边界元法研究了沉 ％2∥夏 式中口为y方向的位移； 积盆地区域的地震场地效应，Dineva，P．[7】利用边界 元法求解在一个多层的无弹性区域里的地震波传播 P是介质密度； 问题。此外，研究SH波传播问题还有许多方法，比 ∥为介质的剪切模量。 如：FEM和BIR相结合方法、Fourier变换和残数积根据边界以及内部介质形状剖分网格(实线)的 局部示意如图1。在本文算法中，图中三角形被称为 分、Green函数方法、H