岩石介质中多裂纹扩展相互作用及其贯通机制的数值模拟.pdfVIP

维普资讯 第 2l卷 第 2期 地 震 V01．2l，No．2 2001年 4月 EARTHQUAKE A 口r．，2001 岩石介质中多裂纹扩展相互作用 及其贯通机制的数值模拟 唐春安 ，黄明利 ，张国民。，焦明若 (I_东北大学岩石破裂与失稳研究 中心，辽宁 沈阳 110006I9中国地震局分析预报中心，Jt~-- 10003B) 摘要 ：运用岩石破裂过程分析 RFPA卯系统t通过对岩石试样 中预置的一组右行右 阶雁列式裂 纹扩展过程的敷值模拟 ，研究了非均匀岩石舟质中多裂纹扩展的相互作用模式及其贯通机制。 数值模拟再现的受压混合型裂纹扩展过程中逐步演变的垒场变形过程 ，以及与细观 非均匀性有 关的声发射和断续扩展模式的 “岩桥 现象，清晰地揭示 出多裂纹扩展的相互作用及其贯通机 制。模拟结果说孵，对于岩石这种地质 舟质而言，忽略其 非均匀性影响，可能会掩盖岩 石变形与 破裂过程中的许多与非均匀性有关的特殊现象，包括声发射或微震模式、岩桥或雁行断裂等 。 关麓词 ：非均匀性；断裂力学 {裂纹相互作用及贯通 岩桥 ；敷值模拟 中圈分类号 ：P313；P584 文献标识码 ：A 文章编号 ：1000—3274(2001)02—0053-06 引言 裂纹扩展中的相互作用及其贯通机制是岩石断裂力学的重要研究课题，已有许多相关研 究的报道 。然而，尽管经典的断裂力学理论为裂纹的萌生、扩展问题的研究奠定了理论基 础，但迄今为止断裂力学理论对于处理岩石介质中多裂纹扩展的相互作用及贯通问题 ·却是极 其困难的。正如李世愚等 所指出的，以往的研究多半只考虑破裂的起始条件和起始的先后， 没有考虑到裂纹之间相互作用在破裂过程中的不断转化，即动态演变。此外 ，由于岩石是一种 典型的非均匀地质材料，它的性质很难真正满足断裂力学理论的基本假设，其中最重要的就是 均匀体和连续体假设。对岩石这种地质材料作为均匀性和连续性假设 ，尽管会带来数学处理上 的方便，但其结果则往往会掩盖岩石在加载过程 中所表现出的一些重要特征，这些特征小到岩 石试样中的声发射、岩桥行为，大到地壳介质中的微震、雁行断裂等。 本文试图运用数值模拟方法，对非均匀岩石介质中裂纹扩展的相互作用及其贯通机制作 一 初步的研究 数值模拟工具采用东北大学岩石破裂与失稳研究中心研究开发的岩石破裂过 - 收稿 日期 2000—08—02I惨改四日期 ：ZO00—1110 基金项 目：国家杰 出青年科学基金项 目(595254091I国隶重 点基础研究发展规划项 目( 作者简介 ：唐春安(1958一)，男．湖南黔阳凡，教授 ．1988年获博士学位 ·1993年博士生导师，主要从事岩石破裂过程 失稳研究。 维普资讯 地 震 21卷 程分 RFPA 系统 ”‘笔者 曾用这一方法研究了均匀介质中裂纹扩展的相互作用及其贯 通机制问题 ，取得了与实验一致的结果 。。本文在此基础上，模拟了单轴加载情况下非均匀岩 石试样 中多个预置裂纹的扩展过程，模拟过程再现了裂纹的扩展、相互作用及其贯通过程 。 1 数值模拟模型 l1 岩石破裂过程分析 l~pA1zD概避 所用的RFPAzD系统，是一个能模拟岩石介质逐渐破坏过程的数值模拟工具 。有关它的 基本思想已在多处介绍 ．本文不再赘述 。与其他 已有的逐渐破坏模型 (ProgressiveFail— ure／Fracture／DamageMode1)一样 ，RFPA 包括二个方面的功能：应力分析和破坏分析。 RFPA 的应力分析采用有限元法进行。破坏分析则是根据一定的破坏准则来检查材料中是 否有单元破坏 。对破坏单元则采用刚度特性退