岩石细观统计损伤模型岩石破裂过程数值试验第一部分第一部分 RFPA系统概述系统概述第二部分第二部分 RFPA系统的基本原理系统的基本原理第三部分第三部分 RFPA-Basic可研究的问题可研究的问题第四部分第四部分 RFPA-Basic软件平台介绍软件平台介绍第五部分第五部分 工程算例模拟演示工程算例模拟演示主要内容一、一、RFPA系统概述系统概述RFPA系统概述系统概述 RFPA是是Realistic Failure Process Analysis的简称。的简称。它是一种基于有限元应力分析和统计损伤理论的材它是一种基于有限元应力分析和统计损伤理论的材料破裂过程分析数值计算方法，是一个能够模拟材料破裂过程分析数值计算方法，是一个能够模拟材料渐进破裂直至失稳全过程的数值试验工具。该方料渐进破裂直至失稳全过程的数值试验工具。该方法的一个重要特色是考虑了材料性质的非均匀性，法的一个重要特色是考虑了材料性质的非均匀性，是一种通过非均匀性模拟非线性、通过连续介质力是一种通过非均匀性模拟非线性、通过连续介质力学方法模拟非连续介质力学问题的材料破裂过程分学方法模拟非连续介质力学问题的材料破裂过程分析新型数值分析方法。析新型数值分析方法。 岩石力学问题，广义的讲包括岩石的破坏问题。岩石岩石力学问题，广义的讲包括岩石的破坏问题。岩石之所以产生非线性变形，是因为岩石在受载过程中其之所以产生非线性变形，是因为岩石在受载过程中其内部不断产生微细破裂的缘故。这种微细破裂的不断内部不断产生微细破裂的缘故。这种微细破裂的不断发展便导致最终的宏观破裂。通常的有限元方法尽管发展便导致最终的宏观破裂。通常的有限元方法尽管可以模拟演示的非线性变形，但只是在宏观上的一种可以模拟演示的非线性变形，但只是在宏观上的一种“形似形似”，而没有模拟出岩石在变形过程中微破裂过程，而没有模拟出岩石在变形过程中微破裂过程，因而不能做到因而不能做到“神似神似”。 为了解决岩石破裂过程的分析，采用有限元法、有限为了解决岩石破裂过程的分析，采用有限元法、有限差分法、边界元法、离散元法等数值模拟方法在全面差分法、边界元法、离散元法等数值模拟方法在全面解决复杂的岩土工程问题，例如岩石材料的非线性问解决复杂的岩土工程问题，例如岩石材料的非线性问题、岩体中节理、裂隙等不连续面对分析计算的影响题、岩体中节理、裂隙等不连续面对分析计算的影响等等方面不同程度的等等方面不同程度的存在缺陷。存在缺陷。 1995年软件系统创始年软件系统创始人唐春安教授针对这些问题提出了基于有限元基本理人唐春安教授针对这些问题提出了基于有限元基本理论，充分考虑岩石破裂过程中伴随的非线性、非均匀论，充分考虑岩石破裂过程中伴随的非线性、非均匀性和各向异性等特点的新的数值模拟方法性和各向异性等特点的新的数值模拟方法“RFPA方法方法”，即真实破裂过程分析方法。，即真实破裂过程分析方法。RFPA系统概述系统概述RFPA方法的主要要点方法的主要要点 将材料的不均质性参数引入到计算单元，宏观破坏是将材料的不均质性参数引入到计算单元，宏观破坏是单元破坏的积累过程。单元破坏的积累过程。 认为单元性质是弹认为单元性质是弹-脆性或弹脆性或弹-塑性的，单元的弹模和塑性的，单元的弹模和强度等其他参数服从某种分布，如正态、韦伯、均匀强度等其他参数服从某种分布，如正态、韦伯、均匀等。等。 认为当单元应力达到破坏的准则将发生破坏，并对破认为当单元应力达到破坏的准则将发生破坏，并对破坏单元进行刚度退化处理，故可以以连续介质力学方坏单元进行刚度退化处理，故可以以连续介质力学方法处理物理非线性介质问题。法处理物理非线性介质问题。 认为岩石的损伤量、声发射同破坏单元数成正比。认为岩石的损伤量、声发射同破坏单元数成正比。二、二、RFPA系统的基本原理系统的基本原理RFPA是一个以是一个以弹性力学性力学为应力分析工具、以力分析工具、以弹性性损伤理理论及其修正后的及其修正后的Coulomb破坏准破坏准则为介介质变形和破坏分析模形和破坏分析模块的岩石破裂的岩石破裂过程分析系程分析系统。其基本思路是：1.岩石介岩石介质模型离散化成由模型离散化成由细观基元基元组成的数成的数值模型，岩石模型，岩石介介质在在细观上是各向同性的上是各向同性的弹-脆性介脆性介质；2.假定离散化后的假定离散化后的细观基元的力学性基元的力学性质服从某种服从某种统计分布分布规律（本律（本书引入引入韦伯分布），由此建立伯分布），由此建立细观与宏与宏观介介质力学力学性能的性能的联系；系；3.按按弹性力学中的基元性力学中的基元线弹性性应力、力、应变求解方法，分析模求解方法，分析模型的型的应力、力、应变状状态。RFPA利用利用线弹性有限元方法作性有限元方法作为应力力计算器；算器；4.引入适当的基元破坏准引入适当的基元破坏准则（相（相变准准则）和）和损伤规律，基元律，基元相相变临界点用修正的界点用修正的CoulombCoulomb准准则和拉伸截断的和拉伸截断的库仑准准则；5.基元的力学性基元的力学性质随演化的随演化的发展是不可逆的；展是不可逆的；6.基元相基元相变前后均前后均为线弹性体；性体；7.岩石介岩石介质中的裂中的裂纹扩展是一个准静展是一个准静态过程，忽略因快速程，忽略因快速扩展引起的展引起的惯性力的影响。性力的影响。非均匀介质单元划分标准非均匀介质单元划分标准 从提高计算速度的角度讲：自然是单元取得大些比较从提高计算速度的角度讲：自然是单元取得大些比较合适，但较大的单元不仅不能反映细观非均匀的力学性合适，但较大的单元不仅不能反映细观非均匀的力学性质，使模型不能很好地反映工程实际，而且较大的单元质，使模型不能很好地反映工程实际，而且较大的单元尺寸本身也会给计算带来较大的误差。尺寸本身也会给计算带来较大的误差。 从计算精度的角度讲：在计算机速度允许的情况下，从计算精度的角度讲：在计算机速度允许的情况下，应尽可能的将单元划分得小些，这样不仅能够使模型更应尽可能的将单元划分得小些，这样不仅能够使模型更加真实地反映实际，而且也有利于提高计算的精度。加真实地反映实际，而且也有利于提高计算的精度。 单元大小的划分标准：根据精度，是否将单元划分的越单元大小的划分标准：根据精度，是否将单元划分的越小越好呢？未必！大家知道，在断裂力学中，为了数学小越好呢？未必！大家知道，在断裂力学中，为了数学处理的方便，裂纹的两端被假设为无穷小的尖端。然而，处理的方便，裂纹的两端被假设为无穷小的尖端。然而，按照这一假设，所计算出的裂纹尖端处的拉应力为无穷按照这一假设，所计算出的裂纹尖端处的拉应力为无穷大。这意味着什么？意味着只要给定任意有限值得加载，大。这意味着什么？意味着只要给定任意有限值得加载，裂纹便会迅速扩展，这显然是与实际情况不相符的。造裂纹便会迅速扩展，这显然是与实际情况不相符的。造成这一错误结论的原因是，在现实的介质特别是岩石介成这一错误结论的原因是，在现实的介质特别是岩石介质中，裂纹尖端实际上是有一定尺度的。这一尺度是与质中，裂纹尖端实际上是有一定尺度的。这一尺度是与介质的基本细观性质有关的，我们称之为介质的细观特介质的基本细观性质有关的，我们称之为介质的细观特征尺度。因此，征尺度。因此， 单元大小的划分标准为：只要数值模型单元大小的划分标准为：只要数值模型中的单元尺寸能够反映或者基本反映这种细观特征尺度，中的单元尺寸能够反映或者基本反映这种细观特征尺度，那么该模型的单元尺寸就是合理或者基本合理的。那么该模型的单元尺寸就是合理或者基本合理的。非均匀介质单元划分标准非均匀介质单元划分标准基元的引入基元的引入 细观力学认为；通过细观单元的变形、破坏的个体行为细观力学认为；通过细观单元的变形、破坏的个体行为的积累来反映宏观行为的演化，为研究介质变形和破裂的积累来反映宏观行为的演化，为研究介质变形和破裂的宏观行为提供了一种新的途径。的宏观行为提供了一种新的途径。所谓的基元，是构成所谓的基元，是构成介质的基本细观尺度单元，是在物理力学性质方面能够介质的基本细观尺度单元，是在物理力学性质方面能够代表介质特征的最小单元。代表介质特征的最小单元。基元是介质破裂研究最基本基元是介质破裂研究最基本的单位，破裂就是基元的破裂，比基元更小的破坏是不的单位，破裂就是基元的破裂，比基元更小的破坏是不存在的。在岩石破裂过程分析存在的。在岩石破裂过程分析RPFA系统中，为了能够系统中，为了能够充分考虑介质力学性能的非均匀性以及由这种非均匀性充分考虑介质力学性能的非均匀性以及由这种非均匀性引起的变形、破裂过程的复杂性，我们引入了三种特性引起的变形、破裂过程的复杂性，我们引入了三种特性的基元，即基质基元、空气基元和接触基元。的基元，即基质基元、空气基元和接触基元。基元的三种形态基元的三种形态 基质基元；是指基元在模型中的当前功能为实体介质。它基质基元；是指基元在模型中的当前功能为实体介质。它的性能由岩石的本构关系来描述。的性能由岩石的本构关系来描述。 空气基元；是指基元在模型中的当前功能为虚体特性。当空气基元；是指基元在模型中的当前功能为虚体特性。当单元介质在拉应力条件下发生断裂后，形成断裂面。就断单元介质在拉应力条件下发生断裂后，形成断裂面。就断裂面的物理本质而言，也就是应力的传递在此出现不连续裂面的物理本质而言，也就是应力的传递在此出现不连续或中断。通常的数值计算方法解决这一问题的方法是将单或中断。通常的数值计算方法解决这一问题的方法是将单元中的节点分开，或者是将单元从模型中去掉。但是，这元中的节点分开，或者是将单元从模型中去掉。但是，这样做的结果使得模型的数学处理变得极其复杂，而且一般样做的结果使得模型的数学处理变得极其复杂，而且一般不适合多裂纹、特别是多裂纹相互交叉的情形。不适合多裂纹、特别是多裂纹相互交叉的情形。RFPA系统系统采用的裂纹处理方法，即空气基元。当基元介质发生断裂采用的裂纹处理方法，即空气基元。当基元介质发生断裂后，我们不是将该单元从模型中去掉，而是用弹模极低的后，我们不是将该单元从模型中去掉，而是用弹模极低的基元性质取代原有的实体基元的性质，由于新的基元弹模基元性质取代原有的实体基元的性质，由于新的基元弹模极低，可以近似的认为实体介质的行为已不存在，这样在极低，可以近似的认为实体介质的行为已不存在，这样在不改变模型数学结构的前提下，却可以使得模型在总体特不改变模型数学结构的前提下，却可以使得模型在总体特性上能够反映出因基元破裂而引起的物理特性的改变。性上能够反映出因基元破裂而引起的物理特性的改变。 接触基元；压、剪破坏后的基元在一定范围内维持残余强度状接触基元；压、剪破坏后的基元在一定范围内维持残余强度状态。但是，对于现实中的介质来说，破坏后的介质在继续受压应态。但是，对于现实中的介质来说，破坏后的介质在继续受压应力、特别是各向均受压应力的条件下，将出现所谓的压密或压实力、特别是各向均受压应力的条件下，将出现所谓的压密或压实现象，其力学表现则是压密后的介质刚度不仅不降低，反而出现现象，其力学表现则是压密后的介质刚度不仅不降低，反而出现上升。对于已经形成的裂纹面而言（即空气基元），当裂隙两面上升。对于已经形成的裂纹面而言（即空气基元），当裂隙两面的介质在压应力作用下产生接触时，则应力仍可以通过接触面传的介质在压应力作用下产生接触时，则应力仍可以通过接触面传递。这时，则可以通过激活空气基元，使其刚度增加，起到传递递。这时，则可以通过激活空气基元，使其刚度增加，起到传递应力的作用。这就是所谓的接触基元特性。在上述两种情况下，应力的作用。这就是所谓的接触基元特性。在上述两种情况下，接触基元只能传递压应力，而不能传递拉应力。即当存在拉应力接触基元只能传递压应力，而不能传递拉应力。即当存在拉应力时，接触基元立刻转化成空气基元。时，接触基元立刻转化成空气基元。基元的三种形态基元的三种形态三种基元在不同条件下的转化三种基元在不同条件下的转化基质基元基质基元空气基元空气基元接触基元接触基元受拉应力破坏受拉应力破坏进行刚度退化处理进行刚度退化处理裂纹两面的介质受压作用产生接触，接触面传递压应力裂纹两面的介质受压作用产生接触，接触面传递压应力激活空气单元，单元刚度增加激活空气单元，单元刚度增加当压应力变为拉应力时，进行刚度退化处理当压应力变为拉应力时，进行刚度退化处理受压应力破坏受压应力破坏进行刚度强化处理进行刚度强化处理基元的相变基元的相变 所谓基元的相变，三种基元在一定条件下，将所谓基元的相变，三种基元在一定条件下，将由一种突然转化为另一种，这种其力学性质全由一种突然转化为另一种，这种其力学性质全然改观，叫做相变。其临界条件即为相变点。然改观，叫做相变。其临界条件即为相变点。分离相残余相弹性相残余相接触相基元及其特性基元及其特性分离点相变点2相变点1接触点图中个参数的意义图中个参数的意义:1类相变阀值类相变阀值1类相变残余阀值类相变残余阀值2类相变阀值类相变阀值2类相变残余阀值类相变残余阀值最大压应变最大压应变极限压应变极限压应变最大拉应变最大拉应变极限拉应变极限拉应变基元相变之前，它代表实际的完好介质，基元相变之前，它代表实际的完好介质，具有线弹性应力具有线弹性应力-应变特性。应变特性。基元相变之后的一定范围内，它代表相变基元相变之后的一定范围内，它代表相变后阶段的弱介质，具有残余阀值特性。后阶段的弱介质，具有残余阀值特性。如果相变后的基元被高度挤压，则基元犹如果相变后的基元被高度挤压，则基元犹如被压密的介质。如果相变后的基元受拉，如被压密的介质。如果相变后的基元受拉，使其拉应变达到一定的值时不再将基元看成使其拉应变达到一定的值时不再将基元看成是它所代表的介质，处理成空气基元。是它所代表的介质，处理成空气基元。RFPA分析过程流程图分析过程流程图开始开始实体建模和网格划分，用统计分布函数，实体建模和网格划分，用统计分布函数，赋每一个基元的刚度和相变值等赋每一个基元的刚度和相变值等施加荷载产生一个新的边界位移或载荷施加荷载产生一个新的边界位移或载荷形成新的刚度矩阵形成新的刚度矩阵将相变基元进将相变基元进行弱化处理行弱化处理线弹性有限线弹性有限元求解器元求解器计算基元节点力和位移计算基元节点力和位移根据相变准则判断根据相变准则判断是否有基元相变是否有基元相变加载是否加载是否需要结束需要结束结束结束是是否否是是否否实体建模实体建模应力分析应力分析相变分析相变分析三、三、RFPA-Basic可研究的问题可研究的问题RFPA-2D基本版基本版 RFPA-2D基本版：集中了基本版：集中了RFPA系列软系列软件的基本功能，可进行岩石、混凝土等脆件的基本功能，可进行岩石、混凝土等脆性材料受载的变形破坏分析，其声发射模性材料受载的变形破坏分析，其声发射模式的探测，亦可进行简单的工程应用分析，式的探测，亦可进行简单的工程应用分析，如地下工程开挖及支护过程中的应力场、如地下工程开挖及支护过程中的应力场、位移场监测及声发射（微震）监控等。位移场监测及声发射（微震）监控等。 请看岩石材料单轴压缩的破坏过程请看岩石材料单轴压缩的破坏过程岩石材料单轴压缩的破坏过程及声发射岩石材料单轴压缩的破坏过程及声发射四、四、RFPA系统平台介绍系统平台介绍RFPA系统平台介绍系统平台介绍、RFPA系统系统计算方面的特点：计算方面的特点： 允许模拟计算由于分步开挖引起的应力重新分布对进一步变形和破坏允许模拟计算由于分步开挖引起的应力重新分布对进一步变形和破坏过程的影响，也就是说在模拟地下开采中伴随的破坏过程时更接近于过程的影响，也就是说在模拟地下开采中伴随的破坏过程时更接近于实际情况。实际情况。 可以通过专门提供的作图工具和参数输入模块，在模型中考虑模拟材可以通过专门提供的作图工具和参数输入模块，在模型中考虑模拟材料的微观料的微观缺陷，也可以考虑节理、裂隙等宏观缺陷。缺陷，也可以考虑节理、裂隙等宏观缺陷。 可以模拟自重引起的破坏过程，新的可以模拟自重引起的破坏过程，新的RFPA2D软件系统增加了对地下软件系统增加了对地下工程开挖，地下破坏、地表沉陷、采动影响下煤岩顶板冒落、边坡失工程开挖，地下破坏、地表沉陷、采动影响下煤岩顶板冒落、边坡失稳等问题的模拟功能。稳等问题的模拟功能。 增加了对多种统计分布函数如韦伯分布、正态分布、均匀分布等在材增加了对多种统计分布函数如韦伯分布、正态分布、均匀分布等在材料特性中的嵌入，来考虑材料力学参数（强度、弹模等）的非均匀性料特性中的嵌入，来考虑材料力学参数（强度、弹模等）的非均匀性分布特征，从而可以从本质上研究岩石变形的非线性特征。分布特征，从而可以从本质上研究岩石变形的非线性特征。 新的新的RFPA2D软件系统增加了对流固耦合（如：水力压裂、底板突水、软件系统增加了对流固耦合（如：水力压裂、底板突水、水工中岩石、体渗流）、气固耦合（媒岩体瓦斯突出）、温度应力场水工中岩石、体渗流）、气固耦合（媒岩体瓦斯突出）、温度应力场耦合问题的模拟分析功能。耦合问题的模拟分析功能。 边坡稳定性分析模块，主要为强度折减法和离心机法。边坡稳定性分析模块，主要为强度折减法和离心机法。、 RFPA系统使用方面的特点：系统使用方面的特点： RFPA2D软件系统由于从设计到具体的实现过程就时刻注意到系统的可视化，所以新系软件系统由于从设计到具体的实现过程就时刻注意到系统的可视化，所以新系统的用户界面不出现任何与数值计算方法有关的术语，用户无需掌握专门的数值计算方统的用户界面不出现任何与数值计算方法有关的术语，用户无需掌握专门的数值计算方法，只需提供与研究对象有关的几何参数和力学参数，就可以用几何作图方式构造力学法，只需提供与研究对象有关的几何参数和力学参数，就可以用几何作图方式构造力学模型；所有操作都是针对图形进行的。用户只需利用内藏的高性能作图工具绘出结构物模型；所有操作都是针对图形进行的。用户只需利用内藏的高性能作图工具绘出结构物的图形，并为组成结构物的各个部分指定材料的力学参数，然后指定边界条件，即可进的图形，并为组成结构物的各个部分指定材料的力学参数，然后指定边界条件，即可进行计算分析。行计算分析。 赋值（前处理）。赋值（前处理）。RFPA对任何一个实际模拟对象的赋值，包括网格大小、力学参数、对任何一个实际模拟对象的赋值，包括网格大小、力学参数、边界条件、加载条件等，都是基于面向用户的对话框。边界条件、加载条件等，都是基于面向用户的对话框。 计算。计算。RFPA可快速进行可快速进行100,000个单元以上的模型分析计算，而且最大容量只取决于计个单元以上的模型分析计算，而且最大容量只取决于计算机的硬件性能。在算机的硬件性能。在RFPA系统中，应力分析和相变分析（破坏分析）是相互独立的，有系统中，应力分析和相变分析（破坏分析）是相互独立的，有限元仅完成应力、应变计算，不参与相变分析。但整个应力分析、相变分析过程是连贯限元仅完成应力、应变计算，不参与相变分析。但整个应力分析、相变分析过程是连贯依此完成的，不需要用户进行任何干预。依此完成的，不需要用户进行任何干预。 显示（后处理）。显示（后处理）。在在RFPA中专门设置了一组显示工具条和菜单栏用于对模拟结果进行中专门设置了一组显示工具条和菜单栏用于对模拟结果进行显示操作。显示操作。RFPA实时显示中间的计算结果，从而可以及时反馈给用户计算进度。实时显示中间的计算结果，从而可以及时反馈给用户计算进度。RFPA还具有对计算结果重画图形的功能，以便获得更为理想的后处理结果。还具有对计算结果重画图形的功能，以便获得更为理想的后处理结果。 图形编辑。图形编辑。RFPA具有对模拟结果图形进行缩放、拷贝等编辑功能，还可以直接运用具有对模拟结果图形进行缩放、拷贝等编辑功能，还可以直接运用windows的剪切、复制、粘贴等功能将模拟结果图形传输到图形处理软件（如的剪切、复制、粘贴等功能将模拟结果图形传输到图形处理软件（如MS-Word，CorelDraw等）进行后处理编辑。等）进行后处理编辑。RFPA-2D工作平台的装入工作平台的装入：通过鼠标双击桌面上：通过鼠标双击桌面上RFPA图标装入图标装入：通过开始菜单装入，开始：通过开始菜单装入，开始程序程序RFPA：通过：通过Windows Explorer中运行中运行RFPAstudio.exe 命令：命令：c：RFPARFPAstudio.exe初始工作平台初始工作平台新建工程文档新建工程文档打开已有的打开已有的*rfp类型文档类型文档设计平台设计平台RFPA 模型模型工作区工作区显示工具条显示工具条基本工具条基本工具条作图工具条作图工具条材料特性类型材料特性类型材料特性分布图材料特性分布图菜单条菜单条调色板调色板状态条状态条工具条介绍工具条介绍 (tool bar)工具条：是为了方便菜单中的主要命令的使用而设置的快捷工具工具条：是为了方便菜单中的主要命令的使用而设置的快捷工具基基基基本本本本工工工工具具具具条条条条New新建新建Open打开打开Save保存保存Print打印打印Copy复制复制Zoom缩放缩放Full Screen全屏显示全屏显示Mesh dividing网格划分网格划分Boundary conditions边界条件边界条件Control conditions控制条件控制条件Run step by step单步运行单步运行Auto run 自动连续运行自动连续运行Stop run强制停止运行强制停止运行作作作作图图图图工工工工具具具具条条条条Draw by mouse鼠标键盘输入切换鼠标键盘输入切换Material property材料类型选择材料类型选择Select area区域选择区域选择Rectangle矩形矩形Circle圆圆Line单直线单直线Double line双直线双直线Polyline多线多线Polygon任意多边形任意多边形Compositematerial复合材料复合材料Closed Arc封闭弧封闭弧Arc圆弧圆弧Anchor支护支护复合材料增强复合材料增强介质颗粒设置介质颗粒设置显示工具条显示工具条显示工具条显示工具条Shear stress剪应力图剪应力图Max. principal stress最大主应力图最大主应力图Min. principal stress最小主应力图最小主应力图Elastic弹性模量图弹性模量图Acoustic Emission声发射图声发射图e/AE弹模和声发射图弹模和声发射图Load-step curve载荷与加载步曲线载荷与加载步曲线AE curve声发射曲线声发射曲线Area AE圈定声发射图圈定声发射图Multi-Element Information多单元信息多单元信息Slide play幻灯播放幻灯播放Displacement vector位移矢量图位移矢量图Stress vector应力矢量图应力矢量图颜色灰度调整颜色灰度调整颜色灰度调整颜色灰度调整(Color Adjustment):用户根据自己的研究和兴趣需要，可进行结果图颜色灰度的调整；用户根据自己的研究和兴趣需要，可进行结果图颜色灰度的调整；单击左边颜色框，调节右边单击左边颜色框，调节右边Minimum 和和Maximum滑动按钮。滑动按钮。网格划分网格划分 Y Length(mm),X length(m): Y方向和方向和X方向的尺寸均为所研究问题的实际尺寸，单位为毫米，建立模型的时方向的尺寸均为所研究问题的实际尺寸，单位为毫米，建立模型的时候需要将实际尺寸换算成毫米，单元尺寸为实际候需要将实际尺寸换算成毫米，单元尺寸为实际 尺寸除单元个数。尺寸除单元个数。Rows(Elements),Cols(Elements): 单元列数和行数，单元尺寸的划分原则前面已讲述。单元列数和行数，单元尺寸的划分原则前面已讲述。Heterogeneity Index行：行：为均质度，为均质度，weibull统计分布函数中的参数统计分布函数中的参数m，反映岩石介质的均质性。若模拟地下，反映岩石介质的均质性。若模拟地下工程等必须考虑模型本身的重力的工程问题时，需要输入自重的参数，注意自重的单位：工程等必须考虑模型本身的重力的工程问题时，需要输入自重的参数，注意自重的单位：N/mm3。Mean Value: 为单元物理力学参数的平均值。为单元物理力学参数的平均值。Weibull分布函数：分布函数：强度特性设置强度特性设置摩擦角摩擦角(Friction)：材料在垂直力作用下发生剪切破坏时错动面的倾角材料在垂直力作用下发生剪切破坏时错动面的倾角 （材料固有）。（材料固有）。压拉比压拉比(C/T ratio)：材料抗压强度与抗拉强度的比值。材料抗压强度与抗拉强度的比值。残余强度百分比残余强度百分比(Residual strength(%)：破坏后的残余强度与破坏前强度的比值。破坏后的残余强度与破坏前强度的比值。残余泊松比百分比残余泊松比百分比(Residual poisson(%)：破坏后的泊松比与破坏前泊松比的比值。破坏后的泊松比与破坏前泊松比的比值。最大应变最大应变(Maximum strain)：极限最大应变与弹性最大应变的比值（理解参见基元与其特性）。极限最大应变与弹性最大应变的比值（理解参见基元与其特性）。力边界条件设置力边界条件设置加载方式加载方式(Load mode)：位移加载、应力加载、围压位移加载、应力加载、围压控制、自由。控制、自由。加载初始步加载初始步(Initial step)：开始施加力边界条件的开始开始施加力边界条件的开始步。步。加载初始值加载初始值(Initial value)：计算循环中加载起始步的计算循环中加载起始步的值（选择不同的加载方式则对应不同的加载值）。值（选择不同的加载方式则对应不同的加载值）。加载增量加载增量(Increment)：两个加载步之间变化量。正值两个加载步之间变化量。正值表示压，负值表示拉。表示压，负值表示拉。加载终值加载终值(Final value)：用于设置加载循环中某方向用于设置加载循环中某方向上的终值。上的终值。加载终止步加载终止步(Final step)：设置加载的终止步。设置加载的终止步。控制条件设置控制条件设置执行完网格划分命令就要对模拟区进行控执行完网格划分命令就要对模拟区进行控制条件设置。制条件设置。总控制步总控制步(Total step)：控制整个程序的循控制整个程序的循环计算步。也就是所计算的循环步最终取环计算步。也就是所计算的循环步最终取决于此处的步数。决于此处的步数。步中步步中步(step In step)：根据计算精度而显根据计算精度而显示的两步之间的小步。示的两步之间的小步。求解类型求解类型(Load type)：可进行平面应变问可进行平面应变问题、平面应力问题和轴对成问题的求解。题、平面应力问题和轴对成问题的求解。重力方向重力方向(Weight direction)：重力方向的重力方向的确定。确定。系统自定义系统自定义 (self define)系统自定义系统自定义(self define)图片类型和尺寸图片类型和尺寸(Picture Type Size)结果图类型结果图类型(Type)：此命令主要是定义此命令主要是定义在模型计算完成后要显示图形的类型和在模型计算完成后要显示图形的类型和尺寸。根据研究问题的需要将相应的结尺寸。根据研究问题的需要将相应的结果图选项上钩上即可。果图选项上钩上即可。图形尺寸图形尺寸(Size)：对系统生成的结果图对系统生成的结果图进行大小设置，选进行大小设置，选Default默认，系统默认，系统根据窗口大小随机生成。根据窗口大小随机生成。矢量单元间隔矢量单元间隔(Vector element intervals)：因为要表示矢量场，用箭因为要表示矢量场，用箭头表示矢量的大小方向，箭头方向为矢头表示矢量的大小方向，箭头方向为矢量方向，箭头长短表示矢量大小。量方向，箭头长短表示矢量大小。位移矢量放大系数位移矢量放大系数 (Displacement vector amplified coeff)：为了计算结为了计算结果显示的我们能更好的看清以说明问题，果显示的我们能更好的看清以说明问题，对原始矢量大小给予适当的放大。对原始矢量大小给予适当的放大。应力矢量放大系数应力矢量放大系数(Stress vector amplified coeff)：解释同位移放大系解释同位移放大系数。数。注：这些参数定义可在计算运行命令之前执行，也可注：这些参数定义可在计算运行命令之前执行，也可在计算命令之后执行。在计算命令之后执行。显示步显示步(Show step)显示步显示步(Show step)：此命令在设置重此命令在设置重画计算结果图时候设置，设置重画计画计算结果图时候设置，设置重画计算步的间隔。也可在算步的间隔。也可在input redraw step 下面直接输入要重画的计算步，用英下面直接输入要重画的计算步，用英文状态下的逗号相隔。文状态下的逗号相隔。重画重画(Redraw)：是系统对计算数据进是系统对计算数据进行重新生成图像的操作，原计算数据行重新生成图像的操作，原计算数据不变化，只是新生成的结果图将覆盖不变化，只是新生成的结果图将覆盖和上次生成的名称相同的图形。不同和上次生成的名称相同的图形。不同于系统重置于系统重置(reset)。系统自定义系统自定义(self define) 声发射图设置声发射图设置(AE Picture Setting)声发射圆大小声发射圆大小(Type size)：设置声设置声发射圆大小（代表相对能量或震发射圆大小（代表相对能量或震级大小），级大小），Amplified radius是指是指在原能量圆大小的基础上放大的在原能量圆大小的基础上放大的倍数。倍数。所有步所有步(For all steps)：显示当前显示当前所有步的声发射圆。所有步的声发射圆。当前步当前步(For current step)：显示当显示当前单步的声发射圆。前单步的声发射圆。发射圆颜色发射圆颜色(AE color)：用声发射用声发射圆颜色进行区别单元的破坏是抗圆颜色进行区别单元的破坏是抗拉破坏、抗压破坏、还是抗剪破拉破坏、抗压破坏、还是抗剪破坏。双击颜色框进行调节。坏。双击颜色框进行调节。系统自定义系统自定义(self define)高级设置高级设置(Advanced Setting)主要对破坏后和空气单元的参数进行调整主要对破坏后和空气单元的参数进行调整(初初学者将词项默认学者将词项默认)。分离系数分离系数(Detached coeff)：主要为了处理破主要为了处理破坏后单元的显示状态，当破坏单元的变形坏后单元的显示状态，当破坏单元的变形(受受拉拉)达到分离系数指定的变形量时，认为单元达到分离系数指定的变形量时，认为单元已彻底分离，将不在显示此单元，值为已彻底分离，将不在显示此单元，值为1-100。设置为。设置为5时，表示当单元的最大变形超过时，表示当单元的最大变形超过5倍单元原始尺寸时，不再显示此单元。倍单元原始尺寸时，不再显示此单元。灰度系数灰度系数(Gray degree coeff)：是通过去掉高：是通过去掉高亮度单元来调整总体的亮度范围，值在亮度单元来调整总体的亮度范围，值在01之间，其意义为：假设模型单元数为之间，其意义为：假设模型单元数为10，000个，灰度系数为个，灰度系数为0.01，则表示在显示时去掉，则表示在显示时去掉10，0000.01=100个最亮的个单元级别。个最亮的个单元级别。大位移大位移(Large displacement)：是在处理自重是在处理自重问题时，当破坏的单元下落的位移超过设置问题时，当破坏的单元下落的位移超过设置值时，则认为破坏单元脱离母体单元，值为值时，则认为破坏单元脱离母体单元，值为单元边长的倍数。单元边长的倍数。系统自定义系统自定义(self define)变形放大系数变形放大系数(Deformation)：对破坏模对破坏模型的尺寸在显示时做放大处理，如有时型的尺寸在显示时做放大处理，如有时计算完成后发现模型没有变化则可能是计算完成后发现模型没有变化则可能是放大系数设置的太小；反之，当计算完放大系数设置的太小；反之，当计算完成后显示的图形出现全黑时，则有可能成后显示的图形出现全黑时，则有可能是因为放大系数设置过大的缘故。是因为放大系数设置过大的缘故。计算精度系数计算精度系数(Accurate of stress)：表示表示控制计算精度，如设置为控制计算精度，如设置为0.0001时，则时，则表示当破坏单元达到总单元数万分之一表示当破坏单元达到总单元数万分之一时，系统不再做当前不的循环计算。时，系统不再做当前不的循环计算。系统自定义系统自定义(self define)其他设置其他设置(Others setting)五、工程算例五、工程算例模拟演模拟演示示岩石的压缩加载数值试验岩石的压缩加载数值试验岩石的剪切加载试验岩石的剪切加载试验 RFPA工程算例工程算例（基本模块）（基本模块）RFPA计算步骤计算步骤（三大步）（三大步）构造模型构造模型 （前处理）（前处理）1建立工程文档2网格划分和参数赋值3插入模型结构4边界条件和控制条件计算计算单步运行连续运行结果分析结果分析（后处理）（后处理）岩石的压缩加载数值试验岩石的压缩加载数值试验试验内容：试验内容： 岩石压缩加载的应力岩石压缩加载的应力-应变曲线应变曲线 破坏模式破坏模式 声发射声发射第第1步：模型建立及参数选择步：模型建立及参数选择 模型：模型： 试样模型尺寸试样模型尺寸150mm100mm ，网个划分为，网个划分为150100个基元个基元 。采用平面应力问。采用平面应力问题，整个加载题，整个加载 过程通过位移加载方式。加载位移量过程通过位移加载方式。加载位移量s=0.002mm.力学性质参数力学性质参数如下表：如下表： 表 1 基元的力学参数、相变准则均质度弹性模量强度Poisson自重1.5600002000.250摩擦角C/T 比Max strain TMax strain CStrength criterion30101.5200M-C CriterionRFPA数数值模型值模型（上下端面与加载板连接）（上下端面与加载板连接）自由边自由边自由边自由边第第2步：计算步：计算几分钟后几分钟后第第3步：计算结果分析步：计算结果分析点击一次观看破裂过程点击一次观看破裂过程从数值试验得到的载荷从数值试验得到的载荷-位移全过程曲线再现了如下基本的岩石力学性质：位移全过程曲线再现了如下基本的岩石力学性质：（1）线性变形阶段。在加载的初期，载荷）线性变形阶段。在加载的初期，载荷-位移曲线几乎是线性的。位移曲线几乎是线性的。（2）非线性变形阶段。当载荷达到试件最大承载能力的）非线性变形阶段。当载荷达到试件最大承载能力的50%左右时，试件的左右时，试件的变形开始偏离线性，部分基元破坏。变形开始偏离线性，部分基元破坏。（3）软化阶段。当达到最大载荷之后，使试件进一步变形的载荷越来越小，）软化阶段。当达到最大载荷之后，使试件进一步变形的载荷越来越小，进入弱化阶段，直至试件产生宏观破坏。进入弱化阶段，直至试件产生宏观破坏。曲线分析曲线分析N-S曲线曲线岩石的剪切加载试验岩石的剪切加载试验试验内容：试验内容： 岩石压缩加载的应力岩石压缩加载的应力-应变曲线应变曲线 破坏模式破坏模式 声发射声发射第第1步：模型建立及参数选择步：模型建立及参数选择 模型：模型： 试样模型尺寸试样模型尺寸100mm100mm ，网个划分为，网个划分为100100个基元个基元 。采用平面应力问题，整个加载过程通过位移加载方式。加载。采用平面应力问题，整个加载过程通过位移加载方式。加载位移量位移量s=0.005mm.力学性质参数如下表：力学性质参数如下表： 均质度弹性模量强度Poisson自重3500002000.250摩擦角C/T 比Max strain TMax strain CStrength criterion30101.5200M-C Criterion表 1 基元的力学参数、相变准则RFPA数值模型数值模型注：预先切出注：预先切出25mm的水平引导裂纹，垂直方向加围压，每步的水平引导裂纹，垂直方向加围压，每步0.5Mpa，终值加到第五步，达到最大值，终值加到第五步，达到最大值2.5Mpa。第第2步：计算步：计算几分钟后几分钟后第第3步：计算结果分析步：计算结果分析剪切破裂过程剪切破裂过程曲线分析曲线分析N-S曲线曲线AE能量图能量图