0 引言 前两天发了【PFC6.0】随机多边形区域划分及颗粒填充 的文章后，公众号后台收到了一些同学的留言。介绍了在PFC高版本软件中，可以直接使用泰森多边形进行rblock的生成。 于是乎我就趁这机会也升了一下目前使用的版本。上个文章采用的是PFC6.0.13，本文章采用PFC6.0.30进行测试。

50ball property "fric" 0.5model solvemodel save "sample"结果如图所示：2 clump替换 这里还是用的上一个案例【【PFC6.0】三维Cluster碎石三轴模拟】的clump模板，这个工况对模板的要求不高，可以自行生成模板。

在具体的实施中，需要在发生耦合作用的连续单元的表面上生成PFC模型中的墙单元，生成的墙单元做为离散元单元和连续单元之间耦合的介质。 在数值试验之前，采用PFC进行岩体三轴试验进行颗粒细观参数标定。三轴试验的样品尺寸为直径1m，高度2.5m。在PFC中生成9627个球体，颗粒的大小分布于0.005-0.1m之间，随机分布）。

一套完整的 PFC 6.0 离散元原创算例，专为研究岩石在复杂应力路径下的力学响应及岩爆（Rockburst）现象而设计。代码实现了从初始围压保载到分级径向加载，再到不同波形动力扰动的全过程模拟，逻辑严密，注释清晰。 代码集成了四种极具科研价值的加载工况，用户可一键切换： 分级静力加载：模拟深部岩体开挖过程中的应力重分布。

0 引言 前两天测试了GPT3.5书写PFC命令，结果是一塌糊涂，给人一种不知悔改且极度傲慢的感觉。于是这两天使用GPT4.0来测试一下PFC的书写能力，结果出人意料的好，不到10句的沟通就完成了，颗粒在墙体中自重沉降的算例。因为GPT4.0更加聪明，并且每三小时只有10句话的沟通机会，我上来直奔主题，提出自己大概的需求。

[-width*0.5] [width*0.5] [-width*0.5] [width*0.5] [-height*0.5] [height*0.5] ball attribute density 2e3 damp 0.7contact cmat default model linear method deform emod 100e6 kratio 1.5model cycle 2000 calm

pos-z [-zhuan_height*1.6] [zhuan_height*1.6] wall-zone create name 'bianjie3' skip-errors range pos-x [-zhuan_length*1.6] [zhuan_length*1.6] ...

box [-sample_width*0.5] [sample_width*0.5] [-sample_hight*0.5] [sample_hight*0.5]cmat default type ball-facet model linear method deform emod 100e6 kratio 1.5 cmat default type ball-ball model linear

property pb_coh 15e6 pb_ten 15e6 pb_fa 50 fric 0.1 lin_mode 1 range fish @isRock1In cmat add 3 model linearpbond method deform emod 6e9 kratio 1.5 pb_deform emod 6e9 kratio 1.5 ...

<p>pfc岩石标定<span style="color: rgb(18, 18, 18);">必备</span>模型试验，可以自行转,PFC6.0版本，抗压强度和巴西劈裂试验,内含微风化石灰岩参数</p><div contenteditable="false" width="100%"><figure class="figure-image" data-img="https://img.jishulink.com

model restore "rblock_slip"rblock export to-geometry "allgeo" cmat default type ball-facet model linear method deform emod 100e6 kratio 1.5 property fric 0.5cmat default type ball-ball model linear method

y_center=(y_max+y_min)*0.5 z_center=(z_max+z_min)*0.5 local gs = geom.set.find(id) loop foreach local gn geom.node.list(gs) geom.node.pos.x(gn)=-x_center+geom.node.pos.x(gn)

6.2盾构隧道掌子面稳定性(命令流+FISH)6.2.1主动失稳模式6.2.2被动失稳模式6.3节理岩体中的硐室开挖稳定性(命令流+FISH)6.3.1节理裂隙岩体的生成6.3.2初始应力状态控制6.3.3 开挖模拟 PFC3D与FLAC3D耦合模拟与分析

一套成熟的PFC 6.0 与 FLAC3D 耦合的数值模拟脚本。重点解决在数值模拟中复杂几何体建模难的问题，实现了通过 CAD直接导入转化为 PFC3D 的几何组（Geometry Group）及块体组（Clump Group），并在此基础上构建真三轴动力冲击耦合模型。

generate id 3 plane dip 90 dd 90 position -0.025 0 0 wall generate id 4 plane dip 90 dd 90 position 0.025 0 0 wall generate id 5 plane dip 90 dd 0 position 0 -0.025 0 wall generate id 6 plane

6.2盾构隧道掌子面稳定性(命令流+FISH)6.2.1主动失稳模式6.2.2被动失稳模式6.3节理岩体中的硐室开挖稳定性(命令流+FISH)6.3.1节理裂隙岩体的生成6.3.2初始应力状态控制6.3.3 开挖模拟 PFC3D与FLAC3D耦合模拟与分析

6.3节理岩体中的硐室开挖稳定性(命令流+FISH)6.3.1节理裂隙岩体的生成6.3.2初始应力状态控制6.3.3 开挖模拟 PFC3D与FLAC3D耦合模拟与分析 7 离散—连续域耦合模拟7.1离散—连续耦合模拟方法Ø 与FLAC3D中一维结构单元耦合Ø 与FLAC3D中二维壳结构单元或三维实体单元的面的耦合Ø 与FLAC3D