rfpa系统介绍

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,岩石细观统计损伤模型,岩石破裂过程数值试验,主讲人：谭志宏,第一部分,RFPA,系统概述,第二部分,RFPA,系统的基本原理,第三部分,RFPA-Basic,可研究的问题,第四部分,RFPA-Basic,软件平台介绍,第五部分 工程算例模拟演示,主要内容,一、,RFPA,系统概述,RFPA,系统概述,RFPA,是,Realistic Failure Process Analysis,的简称。它是一种基于有限元应力分析和统计损伤理论的材料破裂过程分析数值计算方法，是一个能够模拟材料渐进破裂直至失稳全过程的数值试验工具。该方法的一个重要特色是考虑了材料性质的非均匀性，是一种通过非均匀性模拟非线性、通过连续介质力学方法模拟非连续介质力学问题的材料破裂过程分析新型数值分析方法。,岩石力学问题，广义的讲包括岩石的破坏问题。岩石之所以产生非线性变形，是因为岩石在受载过程中其内部不断产生微细破裂的缘故。这种微细破裂的不断发展便导致最终的宏观破裂。通常的有限元方法尽管可以模拟演示的非线性变形，但只是在宏观上的一种“形似”，而没有模拟出岩石在变形过程中微破裂过程，因而不能做到“神似”。,为了解决岩石破裂过程的分析，采用有限元法、有限差分法、边界元法、离散元法等数值模拟方法在全面解决复杂的岩土工程问题，例如岩石材料的非线性问题、岩体中节理、裂隙等不连续面对分析计算的影响等等方面不同程度的,存在缺陷。,1995,年软件系统创始人唐春安教授针对这些问题提出了基于有限元基本理论，充分考虑岩石破裂过程中伴随的非线性、非均匀性和各向异性等特点的新的数值模拟方法,“,RFPA,方法”,，即真实破裂过程分析方法。,RFPA,系统概述,RFPA,方法的主要要点,将材料的不均质性参数引入到计算单元，宏观破坏是单元破坏的积累过程。,认为单元性质是弹,-,脆性或弹,-,塑性的，单元的弹模和强度等其他参数服从某种分布，如正态、韦伯、均匀等。,认为当单元应力达到破坏的准则将发生破坏，并对破坏单元进行刚度退化处理，故可以以连续介质力学方法处理物理非线性介质问题。,认为岩石的损伤量、声发射同破坏单元数成正比。,二、,RFPA,系统的基本原理,RFPA,是一个以弹性力学为应力分析工具、以弹性损伤理论及其修正后的,Coulomb,破坏准则为介质变形和破坏分析模块的岩石破裂过程分析系统。,其基本思路是：,岩石介质模型离散化成由细观基元组成的数值模型，岩石介质在细观上是各向同性的弹,-,脆性介质；,假定离散化后的细观基元的力学性质服从某种统计分布规律（本书引入韦伯分布），由此建立细观与宏观介质力学性能的联系；,按弹性力学中的基元线弹性应力、应变求解方法，分析模型的应力、应变状态。,RFPA,利用线弹性有限元方法作为应力计算器；,引入适当的基元破坏准则（相变准则）和损伤规律，基元相变临界点用修正的,Coulomb,准,则和拉伸截断的库仑准则；,基元的力学性质随演化的发展是不可逆的；,基元相变前后均为线弹性体；,岩石介质中的裂纹扩展是一个准静态过程，忽略因快速扩展引起的惯性力的影响。,非均匀介质单元划分标准,从提高计算速度的角度讲：自然是单元取得大些比较合适，但较大的单元不仅不能反映细观非均匀的力学性质，使模型不能很好地反映工程实际，而且较大的单元尺寸本身也会给计算带来较大的误差。,从计算精度的角度讲：在计算机速度允许的情况下，应尽可能的将单元划分得小些，这样不仅能够使模型更加真实地反映实际，而且也有利于提高计算的精度。,单元大小的划分标准：根据精度，是否将单元划分的越小越好呢？未必！大家知道，在断裂力学中，为了数学处理的方便，裂纹的两端被假设为无穷小的尖端。然而，按照这一假设，所计算出的裂纹尖端处的拉应力为无穷大。这意味着什么？意味着只要给定任意有限值得加载，裂纹便会迅速扩展，这显然是与实际情况不相符的。造成这一错误结论的原因是，在现实的介质特别是岩石介质中，裂纹尖端实际上是有一定尺度的。这一尺度是与介质的基本细观性质有关的，我们称之为介质的细观特征尺度。因此， 单元大小的划分标准为：只要数值模型中的单元尺寸能够反映或者基本反映这种细观特征尺度，那么该模型的单元尺寸就是合理或者基本合理的。,非均匀介质单元划分标准,基元的引入,细观力学认为；通过细观单元的变形、破坏的个体行为的积累来反映宏观行为的演化，为研究介质变形和破裂的宏观行为提供了一种新的途径。,所谓的基元，是构成介质的基本细观尺度单元，是在物理力学性质方面能够代表介质特征的最小单元。,基元是介质破裂研究最基本的单位，破裂就是基元的破裂，比基元更小的破坏是不存在的。在岩石破裂过程分析,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/mm,3,。,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.5,60000,200,0.25,0,摩擦角,C/T,比,Max strain T,Max strain C,Strength criterion,30,10,1.5,200,M-C Criterion,RFPA,数,值模型,（上下端面与加载板连接）,自由边,自由边,第,2,步：计算,几分钟后,第,3,步：计算结果分析,点击一次观看破裂过程,从数值试验得到的载荷,-,位移全过程曲线再现了如下基本的岩石力学性质：,（,1,）线性变形阶段。在加载的初期，载荷,-,位移曲线几乎是线性的。,（,2,）非线性变形阶段。当载荷达到试件最大承载能力的,50%,左右时，试件的变形开始偏离线性，部分基元破坏。,（,3,）软化阶段。当达到最大载荷之后，使试件进一步变形的载荷越来越小，进入弱化阶段，直至试件产生宏观破坏。,曲线分析,N-S,曲线,岩石的剪切加载试验,试验内容：,岩石压缩加载的应力,-,应变曲线,破坏模式,声发射,第,1,步：模型建立及参数选择,模型：,试样模型尺寸,100mm100mm,，网个划分为,100100,个基元 。采用平面应力问题，整个加载过程通过位移加载方式。加载位移量,s=0.005mm.,力学性质参数如下表：,均质度,弹性模量,强度,Poisson,自重,3,50000,200,0.25,0,摩擦角,C/T,比,Max strain T,Max strain C,Strength criterion,30,10,1.5,200,M-C Criterion,表,1,基元的力学参数、相变准则,RFPA,数值模型,注：预先切出,25mm,的水平引导裂纹，垂直方向加围压，每步,0.5Mpa,，终值加到第五步，达到最大值,2.5Mpa,。,第,2,步：计算,几分钟后,第,3,步：计算结果分析,剪切破裂过程,曲线分析,N-S,曲线,AE,能量图,