RFPA软件的基本原理说课讲解.doc

RFPA软件的基本原理 精品文档 精品文档 收集于网络，如有侵权请联系管理员删除 收集于网络，如有侵权请联系管理员删除 精品文档 收集于网络，如有侵权请联系管理员删除 RFPA软件的基本原理 1 RFPA系统简介 本文所运用的RFPA系统，即岩石破裂过程分析系统（Rock Failure Process Analysis System），是基于弹性损伤模型的一个工具，1991年唐春安教授提出了岩石细观单元强度满足某个正态统计分布的假设，认为细观非均匀性是造成准脆性材料宏观非线性的根本原因，用统计损伤的本构关系考虑了岩石材料的非均匀性和缺陷分布的随机性。此后，又把这种材料性质的统计分布假设结合到数值计算方法(如有限元法)中，并对满足给定强度准则的单元进行破坏处理，实现了对岩石破坏过程中的起裂、变形局部化、成核，最终导致失稳破坏的全程模拟，使有限元技术的应用发展到直接模拟岩体破裂过程的新阶段。 在非均匀介质中单元划分引入了基元这一标准，我们将反映介质细观尺度的基本单元称为基元。基元是构成介质的基本细观尺度单元，是在物理力学性质方面能够代表介质特征的最小单元。基元是构成介质的基本细观尺度单元，是在物理力学性质方面能够代表介质特征的最小单元。基元是介质破裂的研究最基本的单位，破裂就是基元的破裂。 1.1 RFPA基本原理 RFPA是一个以弹性力学为应力分析工具、以弹性损伤理论及其修正后的Coulomb破坏准则为介质变形和破坏分析模块的岩石破裂过程分析系统。其基本思路是： （1）岩石介质模型离散化成有细观基元组成的数值模型，岩石介质在细观上市各向同性的弹-脆性介质。 （2）假定离散化后的细观基元的力学性质服从某种统计分布规律（Weibull分布），由此建立细观与宏观介质力学性能的联系。 （3）按弹性力学中的基元线弹性应力、应变求解方法，分析模型的应力、应变状态。RFPA利用线弹性有限元方法作为应变计算器。 （4）引入适当的基元破坏准则（相变准则）和损伤规律，基元的相变临界点用修正的Coulomb准则。 （5）基元的力学性质随演化的发展是不可逆的。 （6）基元相变前后均为线弹性体。 （7）岩石介质中的裂纹扩展是一个准静态过程，忽略因快速扩展引起的惯性力的影响。 1.2 网格剖分 RFPA选取等面积四截点的四边形单元剖分计算机。为了使问题的解答足够精确，RFPA方法要求模型中的单元能足够小（相对于宏观介质），以能足够精确地反映介质的非均匀性。但它又必须是足够大（包含一定数量的矿物和胶结物颗粒，以及微裂隙、空洞等细小缺陷），因为作为子系统的单元实际上仍是一个自由度足够大的系统，它具有远大于微观尺度的细观尺度，。这一要求正是为了保证使剖分后的单元性质尽量接近基元的性质。尽管这样会增加计算量，但是问题的处理变得简单，而且随着计算机技术的高速发展，计算能力瓶颈的影响将会被逐渐消除。 由于模型中的基元数量足够多，宏观的力学行为，本质上是介质大量基元力学行为的集体效应。但是每个基元的个体行为对宏观性能的影响却是有限的。正如夏梦芬等说指出的：“对单个个体的力学性能做出详尽无遗的描述不仅不可能，而且也不必要，只需给出一种详略得当的描写即可。这个描述即包含与问题有关的最重要信息，又不因为过于繁杂而失去可操作性。” RFPA系统正是基于这种原则对基元的力学行为进行描述的。 1.3 基元赋值 本文采用Monte-Carlo方法和统计描述相结合对基元进行初始化赋值。设模型中所有基元的弹性模量平均值为，代表了具有某弹性模量基元的分布值弹性模量Weibull分布函数的积分为： (4-1） 式中，为具有弹性模量的基元统计数量，由式（4-2）统计分布构成的基元组成一个样本空间，在均值不变的情况下，由于值的差别，积分空间分布却不完全一样。这些基元构成的岩石类介质的细观平均性质可能大体一致（相同），但是由于细观结构的无序性，使得基元的空间排列方式有显著的不同。这种细观上的无序性正好体现了岩石类介质独特的离散性特征。 1.4 应力计算 在RFPA系统中，整个分析对象被离散成若干具有不同物理力学性质的基元。为了求解各个基元的应力、应变状态，各基元之间需要满足力的平衡、变形协调方程和一定的应力、应变关系（物理方程）。在众多有关应力、应变的数值计算方法中，有限元是最理想的一种数值计算方法之一。它是将一个连续的介质离散成由诸多有限大小的单元组成的结构物体，然后通过力的平衡方程、几何方程和物理方程求解各个离散体的力学状态。因此，在RFPA系统中利用有限元作为应力分析求解器，它完成外载荷作用下对象内部各基元的应力、应变状态的计算工作。在RFPA系统中，应力分析求解器和相变分析相互独立，应力计算器仅完成应力、应变计算，不参与相变分析。 1.5 RFPA分