断层断裂过程的数值模拟.pdfVIP

2011年第 9期 西部探矿工程 19 断层断裂过程的数值模拟 杨 杰 ，陈新 民 ，沈 建。 (1．中国公路工程咨询集团有限公司，北京 100097；2．南京工业大学地质工程研究所， 江苏 南京 210009；3．河海大学，江苏 南京 210098) 摘 要：逆断层的错动角度、错动的位移量以及土体性质等都会对上覆土层造成不同的影响。利用美 国Cornell大学断裂力学小组开发的Franc2d有限元软件，对断层地表露头位置、断层在土体中传播 模式以及断层的不同厚度与断层水平扩展距离间的关系，断层不同厚度与断层抬升的垂直高度间的 关系进行断层错动的模拟，得出一些有意义的结论。 关键词：断层 ；断裂；Franc2d数值模拟 中图分类号：P542．3 文献标识码 ：A 文章编号：1O04—5716(2O11)O9一OO19—06 我国是一个活断层分布比较广泛的国家[1屯]，也是 的模型中添加初始裂纹，软件就会 自动将裂纹尖端单元 一 个多震 国家。因此，我国在加紧大城市活断层探 奇异化，奇异化单元会随着裂纹尖端一起扩展。像裂缝 测 ¨3与地震危险性评价的同时，深入开展断层在土体 的传播这样，自动重新网格划分的计算法则消除了本属 中传播规律[5．6]的研究是必要的，弄清断层在地表的露 于裂纹尖端的网格，延伸裂缝，并且在新的裂缝尖端建 头位置对地下、地表结构物的抗震设计及对 已存在结构 起新的网格 。如图 1所示 。 物的影响分析都具有重要意义。很多学者在这方面也 2．2 数值模拟的分析过程 做了不少研究，利用模型试验[7．。]以及有限元模型，对断 一 个完整的FRANC2D模拟裂纹的扩展过程一般 层地表露头位置、断层在土体中传播模式以及断层到达 可以分为3个步骤：前处理、模拟计算和后处理。这 3 地表时断层位移与土体厚度的关系进行了模拟。 个步骤之间的联系与作用如图2所示。 最近几年来，像 FLAC3D、ANSYS、ABAQUS等 大型的有限元分析软件得到了广泛的应用，但是在使用 表 1 砂体厚度H=25cm时的试验结果 它们对疲劳裂纹问题进行计算时，往往需要进行 比较复 杂的操作，因此在实际应用上受到一定的限制。由美 国 Cornell大学断裂力学组开发的 FRANC2D(Fracture AnalysisCodein2Dimensions)L9- 有限元分析软件 ， 为解决断层的断裂模拟方面的问题提供了很好的平台。 该软件操作简单、兼容性好以及计算精度高等优点。 1 砂体模型试验结果 通过 45。断层错动 的模型箱试验，得出结果：如 表 2 砂体厚度为H=30cm的试验结果 表 1～3所示。 2 FRANC2D软件说明[9_埔] 2．1 FRANC2D模拟的原理 利用有限元软件模拟裂纹的扩展过程，基本上有以 下4种方法：(1)重分网格法；(2)增量法；(3)松解技术 法；(4)通过 特 殊 网格 移 动来 实 现裂 纹 扩 展 法 。 FRANC2D就是利用第 4种方法进行模拟，即在前处理 * 收稿 日期：2010~1卜11 修 回日期：2010—11—16 第一作者简介：杨杰(1981一)，男 (汉族)，江西上饶人，工程师，现从事岩土工程方面的设计研究工作。 20 西部探矿工程 2011年第 9期 表 3 砂体厚度为 H=35cm 的试验结果 比较典型的表达是：应力在任何线弹性裂纹体 内，都能