RFPA2D渗流版培训教程PPT课件.ppt

初级培训教程 RFPA2D Flow渗流版 2020 3 26 大连力软科技有限公司 1 主要内容 RFPA Flow渗流版简介 岩石 岩体 渗流力学特性及基本原理介绍 RFPA Flow渗流版功能按钮与工具条介绍 RFPA Flow渗流版分析问题的流程 工程中有关渗流问题的RFPA模拟计算 2020 3 26 大连力软科技有限公司 2 一 RFPA Flow渗流版简介 2020 3 26 大连力软科技有限公司 3 RFPA Flow渗流版简介 RFPA Flow渗流版 是继RFPA基本版推出的 利用该版本可进行岩石 体 基本渗流特性的模拟研究 亦可进行水工中岩石 体 流固耦合问题的数值计算分析 如可开展水力水电工程中各类坝体渗流场 孔隙压力场 微震时空展布等特征的计算分析 对承压水煤层开采等工程中的突水等突发性地质灾害的预测预报 观看岩石水压致裂全过程 2020 3 26 大连力软科技有限公司 4 水压致裂演化全过程 2020 3 26 大连力软科技有限公司 5 几个有关渗流的事例 大坝基础破坏与相应的渗流场 破坏模式 声发射 水流矢量场 加载过程中的载荷 渗透性及声发射演化 2020 3 26 大连力软科技有限公司 6 二 岩石 岩体 渗流力学特性及基本原理介绍 2020 3 26 大连力软科技有限公司 7 岩石 岩体 渗流力学特性 岩石作为天然地质材料 内部赋存着大量的孔隙 裂隙 这些缺陷结构的存在不但改变了岩石的力学性质 而且严重影响着岩体的渗透特性 岩石是低渗透性的材料 可以认为其不透水 但破裂失稳过程中随着裂纹的萌生 扩展和贯通 其渗透率发生重大的变化 一些工程事故中 如矿井突水 坝体溃决等 还有水压致裂问题 岩体破坏前是不透水的 或低渗透的 渗流应力耦合效应不明显 但应力和水压诱发新的损伤过程中 随着裂纹的扩展 应力场 渗流场以及损伤演化过程调整 迁移以及相互作用十分显著 并将导致岩体失稳 由此可见 研究岩石块体的力学性质 渗透性质对于工程岩体的渗流稳定性分析 设计具有重要的理论意义和实用价值 2020 3 26 大连力软科技有限公司 8 岩石 岩体 的孔隙特征 a 良好顺序的沉积矿床的高孔隙结构 b 由矿物碎屑充满的低孔隙结构 岩石的孔隙性通常用孔隙率n表示 对于大裂隙不发育的砂岩 蚀变花岗岩 风化带岩体及超固结的岩体 孔隙性的描述一般与土体大致相同 可按多孔介质处理 描述孔隙性的指标为孔隙率 即式中 n为岩石的孔隙率 Vp为岩石中孔隙 裂隙所占的体积 V为岩石的总体积 孔隙率是衡量岩石工程质量的重要物理性质指标之一 岩石的孔隙率反映了孔隙和裂隙在岩石中所占的百分率 孔隙率越大 岩石中的孔隙和裂隙就越多 岩石的力学性能就越差 2020 3 26 大连力软科技有限公司 9 岩石 岩体 的渗透特征 岩体中的水只能沿连通孔隙或裂缝渗透 岩石的渗透性大小可用渗透系数衡量 它主要决定于岩石孔隙的大小 方向及其相互连通情况 渗透率 permeability 是岩体介质特征的函数 它描述了岩体介质的一种平均性质 表示岩体介质传导流体的能力 对于均匀各向同性多孔介质而言 其渗透率为 式中 k为多孔介质渗透率 d为岩土体颗粒的有效粒径 c为比例常数 2020 3 26 大连力软科技有限公司 10 渗流 损伤耦合方程 当单元的应力状态或者应变状态将满足某个给定的损伤阈值时 单元开始损伤 损伤单元的弹性模量为 式中 D为损伤变量 E和Eo分别为损伤单元和无损伤单元的弹性模量 这些参数假定都是标量 下面以单轴压缩和拉伸本构关系为例 介绍单元的渗透率 损伤耦合方程 2020 3 26 大连力软科技有限公司 11 单元渗透率 损伤耦合方程 2020 3 26 大连力软科技有限公司 12 单元渗透率 损伤耦合方程 2020 3 26 大连力软科技有限公司 13 三 RFPA Flow渗流版功能按钮与工具条介绍 2020 3 26 大连力软科技有限公司 14 渗流边界条件 动态问题 2020 3 26 大连力软科技有限公司 15 基本工具条 New新建 Open打开 Save保存 Print打印 Copy复制 Zoom缩放 FullScreen全屏显示 Meshdividing网格划分 Boundaryconditions边界条件 Controlconditions控制条件 Selfdefining系统自定义 Seepageboundaryconditions渗流边界条件 Runstepbystep单步运行 Autorun自动连续运行 Stoprun强制停止运行 2020 3 26 大连力软科技有限公司 16 作图工具条 Drawbymouse鼠标键盘输入切换 Materialproperty材料类型选择 Selectarea区域选择 Rectangle矩形 Circle圆 Line单直线 Doubleline双直线 Polyline多线 Polygon任意多边形 Compositematerial复合材料 ClosedArc封闭弧 Arc圆弧 Anchor支护 复合材料增强介质颗粒设置 2020 3 26 大连力软科技有限公司 17 显示工具条 Shearstress剪应力图 Max principalstress最大主应力图 Min principalstress最小主应力图 Elastic弹性模量图 AcousticEmission声发射图 e AE弹模和声发射图 Load stepcurve载荷与加载步曲线 AEcurve声发射曲线 AreaAE圈定声发射图 Multi ElementInformation多单元信息 Slideplay幻灯播放 Displacementvector位移矢量图 Stressvector应力矢量图 Seepageinformation渗流多单元信息 Flowvolume流量分布图 Waterhead水位等值线分布图 Securitycurve安全系数曲线 2020 3 26 大连力软科技有限公司 18 Weibull分布函数 网格划分 YLength mm Xlength m Y方向和X方向的尺寸均为所研究问题的实际尺寸 单位为毫米 建立模型的时候需要将实际尺寸换算成毫米 单元尺寸为实际尺寸除单元个数 Rows Elements Cols Elements 单元列数和行数 单元尺寸的划分原则前面已讲述 HeterogeneityIndex行 为均质度 weibull统计分布函数中的参数m 反映岩石介质的均质性 若模拟地下工程等必须考虑模型本身的重力的工程问题时 需要输入自重的参数 注意自重的单位 N mm3 MeanValue 为单元物理力学参数的平均值 2020 3 26 大连力软科技有限公司 19 渗透特性设置 孔隙压力系数 poro waterpressure 是岩石本身的固有特性 值的大小为0到1之间 当单元完全破坏 成为空气元的时候 该数变成1 即孔隙压力终值 损伤突变系数 damagemutaioncoeff 是指基元体一旦有损伤 则其渗透性要提高的倍数 分离系数 separationcoeff 是指基元体在完全破坏后 其渗透性提高的倍数 耦合系数 couplingcoeff 表征的是应力对材料渗透特性的影响程度 渗透角度暂时无效 水平渗透系数 垂直渗透系数 孔隙率 孔隙压力系数 损伤突变系数 分离系数 耦合系数 孔隙压力终值 2020 3 26 大连力软科技有限公司 20 渗流边界条件设置 如果涉及渗流计算问题 则必须初始化模型的渗流边界条件 主要操作 先按下 再按下渗流边界类型 选择渗流边界条件 然后用鼠标在模型边界拉线 鼠标单击结束后 上述对话框中所设置的值就赋予了所选中的边界单元 既然是边界条件 所以注意不要将其设如模型内部 注意 压力用水头表示 在计算时候要注意单位的换算 1MPa约等于100米水头高 Increment增量正值为增加 负值为减小 2020 3 26 大连力软科技有限公司 21 渗流时间选项 如果计算需要考虑渗流时间问题 需要进行此项设置 时间选项 InterativeTimes 是计算的迭代次数 这个你保持默认值就行 默认值是2 TimesSpan Day 是时间尺度 即在非稳态计算的时候 每一个计算步所代表的时间 这个计算步就是rfpa控制中的步 我们非稳态计算采用的时间差分法 InitialPressure MPa 是整个模型中初始压力 各个单元的 2020 3 26 大连力软科技有限公司 22 控制边界条件设置 执行完网格划分命令就要对用户的模拟区进行控制条件赋值 控制条件涉及的是一些基本条件如 求解的类型是平面应力还是平面应变问题 加载总步数 是否为流固耦合等问题以及重力方向等如下图所示 其中Cavityelement和Equivalentforce两种计算方法 一般计算采用空单元法 如果考虑的问题和分步开挖加支护的话 选择等效节点力法 Pressurewater承压水和Phreatic潜水 是地下水两种不同的埋藏类型 承压水 pressurewater 埋藏在上下两个稳定隔水层之间 潜水 phreaticwater 埋藏在第一个稳定隔水层之上 2020 3 26 大连力软科技有限公司 23 系统自定义 selfdefine 2020 3 26 大连力软科技有限公司 24 系统自定义 selfdefine 图片类型和尺寸 PictureTypeSize 结果图类型 Type 此命令主要是定义在模型计算完成后要显示图形的类型和尺寸 根据研究问题的需要将相应的结果图选项上钩上即可 图形尺寸 Size 对系统生成的结果图进行大小设置 选Default默认 系统根据窗口大小随机生成 矢量单元间隔 Vectorelementintervals 因为要表示矢量场 用箭头表示矢量的大小方向 箭头方向为矢量方向 箭头长短表示矢量大小 位移矢量放大系数 Displacementvectoramplifiedcoeff 为了计算结果显示的我们能更好的看清以说明问题 对原始矢量大小给予适当的放大 应力矢量放大系数 Stressvectoramplifiedcoeff 解释同位移放大系数 注 这些参数定义可在计算运行命令之前执行 也可在计算命令之后执行 2020 3 26 大连力软科技有限公司 25 显示步 Showstep 显示步 Showstep 此命令在设置重画计算结果图时候设置 设置重画计算步的间隔 也可在inputredrawstep下面直接输入要重画的计算步 用英文状态下的逗号相隔 重画 Redraw 是系统对计算数据进行重新生成图像的操作 原计算数据不变化 只是新生成的结果图将覆盖和上次生成的名称相同的图形 不同于系统重置 reset 系统自定义 selfdefine 2020 3 26 大连力软科技有限公司 26 声发射图设置 AEPictureSetting 声发射圆大小 Typesize 设置声发射圆大小 代表相对能量或震级大小 Amplifiedradius是指在原能量圆大小的基础上放大的倍数 所有步 Forallsteps 显示当前所有步的声发射圆 当前步 Forcurrentstep 显示当前单步的声发射圆 发射圆颜色 AEcolor 用声发射圆颜色进行区别单元的破坏是抗拉破坏 抗压破坏 还是抗剪破坏 双击颜色框进行调节 系统自定义 selfdefine 2020 3 26 大连力软科技有限公司 27 高级设置 AdvancedSetting 主要对破坏后和空气单元的参数进行调整 初学者将词项默认 分离系数 Detachedcoeff 主要为了处理破坏后单元的显示状态 当破坏单元的变形 受拉 达到分离系数指定的变形量时 认为单元已彻底分离 将不在显示此单元 值为1 100 设置为5时 表示当单元的最大变形超过5倍单元原始尺寸时 不再显示此单元 灰度系数 Graydegreecoeff 是通过去掉高亮度单元来调整总体的亮度范围 值在0 1之间 其意义为 假设模型单元数为10 000个 灰度系数为0 01 则表示在显示时去掉10 000 0 01 100个最亮的个单元级别 大位移 Largedisplacement 是在处理自重问题时 当破坏的单元下落的位移超过设置值时 则认为破坏单元脱离母体单元 值为单元边长的倍数 系统自定义 selfdefine 2020 3 26 大连力软科技有限公司 28 变形放大系数 Deformation 对破模型的尺寸在显示时做放大处理 如有时计算完成后发现模型没有变化则可能是放大系数设置的太小 反之 当计算完成后显示的图形出现全黑时 则有可能是因为放大系数设置过大的缘故 计算精度系数 Accurateofstress 表示控制计算精度 如设置为0 0001时 则表示当破坏单元达到总单元数万分之一时 系统不再做当前不的循环计算 系统自定义 selfdefine 其他设置 Otherssetting 2020 3 26 大连力软科技有限公司 29 四 RFPA Flow渗流版分析问题流程 2020 3 26 大连力软科技有限公司 30 RFPA计算步骤 三大步 构造模型 前处理 1建立工程文档2网格划分和参数赋值3插入模型结构4边界条件 力 位移 渗流 5控制条件 选择渗流问题 计算单步运行连续运行 结果分析 后处理 选择渗流图 2020 3 26 大连力软科技有限公司 31 五 工程中几个渗流问题的RFPA模拟 2020 3 26 大连力软科技有限公司 32 孔隙水压力作用下岩石破裂过程RFPA模拟 水压致裂岩石过程的RFPA模拟 RFPA Flow工程算例 渗流模块 2020 3 26 大连力软科技有限公司 33 空隙水压下岩石破裂过程RFPA模拟 试验内容 考虑孔压作用下岩石的应力 应变曲线 岩石损伤破坏中渗透性的演化规律 2020 3 26 大连力软科技有限公司 34 孔隙压力对岩石破裂的影响 对于孔隙和微裂隙中含有重力水的岩石 当其突然受载水来不及排出时 岩石孔隙或裂隙中将产生很高的孔隙压力 这种孔隙压力 减小了岩石颗粒之间的压应力 从而降低了岩石的抗剪强度 甚至使岩石的微裂隙端处由于受拉状态从而破坏岩石的连接 2020 3 26 大连力软科技有限公司 35 第1步 模型建立及参数选择 模型 数值模型采用二维平面应力薄板模型 试样模型尺寸80mm 50mm 网个划分为160 100个基元 整个加载过程通过位移加载方式 轴向加载位移增量为 s 0 01mm 侧压p2 4Mpa 上下边界孔隙压力p3 p4分别为2 3Mpa 3 8Mpa 左右边界孔隙压力为0 控制步数为100步 力学性质参数如下表 表1强度准则和相变准则 2020 3 26 大连力软科技有限公司 36 RFPA数值模型 孔隙压力p4 3 8Mpa 孔隙压力p3 2 3Mpa 侧压p2 4Mpa 侧压p2 4Mpa 轴向位移加载 s 0 01mm 2020 3 26 大连力软科技有限公司 37 第2步 计算 几分钟后 2020 3 26 大连力软科技有限公司 38 第3步 计算结果分析 2020 3 26 大连力软科技有限公司 39 2020 3 26 大连力软科技有限公司 40 2020 3 26 大连力软科技有限公司 41 2020 3 26 大连力软科技有限公司 42 2020 3 26 大连力软科技有限公司 43 2020 3 26 大连力软科技有限公司 44 2020 3 26 大连力软科技有限公司 45 2020 3 26 大连力软科技有限公司 46 2020 3 26 大连力软科技有限公司 47 2020 3 26 大连力软科技有限公司 48 N S曲线 Seepage曲线 AE曲线 2020 3 26 大连力软科技有限公司 49 岩石水压致裂过程RFPA数值模拟 试验内容 含单孔围压作用下岩石水压致裂