RFPA2D渗流版培训教程.pptx

主要内容主要内容 RFPA-FlowRFPA-Flow渗流版简介渗流版简介渗流版简介渗流版简介 岩石岩石岩石岩石（岩体）（岩体）（岩体）（岩体）渗流力学特性及基本原理介绍渗流力学特性及基本原理介绍渗流力学特性及基本原理介绍渗流力学特性及基本原理介绍 RFPA-FlowRFPA-Flow渗流版功能按钮与工具条介绍渗流版功能按钮与工具条介绍渗流版功能按钮与工具条介绍渗流版功能按钮与工具条介绍 RFPA-FlowRFPA-Flow渗流版分析问题的流程渗流版分析问题的流程渗流版分析问题的流程渗流版分析问题的流程 工程中有关渗流问题的工程中有关渗流问题的工程中有关渗流问题的工程中有关渗流问题的RFPARFPA模拟计算模拟计算模拟计算模拟计算一、一、RFPA-Flow渗流版简介渗流版简介RFPA-Flow渗流渗流版简介版简介 RFPA-Flow渗流版渗流版:是继是继是继是继RFPARFPA基本版推出的，利用该版本可进行基本版推出的，利用该版本可进行基本版推出的，利用该版本可进行基本版推出的，利用该版本可进行岩石岩石岩石岩石（体）（体）（体）（体）基本渗流特性的模拟研究，亦可进行水工中岩石基本渗流特性的模拟研究，亦可进行水工中岩石基本渗流特性的模拟研究，亦可进行水工中岩石基本渗流特性的模拟研究，亦可进行水工中岩石（体）（体）（体）（体）流流流流固耦合问题的数值计算分析，如可开展水力水电工程中各类坝体固耦合问题的数值计算分析，如可开展水力水电工程中各类坝体固耦合问题的数值计算分析，如可开展水力水电工程中各类坝体固耦合问题的数值计算分析，如可开展水力水电工程中各类坝体渗流场、孔隙压力场、微震时空展布等特征的计算分析，对承压渗流场、孔隙压力场、微震时空展布等特征的计算分析，对承压渗流场、孔隙压力场、微震时空展布等特征的计算分析，对承压渗流场、孔隙压力场、微震时空展布等特征的计算分析，对承压水煤层开采等工程中的突水等突发性地质灾害的预测预报。水煤层开采等工程中的突水等突发性地质灾害的预测预报。水煤层开采等工程中的突水等突发性地质灾害的预测预报。水煤层开采等工程中的突水等突发性地质灾害的预测预报。观看岩石水压致裂全过程观看岩石水压致裂全过程水压致裂演化全过程水压致裂演化全过程几个有关渗流的事例几个有关渗流的事例大坝基础破坏与相应的渗流场大坝基础破坏与相应的渗流场破坏模式破坏模式声发射声发射水流矢量场水流矢量场加载过程中的载荷、渗透性及声发射演化加载过程中的载荷、渗透性及声发射演化二、岩石二、岩石（岩体）（岩体）渗流力学特性渗流力学特性 及基本原理介绍及基本原理介绍岩石岩石（岩体）（岩体）渗流力学特性渗流力学特性岩石作为天然地质材料，内部赋存着大量的孔隙、裂隙，这些缺陷结岩石作为天然地质材料，内部赋存着大量的孔隙、裂隙，这些缺陷结构的存在不但改变了岩石的力学性质，而且严重影响着岩体的构的存在不但改变了岩石的力学性质，而且严重影响着岩体的 渗透特渗透特性。性。岩石是低渗透性的材料，可以认为其不透水，但破裂失稳过程中随着岩石是低渗透性的材料，可以认为其不透水，但破裂失稳过程中随着裂纹的萌生、扩展和贯通，其渗透率发生重大的变化。一些工程事裂纹的萌生、扩展和贯通，其渗透率发生重大的变化。一些工程事 故故中，如矿井突水、坝体溃决等，还有水压致裂问题，岩体破坏前是不透中，如矿井突水、坝体溃决等，还有水压致裂问题，岩体破坏前是不透水的水的（或低渗透的），渗流应力耦合效应不明显。但应力和水压诱发新（或低渗透的），渗流应力耦合效应不明显。但应力和水压诱发新的损伤过程中，随着裂纹的扩展，应力场、渗流场以及损伤演化过程调的损伤过程中，随着裂纹的扩展，应力场、渗流场以及损伤演化过程调整、迁移以及相互作用十分显著，并将导致岩体失稳。由此可见，研究整、迁移以及相互作用十分显著，并将导致岩体失稳。由此可见，研究岩石块体的力学性质、渗透性质对于工程岩体的渗流稳定性分析、设计岩石块体的力学性质、渗透性质对于工程岩体的渗流稳定性分析、设计具有重要的理论意义和实用价值。具有重要的理论意义和实用价值。岩石岩石（岩体）（岩体）的孔隙特征的孔隙特征(a)良好顺序的沉积矿床的高孔隙结构良好顺序的沉积矿床的高孔隙结构(b)由矿物碎屑充满的低孔隙结构由矿物碎屑充满的低孔隙结构 岩石的孔隙性通常用孔隙率岩石的孔隙性通常用孔隙率n表示。对于大裂隙不表示。对于大裂隙不发育的砂岩、蚀变花岗岩、风化带岩体及超固结的岩体，发育的砂岩、蚀变花岗岩、风化带岩体及超固结的岩体，孔隙性的描述一般与土体大致相同，可按多孔介质处理，孔隙性的描述一般与土体大致相同，可按多孔介质处理，描述孔隙性的指标为孔隙率，即描述孔隙性的指标为孔隙率，即式中：式中：n为岩石的孔隙率；为岩石的孔隙率；Vp为岩石中孔隙、为岩石中孔隙、裂隙所裂隙所占的体积；占的体积；V为岩石的总体积。为岩石的总体积。孔隙率是衡量岩石工程质量的重要物理性质指标之孔隙率是衡量岩石工程质量的重要物理性质指标之一。岩石的孔隙率反映了孔隙和裂隙在岩石中所占的百一。岩石的孔隙率反映了孔隙和裂隙在岩石中所占的百分率，孔隙率越大，岩石中的孔隙和裂隙就越多，岩石分率，孔隙率越大，岩石中的孔隙和裂隙就越多，岩石的力学性能就越差。的力学性能就越差。岩石岩石（岩体）（岩体）的渗透特征的渗透特征 岩体中的水只能沿连通孔隙或裂缝渗透。岩石的渗透性大小可用渗透系数衡量，它主要决定于岩体中的水只能沿连通孔隙或裂缝渗透。岩石的渗透性大小可用渗透系数衡量，它主要决定于岩石孔隙的大小、方向及其相互连通情况。岩石孔隙的大小、方向及其相互连通情况。渗透率渗透率(permeability)是岩体介质特征的函数，它描述了岩体介质的一种平均性质，表示岩体介是岩体介质特征的函数，它描述了岩体介质的一种平均性质，表示岩体介质传导流体的能力。对于均匀各向同性多孔介质而言，其渗透率为：质传导流体的能力。对于均匀各向同性多孔介质而言，其渗透率为：式中：式中：k为多孔介质渗透率；为多孔介质渗透率；d为岩土体颗粒的有效粒径；为岩土体颗粒的有效粒径；c为比例常数。为比例常数。渗流渗流-损伤耦合方程损伤耦合方程 当单元的应力状态或者应变状态将满足某个给定的损伤阈值时，单元开始损当单元的应力状态或者应变状态将满足某个给定的损伤阈值时，单元开始损伤，损伤单元的弹性模量为：伤，损伤单元的弹性模量为：式中：式中：D为损伤变量；为损伤变量；E和和Eo分分别为损伤单元和无损伤单元的弹性模量，这些参数假定都是别为损伤单元和无损伤单元的弹性模量，这些参数假定都是标量。标量。下面以单轴压缩和拉伸本构关系为例，介绍单元的渗透率下面以单轴压缩和拉伸本构关系为例，介绍单元的渗透率-损伤耦合方程。损伤耦合方程。单元渗透率单元渗透率-损伤耦合方程损伤耦合方程单元渗透率单元渗透率-损伤耦合方程损伤耦合方程三、三、RFPA-Flow渗流版功能按钮渗流版功能按钮与工具条介绍与工具条介绍渗流边界条件渗流边界条件动态问题动态问题基基本本工工具具条条New新建新建Open打开打开Save保存保存Print打印打印Copy复制复制Zoom缩放缩放Full Screen全屏显示全屏显示Mesh dividing网格划分网格划分Boundary conditions边界条件边界条件Control conditions控制条件控制条件Self defining系统自定义系统自定义Seepage boundary 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应力矢量图应力矢量图Seepage information 渗流多单元信息渗流多单元信息Flow volume流量分布图流量分布图Water head水位等值线分布图水位等值线分布图Security curve安全系数曲线安全系数曲线Weibull分布函数：分布函数：网格划分网格划分 Y Length(mm),X length(m):Y方向和方向和X方向的尺寸均为所研究问题的实际尺寸，单位为毫米，建立模型的时方向的尺寸均为所研究问题的实际尺寸，单位为毫米，建立模型的时候需要将实际尺寸换算成毫米，单元尺寸为实际候需要将实际尺寸换算成毫米，单元尺寸为实际 尺寸除单元个数。尺寸除单元个数。Rows(Elements),Cols(Elements):单元列数和行数，单元尺寸的划分原则前面已讲述。单元列数和行数，单元尺寸的划分原则前面已讲述。Heterogeneity Index行：行：为均质度，为均质度，weibull统计分布函数中的参数统计分布函数中的参数m，反映岩石介质的均质性。若模拟地下，反映岩石介质的均质性。若模拟地下工程等必须考虑模型本身的重力的工程问题时，需要输入自重的参数，注意自重的单位：工程等必须考虑模型本身的重力的工程问题时，需要输入自重的参数，注意自重的单位：N/mm3。Mean Value:为单元物理力学参数的平均值。为单元物理力学参数的平均值。渗透特性设置渗透特性设置孔隙压力系数孔隙压力系数(poro-water pressure)：是岩石本身的固有特性，值的大小为是岩石本身的固有特性，值的大小为0到到1之间；当单元之间；当单元完全破坏，成为空气元的时候，该数变成完全破坏，成为空气元的时候，该数变成1，即孔隙压力终值。，即孔隙压力终值。损伤突变系数损伤突变系数(damage mutaion coeff)：是指基元体一旦有损伤，则其渗透性要提高的倍数。是指基元体一旦有损伤，则其渗透性要提高的倍数。分离系数分离系数(separation coeff)：是指基元体在完全破坏后，其渗透性提高的倍数。是指基元体在完全破坏后，其渗透性提高的倍数。耦合系数耦合系数(coupling coeff)：表征的是应力对材料渗透特性的影响程度。（渗透角度暂时无效）表征的是应力对材料渗透特性的影响程度。（渗透角度暂时无效）水平渗透系数垂直渗透系数孔隙率孔隙压力系数损伤突变系数分离系数耦合系数孔隙压力终值渗流边界条件设置渗流边界条件设置如果涉及渗流计算问题，则必须初始化模型的渗流边界条件，主要操作：先按下如果涉及渗流计算问题，则必须初始化模型的渗流边界条件，主要操作：先按下再按下渗流边界类型再按下渗流边界类型，选择渗流边界条件：，选择渗流边界条件：然后用鼠标在模型边界拉线，鼠标单击结束后，上述对话框中所设置的值就赋予了所选中的边界单元（既然后用鼠标在模型边界拉线，鼠标单击结束后，上述对话框中所设置的值就赋予了所选中的边界单元（既然是边界条件，所以注意不要将其设如模型内部）。注意：压力用水头表示，在计算时候要注意单位的换然是边界条件，所以注意不要将其设如模型内部）。注意：压力用水头表示，在计算时候要注意单位的换算，算，1MPa约等于约等于100米水头高。米水头高。Increment增量正值为增加，负值为减小。增量正值为增加，负值为减小。渗流时间选项渗流时间选项:如果计算需要考虑渗流时间问题如果计算需要考虑渗流时间问题,需要进行此项设置：需要进行此项设置：时间选项时间选项Interative Times：是计算的迭代次数，这个你保持默认值就行，默认值是是计算的迭代次数，这个你保持默认值就行，默认值是2。Times Span(Day)：是时间尺度，即在非稳态计算的时候，每一个计算步所代表的是时间尺度，即在非稳态计算的时候，每一个计算步所代表的时间，这个计算步就是时间，这个计算步就是rfpa控制中的步，我们非稳态计算采用的时间差分法。控制中的步，我们非稳态计算采用的时间差分法。Initial Pressure(MPa)：是整个模型中初始压力，各个单元的。是整个模型中初始压力，各个单元的。控制边界条件设置控制边界条件设置执行完网格划分命令就要对用户的模拟区进行控制条件赋值。控制条件涉及的是一些基本条件如：求解的执行完网格划分命令就要对用户的模拟区进行控制条件赋值。控制条件涉及的是一些基本条件如：求解的类型是平面应力还是平面应变问题；加载总步数、是否为流固耦合等问题以及重力方向等如下图所示：类型是平面应力还是平面应变问题；加载总步数、是否为流固耦合等问题以及重力方向等如下图所示：其中其中Cavity element 和和 Equivalent force两种计算方法，一般计算采用空单元法，如果考虑两种计算方法，一般计算采用空单元法，如果考虑的问题和分步开挖加支护的话，选择等效节点力法。的问题和分步开挖加支护的话，选择等效节点力法。Pressure water承压水和承压水和Phreatic潜水，是地下潜水，是地下水两种不同的埋藏类型水两种不同的埋藏类型。承压水承压水(pressure water)：埋藏在上下两个稳定埋藏在上下两个稳定隔水层之间隔水层之间。潜水潜水(phreatic water)：埋藏在第一个稳定隔水埋藏在第一个稳定隔水层之上层之上。系统自定义系统自定义 (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-Flow渗流版分析问题流程渗流版分析问题流程RFPA计算步骤计算步骤（三大步）（三大步）构造模型构造模型 （前处理）（前处理）1建立工程文档建立工程文档2网格划分和参数赋值网格划分和参数赋值3插入模型结构插入模型结构4边界条件边界条件(力、位移、渗流力、位移、渗流)5控制条件控制条件(选择渗流问题）选择渗流问题）计算计算单步运行连续运行结果分析结果分析（后处理）（后处理）选择渗流图选择渗流图五、工程中几个渗流问题的五、工程中几个渗流问题的RFPA模拟模拟 孔隙水压力作用下岩石破裂过程孔隙水压力作用下岩石破裂过程孔隙水压力作用下岩石破裂过程孔隙水压力作用下岩石破裂过程RFPARFPA模拟模拟模拟模拟 水压致裂岩石过程的水压致裂岩石过程的水压致裂岩石过程的水压致裂岩石过程的RFPARFPA模拟模拟模拟模拟RFPA-Flow工程算例工程算例（渗流模块）（渗流模块）空隙水压下岩石破裂过程空隙水压下岩石破裂过程RFPA模拟模拟试验内容：试验内容：考虑孔压作用下岩石的应力考虑孔压作用下岩石的应力-应变曲线应变曲线 岩石损伤破坏中渗透性的演化规律岩石损伤破坏中渗透性的演化规律孔隙压力对岩石破裂的影响孔隙压力对岩石破裂的影响 对于孔隙和微裂隙中含有重力水的岩石，当其突然受载水来不及对于孔隙和微裂隙中含有重力水的岩石，当其突然受载水来不及排出时，岩石孔隙或裂隙中将产生很高的孔隙压力。这种孔隙压排出时，岩石孔隙或裂隙中将产生很高的孔隙压力。这种孔隙压力，减小了岩石颗粒之间的压应力，从而降低了岩石的抗剪强度，力，减小了岩石颗粒之间的压应力，从而降低了岩石的抗剪强度，甚至使岩石的微裂隙端处由于受拉状态从而破坏岩石的连接。甚至使岩石的微裂隙端处由于受拉状态从而破坏岩石的连接。第第1步：步：模型建立及参数选择模型建立及参数选择 模型：模型：数值模型采用二维平面应力薄板模型。试样模型尺寸数值模型采用二维平面应力薄板模型。试样模型尺寸80mm50mm，网个划分为，网个划分为160100个基元个基元。整个加载过程通过位移加载方式。轴向加载位移增量为。整个加载过程通过位移加载方式。轴向加载位移增量为s=0.01mm。侧压。侧压p2=4Mpa，上下边界孔隙压力，上下边界孔隙压力p3,p4分别为分别为2.3Mpa,3.8Mpa。左。左右边界孔隙压力为右边界孔隙压力为0。控制步数为。控制步数为100步。步。力学性质参数如下表：力学性质参数如下表：均质度弹性模量强度PoissonC/T 2100001000.2510摩擦角孔隙率渗透系数孔隙水压系数耦合系数300.10.10.60.1表表 1 强度准则和相变准则强度准则和相变准则RFPA数值模型数值模型孔隙压力孔隙压力p4=3.8Mpa孔隙压力孔隙压力p3=2.3Mpa侧压侧压p2=4Mpa侧压侧压p2=4Mpa轴向位移加载轴向位移加载s=0.01mm第第2步：步：计算计算几分钟后几分钟后第第3步：步：计算结果分析计算结果分析N-S曲线、曲线、Seepage曲线、曲线、AE曲线曲线 岩石水压致裂过程岩石水压致裂过程RFPA数值模拟数值模拟试验内容：试验内容：含单孔围压作用下岩石水压致裂过程含单孔围压作用下岩石水压致裂过程 从孔心到边界应力分布情况从孔心到边界应力分布情况 不同均质度时孔边应力分布的指纹图不同均质度时孔边应力分布的指纹图 不同边界压力时破坏裂纹模式、应力分布、流量分布情况不同边界压力时破坏裂纹模式、应力分布、流量分布情况第第1步：步：模型建立及参数选择模型建立及参数选择 模型：模型：为了研究不同均质度时孔边应力分布的光弹图，首先建立一个均匀材料的模型（均质度为为了研究不同均质度时孔边应力分布的光弹图，首先建立一个均匀材料的模型（均质度为100）材）材料力学参数如下表，再建立一个非均匀材料的模型（均质度为料力学参数如下表，再建立一个非均匀材料的模型（均质度为3）。数值模型采用二维平面应力薄）。数值模型采用二维平面应力薄板模型。试样模型尺寸板模型。试样模型尺寸100mm100mm，孔径为，孔径为5mm。网个划分为。网个划分为100100个基元个基元。整个加载过。整个加载过程通过压力加载方式。孔压增量为程通过压力加载方式。孔压增量为p=1Mpa。水平压力。水平压力ph=1.0Mpa，垂直压力，垂直压力 pv=4.0Mpa，孔内，孔内初始压力为初始压力为1Mpa.控制步数为控制步数为50步。力学性质参数如下表：步。力学性质参数如下表：均质度弹性模量强度PoissonC/T 3（100）40000600.2510摩擦角孔隙率渗透系数孔隙水压系数耦合系数300.10.010.60.01表表 1 强度准则和相变准则强度准则和相变准则RFPA数值模型数值模型Ph=1.0MpaPv=4.0Mpa孔内初始压力为孔内初始压力为1.0Mpa第第2步：步：计算计算几分钟后几分钟后第第3步：步：计算结果分析计算结果分析应力分布指纹图应力分布指纹图水压等值线图水压等值线图水流量图水流量图 m=10，计算到第，计算到第10步的模拟结果图步的模拟结果图从孔心到边界应力分布情况从孔心到边界应力分布情况