PFC(6.0)单轴模拟花岗岩等矿物晶体岩石：基于Potyondy D.2010的GBM模型与...

PFC(6.0)模拟:GBM模型grain- based model 参考文献Potyondy. D. 2010) pb-sj或pb-pb 单轴。 模拟花岗岩等矿物晶体岩石多种矿物晶体模型其中矿物种类 数量分布可以自定义。可以监测sj裂纹和各矿物内裂纹。一、文档概述本文档针对岩石矿物GBMGrain-Based Model基于颗粒的模型的PFCParticle Flow Code模拟代码进行功能解析。该代码体系共包含7个核心文件遵循Potyondy. D. 2010提出的GBM模型理论通过多步骤的颗粒生成、几何建模、力学参数赋值与加载模拟实现岩石矿物在受力条件下的力学响应与断裂行为模拟为岩石力学领域的数值分析提供标准化的计算流程。二、核心功能模块拆解一基础模型初始化模块01.dat该模块是整个模拟的“地基”核心功能是构建初始的颗粒系统与力学环境为后续GBM模型的颗粒分组与几何转化奠定基础具体功能包含以下三方面计算域与边界定义通过domain extent -1 1 condition destroy设定计算域范围并采用wall generate命令生成边界墙体墙体覆盖范围在-0.05 0.05x轴与-0.1 0.1y轴之间且通过expand 1.5参数扩展边界避免颗粒在运动过程中超出计算域。接触模型与参数预设采用线性接触模型linear method deform作为默认接触模型设置弹性模量emod1.0e9 Pa、刚度比kratio0.5同时定义法向刚度比dp_nratio0.5确保颗粒间接触力学行为符合线性弹性规律。多矿物颗粒生成与分组按照岩石矿物的实际组分比例通过ball distribute命令生成四种不同矿物的颗粒各矿物颗粒的半径范围、体积占比与分组对应关系明确且所有颗粒均分布在-0.05 0.05x轴与-0.1 0.1y轴的指定区域内同时设置颗粒密度为2500 kg/m³、阻尼系数为0.7模拟真实岩石的物理属性。系统平衡与模型保存通过model cycle 1000 calm 10完成系统初始平静化消除颗粒初始叠加产生的不稳定力再通过model solve ratio-average 1e-2将系统力比收敛至1e-2确保颗粒间受力平衡最终以model save model1保存初始模型为后续模块调用提供基础数据。二颗粒-几何节点映射模块02.dat该模块是连接“颗粒系统”与“几何模型”的关键桥梁核心功能是将01.dat中生成的分组颗粒转化为几何节点为后续Voronoi多面体模拟岩石晶粒的构建提供几何依据功能流程如下模型恢复与几何集合创建通过model restore model1调用初始颗粒模型同时创建名为“rock”的几何集合geom.set.create(rock)作为后续几何节点的统一存储容器。颗粒信息提取与节点转化通过循环遍历所有颗粒loop foreach bp ball.list提取每个颗粒的坐标ball.pos.x(bp)、ball.pos.y(bp)与所属矿物分组ball.group(bp)并基于颗粒坐标创建几何节点geom.node.create同时将颗粒的矿物分组信息同步赋给对应几何节点确保“颗粒分组”与“几何节点分组”一一对应。Voronoi多面体构建与模型保存基于“rock”几何集合中的节点通过rblock construct from-geometry rock voronoi命令构建Voronoi多面体这些多面体将作为模拟岩石中“晶粒”的几何形态最后以model save model2保存几何模型完成从“离散颗粒”到“连续晶粒”的形态转化。三Voronoi多面体矿物分组模块03.dat该模块的核心功能是为Voronoi多面体即“晶粒”赋予矿物属性通过颗粒与多面体的空间关联实现多面体的矿物分组具体功能包括模型恢复与多面体-颗粒关联调用model restore model2加载Voronoi多面体模型通过循环遍历所有多面体loop foreach rp rblock.list提取每个多面体的中心坐标rblock.ball.pos(rp)并找到该坐标附近最近的颗粒ball.near(rv)以颗粒的矿物分组作为多面体分组的依据。多面体矿物分组赋值根据最近颗粒的分组石英、钾长石、斜长石、云母将对应矿物分组名称赋给多面体rblock.group(rp) 矿物名称同时通过计数器minerala/mineralb/mineralc/minerald统计各矿物多面体的数量为后续组分比例验证提供数据支持。几何导出与冗余元素清理将不同矿物分组的多面体分别导出为对应的几何集合如rblock export to-geometry shiying随后删除原始的“rock”几何集合、Voronoi多面体与初始颗粒仅保留矿物分组后的几何集合以model save model3保存为后续细观颗粒填充做准备。四细观颗粒填充与分组模块04.dat该模块是GBM模型“细观化”的核心通过在矿物几何集合内填充细颗粒模拟真实岩石中“晶粒由更小颗粒组成”的结构具体功能如下模型恢复与接触参数重置加载model3中的矿物几何集合重新设置接触模型为线性接触参数与01.dat保持一致确保力学模型的连贯性同时重新生成边界墙体覆盖范围调整为-0.025 0.025x轴与-0.05 0.05y轴适配细颗粒填充区域。细颗粒生成与分布采用ball distribute命令生成细颗粒颗粒半径范围为0.2e-3~0.3e-3 m孔隙率控制在0.02接近真实岩石的致密结构且仅分布在上述指定区域内颗粒密度与阻尼系数延续01.dat的设置2500 kg/m³、0.7。细颗粒矿物分组赋值通过geometry-distance命令根据细颗粒与各矿物几何集合的空间距离将细颗粒分配至对应矿物分组如ball group shiying range geometry-distance shiying实现“细颗粒-矿物分组”的精准关联最后完成系统平衡力比收敛至1e-3并保存模型为model4为后续力学参数赋值做好准备。五力学参数赋值与系统平衡模块05.dat该模块是模拟岩石力学行为的“核心引擎”通过为不同矿物颗粒的接触与界面接触设置差异化力学参数还原真实岩石的非均质力学特性功能拆解如下模型恢复与边界优化加载model4的细颗粒模型删除多余的边界墙体仅保留关键承载墙体减少非必要的边界约束对模拟结果的干扰。接触分组定义通过循环遍历所有接触loop foreach local cp contact.list区分“颗粒-颗粒”接触的类型若接触的两个颗粒属于同一矿物分组将接触分组定义为“rock”代表晶粒内部接触若属于不同矿物分组定义为“jiemian”代表晶粒间界面接触为差异化参数赋值提供分类依据。界面与晶粒内部参数赋值针对“jiemian”接触晶粒界面采用光滑节理模型smoothjoint设置节理刚度、内摩擦角、抗拉强度与黏聚力等参数模拟界面较弱的力学特性针对“rock”接触晶粒内部按矿物类型分别设置力学参数不同矿物的弹性模量、抗拉强度、黏聚力等参数差异显著且均采用平节理模型flatjoint模拟晶粒内部较强的力学性能。系统平衡与初始状态清零通过model cycle 1000 calm 10与model solve ratio-average 1e-3完成系统二次平衡消除参数赋值后产生的局部应力集中随后将颗粒的位移、接触力与接触力矩清零multiply 0.0确保模拟起始状态为“无初始应力”最终保存模型为model5为加载试验模块提供“洁净”的力学模型。六加载试验与数据监测模块06.dat该模块是模拟岩石受力破坏过程的“执行端”通过施加轴向荷载监测岩石的应力-应变响应与断裂发展具体功能包括模型恢复与加载系统初始化加载model5的力学模型通过自定义函数setupwall识别上下承载墙体wptop、wpbottom计算初始试样高度sampleheight与横截面面积crosssectionalarea为应力计算提供几何参数。轴向加载控制采用位移控制方式通过wall attribute velocity-y命令为上下墙体施加反向速度上墙体向下、下墙体向上实现轴向压缩加载同时设置颗粒阻尼系数为0.1抑制加载过程中产生的数值振荡。关键数据监测通过fish history命令记录两个核心指标一是轴向应力由axialstresswall函数基于墙体接触力与横截面面积计算二是轴向应变由axialstrainwall函数基于墙体位移与试样高度计算形成应力-应变曲线的基础数据同时调用断裂监测脚本fracture.p2fis开启裂纹数量crack_num的实时监测。加载终止条件与模型保存设置加载终止条件为“应力降至峰值应力的30%”peakfraction0.3当满足该条件时通过model solve fish-halt自动终止加载模拟岩石从“峰值强度”到“残余强度”的全破坏过程最后输出峰值应力peakstress并保存模型为model6完成整个加载模拟流程。七断裂行为监测模块fracture.p2fis该模块是岩石断裂机制分析的“眼睛”通过监测颗粒间接触键的断裂事件实现裂纹的动态捕捉、分类与记录具体功能如下断裂监测初始化通过trackinit函数完成初始化工作包括删除历史断裂数据、清空碎片记录、注册颗粒-颗粒接触的断裂监测事件同时设置裂纹累积计数器crackaccum与裂纹总数计数器crack_num确保监测起点“归零”。裂纹捕捉与分类当颗粒间接触键发生断裂时触发bondbreak事件addcrack函数自动响应提取断裂接触的位置fracpos、法线方向norm与断裂模式拉伸或剪切根据断裂模式将裂纹分为“拉伸裂纹cracktension”与“剪切裂纹crackshear”并基于颗粒半径确定裂纹长度fracsize。裂纹几何与时间记录为每一条新生成的裂纹创建几何实体fracture.create记录裂纹的生成时间mech.time.total、关联颗粒编号fracture.extra并将裂纹归入对应模式的断裂集合dfn_label同时每累积50条裂纹自动更新一次裂纹位置确保裂纹几何形态与颗粒运动同步。碎片计算触发当裂纹累积达到50条时触发fragment compute命令计算当前试样的碎片数量与形态为岩石破坏后的碎片分布分析提供数据支持完整还原岩石从“连续体”到“离散碎片”的断裂演化过程。三、整体工作流程与逻辑关联整个代码体系遵循“从几何建模到力学模拟从静态构建到动态加载”的逻辑各模块间通过模型文件model1~model6实现数据传递形成闭环工作流程几何构建阶段01.dat→02.dat→03.dat从“离散颗粒生成”到“几何节点映射”再到“Voronoi多面体分组”逐步将岩石的矿物组分转化为几何形态完成GBM模型的“骨架”构建。细观填充阶段04.dat在几何骨架内填充细颗粒实现“宏观晶粒”到“细观颗粒”的尺度下探还原岩石的细观结构特征。力学赋能阶段05.dat为不同矿物的细颗粒接触与界面接触赋予差异化力学参数让模型具备“非均质力学特性”接近真实岩石的物理本质。加载与监测阶段06.dat→fracture.p2fis通过轴向加载模拟岩石受力过程同步监测应力-应变响应与断裂行为输出岩石力学性能峰值强度、弹性模量与断裂机制拉伸/剪切裂纹比例的核心数据为岩石力学研究提供数值支撑。四、应用价值与扩展方向一应用价值岩石力学研究可用于模拟不同矿物组分、不同细观结构的岩石在压缩、拉伸等荷载下的力学响应为岩石强度预测、破坏机制分析提供数值工具。工程设计支撑可为隧道开挖、矿岩开采等工程提供岩石稳定性分析依据通过模拟不同工程条件下岩石的破坏过程优化工程设计参数。二扩展方向多荷载条件扩展当前代码仅支持轴向压缩加载可进一步扩展为拉伸、弯曲、剪切等多荷载模式覆盖更全面的岩石力学试验场景。多物理场耦合可引入温度、渗流等物理场参数构建热力耦合、水力耦合的GBM模型模拟高温、高渗环境下岩石的力学行为。参数敏感性分析可基于该代码体系开展矿物弹性模量、界面黏聚力等参数的敏感性分析明确各参数对岩石力学性能的影响权重为模型参数优化提供依据。PFC(6.0)模拟:GBM模型grain- based model 参考文献Potyondy. D. 2010) pb-sj或pb-pb 单轴。 模拟花岗岩等矿物晶体岩石多种矿物晶体模型其中矿物种类 数量分布可以自定义。可以监测sj裂纹和各矿物内裂纹。