[论文解读] Normal-fault stress and displacement through finite-element analysis

Antonietta Megna, S. Barba+1 more|ArXiv.org|Oct 17, 2005

earthquake and tectonic studies参考文献 15被引用 8

一句话总结

本文提出一种有限元方法，通过迭代力校准方案模拟正断层引起的应力场和位移场。其结果与解析解高度一致——位移差异小于5%，方向误差小于15°，应力误差小于35%，证明该方法在地质力学研究中模拟平滑滑动分布方面具有有效性。

ABSTRACT

We compute displacement and stress due to a normal fault by means of two-dimensional plane-strain finite-element analysis. To do so, we apply a system of forces to the fault nodes and develop an iterative algorithm serving to determine the force magnitudes for any slip distribution. As a sample case, we compute the force magnitudes assuming uniform slip on a 10-km two-dimensional normal fault. The numerical model generates displacement and stress fields that compare well to the analytical solution. In fact, we find little difference in displacements (<5%), displacement orientation (<15 DEG), and stress components (<35%, half of which due to slip tolerance). We analyze such misfit, and discuss how the error propagates from displacement to stress. Our scheme provides a convenient way to use the finite-elements direct method in a trial-and-error procedure to reproduce any smooth slip distribution.

研究动机与目标

• 开发一种数值方法，利用有限元分析精确模拟正断层滑动引起的应力场和位移场。
• 解决在数值模型中对断层面施加合适剪切力以重现所需滑动分布的挑战。
• 通过二维正断层均匀滑动的解析解验证有限元模型的准确性。
• 量化并分析数值结果与解析解之间的差异，特别关注从位移到应力的误差传播。
• 提供一种实用的迭代框架，利用直接有限元法重现任何平滑滑动分布的断层建模。

提出的方法

• 采用二维平面应变有限元分析，模拟长度为10公里的正断层及其均匀滑动。
• 在断层面施加节点力系以模拟滑动，力的大小通过迭代确定。
• 使用迭代算法调整力的幅值，直至所得滑动分布与目标（如均匀滑动）相匹配。
• 通过与解析解比较计算得到的位移场和应力场以进行验证。
• 通过分析位移幅值、方向及应力分量的差异来评估误差来源。
• 评估位移误差如何传播至应力计算，尤其考虑滑动容差的影响。

实验结果

研究问题

• RQ1 有限元模型在模拟具有均匀滑动的正断层时，其位移场和应力场的重现精度如何？
• RQ2 数值解与解析解在位移、位移方向及应力分量上的误差大小是多少？
• RQ3 位移误差如何传播至计算得到的应力场？
• RQ4 迭代力校准方法是否能有效重现断层上期望的平滑滑动分布？
• RQ5 在该建模框架中，数值结果与解析结果之间差异的主要来源是什么？

主要发现

• 有限元模型生成的位移场与解析解相比，幅值差异小于5%。
• 位移方向误差低于15°，表明数值结果与解析结果在角度上具有良好的一致性。
• 应力分量差异小于35%，其中约一半误差可归因于数值设置中的滑动容差。
• 迭代力校准方案成功重现了平滑滑动分布，使直接有限元法可用于断层建模。
• 误差传播分析表明，位移不准确是导致应力场差异的主要原因，凸显了应力计算对位移误差的高度敏感性。
• 该模型与解析解表现出高度一致性，验证了其在正断层系统地质力学模拟中的适用性。