地震危险性分析的波动方法研究及其应用-固体力学专业论文.docx

硕士研究生学位论文 摘要 I 摘 要 场地的地震危险性分析，为采取各种工程措施防止或减少地震灾害提供依据，对抗 震减灾至关重要。地震时应力波引起地面震动。波的传播问题是力学领域极为活跃的前 沿课题之一，在理论研究和工程应用中都有重要意义。地震波在地球介质中传播时将产 生能量损耗，用粘弹性介质模型描述地球介质比弹性介质模型更接近于实际情况。在粘 弹性介质模型中，标准线性固体模型能较全面地反映岩石地层的粘弹性特征。粘弹性波 的解析解不易求得，而数值方法随着计算机及相关技术发展起来，不仅适用于求解线性 问题，也适用于求解非线性问题，因此需要发展行之有效的数值计算方法。有限差分法 和有限元法是研究粘弹性波动问题的主要数值方法，有限差分法较难处理任意起伏自由 表面问题和内部交界面问题，有限元法计算量大，效率低，处理大范围区域中的波动问 题时有一定难度。因此，探索简单实用的波动数值计算方法具有重要意义。 本文提出一种能够处理考虑体积应变率的标准线性固体粘弹性非均匀介质中波动传 播的三角形格子法。将三维的考虑体积应变率标准线性固体的本构方程退化，给出二维 平面应变和平面应力问题的本构方程。本文算法将实际介质计算区域剖分成三角形网格， 由此构造出控制体，给出积分形式的控制体运动方程；交替运用控制体运动方程和考虑 体积应变率标准线性固体的微分型本构方程，在时间域上递推计算给出被研究区域内各 节点的加速度、速度、位移和各网格的应力。 对考虑体积应变率的标准线性固体介质中的波动传播进行研究，给出数值算例。数 值结果与解析解和特征线解进行比较，验证了方法的计算精度和有效性。 针对通过北川县城、考虑地形起伏的四个地质剖面，结合地质参数和 Taft 地震波数 据，对这四个地质剖面内的波动传播进行了数值模拟。研究了最大主应力、水平加速度、 水平位移以及最大剪应力的动态集中现象。利用地表处水平加速度的时程响应曲线，给 出频谱分析和加速度反应谱，确定地震波入射时地表附近的危险区域，进而对地震诱发 的山体和实际场地的灾害性分布进行了研究。 本文提出的粘弹性波动数值模拟方法，为研究大规模、不规则自由边界非均匀介质 中的波动传播问题提供了一种新的数值方法，为地震危险性分析和地震小区划的研究提 供了一种工具。 关键词：地震危险性；波动；体积应变率；标准线性固体；三角形网格；动态集中 硕士研究生学位论文 Abstract III Abstract Seismic risk analysis is very important for mitigating earthquake hazard. The ground motion is caused by stress wave during an earthquake. Wave propagation is one of the important researching fields in mechanics and is of significance in theory and engineering. The viscoelastic medium model is appropriate for seismic wave propagation in earth medium. Among all the viscoelastic medium models, standard linear solid model can give relatively comprehensive viscoelastic properties of lithosphere. The analytical solution of viscoelastic wave is difficult to obtain. N umerical methods are used widely to solve both linear problems and nonlinear problems with the developing of computer and calculating technology. Finite difference method and finite element method are two main numerical methods for simulating viscoelastic wave problems. Finite difference method can not deal well with the irregular free surface and internal interface. Finite element method requires large co